CN115808909A - Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment - Google Patents

Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment Download PDF

Info

Publication number
CN115808909A
CN115808909A CN202211474876.0A CN202211474876A CN115808909A CN 115808909 A CN115808909 A CN 115808909A CN 202211474876 A CN202211474876 A CN 202211474876A CN 115808909 A CN115808909 A CN 115808909A
Authority
CN
China
Prior art keywords
workpiece
heat treatment
lab
scheduling
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211474876.0A
Other languages
Chinese (zh)
Inventor
廖勇
彭乘风
李翔
钟宏扬
蒋纯志
雷大军
黄健全
谢光奇
段凌飞
张宏桥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xiangnan University
Original Assignee
Xiangnan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xiangnan University filed Critical Xiangnan University
Priority to CN202211474876.0A priority Critical patent/CN115808909A/en
Publication of CN115808909A publication Critical patent/CN115808909A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

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

Abstract

A dynamic batch scheduling method for on-time and energy-saving production of die heat treatment comprises the following steps: periodically acquiring information of a workpiece to be heat-treated, and taking part of the workpiece in the rolling period as an external cooperation workpiece according to the rolling period in the information of the workpiece to be heat-treated; and performing on-line dispatching on the rest workpieces on a hot workshop. The invention provides a proactive periodical rolling scheduling, which decomposes a dynamic scheduling problem into a series of deterministic subproblems according to a time axis in sequence, and formulates a workpiece external cooperation plan and a pre-scheduling scheme in a period, thereby balancing the production load of a short-term heat treatment workshop; therefore, a heat treatment production plan in a future period is obtained, a reliable outsourcing decision is made, and the proactive of the production plan is realized.

Description

Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment
Technical Field
The invention relates to the technical field of production scheduling, in particular to a dynamic batch scheduling method for on-time and energy-saving production of die heat treatment.
Background
Heat treatment is also a typical batch process, which is an important process in the manufacture of molds. The heat treatment production of the die is a key process which seriously influences the production progress and the manufacturing cost of the die. Most of the dies need heat treatment between rough machining and finish machining, which is an important process for ensuring the performance of the dies and directly influences the precision, the strength, the service life and the manufacturing cost of the dies. In actual production, the die manufacturing is often delayed in the heat treatment process, so the heat treatment often becomes the bottleneck of die manufacturing, and the punctuality and energy conservation of die production are greatly influenced. From the perspective of production management control (different from improving the heat treatment process or purchasing new equipment), it is of great practical significance to improve the operation efficiency of the mold production line by optimizing the batch scheduling scheme of mold heat treatment.
However, the mold heat treatment production is difficult to control in practice, and there are various problems: 1) There is a large variation between processing tasks. During the heat treatment, parts of different materials need to be heated to different critical temperatures, and the cooling time is different. This difference in the requirements of the heat treatment process results in an inability to batch at will between workpieces, and thus the mold heat treatment tasks are of the type typically incompatible with multiple workpiece families, and the tasks within the same workpiece family also differ in weight, lead time, and importance.
2) The arrival time of the processing task has high dynamic uncertainty. Different from the static batch scheduling problem, the processing tasks do not arrive at the same time at zero time, so that the processing tasks of the die heat treatment unit have dynamic arrival, therefore, the die heat treatment production control is also a dynamic batch scheduling problem, and not only the existing tasks but also the tasks to be arrived need to be considered at the scheduling time.
3) Enterprises often use an outsourcing processing mode to relieve the contradiction between production and demand. Because the order of the die enterprises often does not follow the uniform distribution, in order to ensure that the delivery date of the die is met as much as possible, when the actual production of a workshop exceeds the load, an outsourcing decision is often adopted, and part of workpieces which are possibly delayed are delivered to other enterprises for production, and the complexity of heat treatment management and control is further increased by the outsourcing decision.
4) The heat treatment production of the die is a key process for restricting delivery date and a high-energy consumption process, so that the scheduling optimization target is not only required to meet the optimization completion time, but also required to meet the energy-saving production target, and the heat treatment batch scheduling research of the die is required to be oriented to energy-saving and on-time manufacturing. It is because of these four main realistic features that the mold heat treatment production is often difficult to control effectively.
Disclosure of Invention
In view of the above-mentioned drawbacks, the present invention provides a method for dynamic batch scheduling of on-time and energy-saving production of mold heat treatment, so as to solve the problems of difficult on-schedule delivery, difficult external processing and difficult internal processing balance in the existing heat treatment workshop.
In order to achieve the purpose, the invention adopts the following technical scheme: a dynamic batch scheduling method for on-time and energy-saving production of die heat treatment comprises the following steps:
periodically acquiring information of a workpiece to be heat treated;
according to the rolling period in the information of the workpiece to be heat-treated, taking part of the workpiece in the rolling period as an external cooperation workpiece;
and performing on-line dispatching on the rest workpieces on a hot workshop.
Preferably, the step of obtaining the outside cooperation workpiece is as follows:
establishing a rolling scheduling mathematical model through the information of the workpiece to be heat-treated;
and solving the rolling scheduling mathematical model through an M-GEP algorithm to obtain a scheduling scheme of the outsourced workpieces.
Preferably, the specific steps of constructing the rolling scheduling mathematical model are as follows:
obtaining the production cost PI of the die workpiece j in the workpiece information to be heat-treated in the enterprise from the erp of the enterprise j Exterior cost PO j Calculating the total production cost in the rolling period and the extra cost for the outside cooperation;
setting the drag penalty of the workpiece according to the drag amount and the drag penalty coefficient of the workpiece;
predicting a stall penalty according to the total production cost, and establishing an objective function of a minimum stall penalty coefficient and an extra cost coefficient;
obtaining a scheduling scheme of the outside cooperation workpiece by analyzing the objective function;
the total production cost is as follows:
Figure BDA0003959509500000031
the additional cost in the scroll cycle is as follows:
Figure BDA0003959509500000032
where Nb is the number of workpieces in the rolling period, n o Number n of outer cooperation workpieces Em For number of emergency outsourcing workpieces, PI j For the cost of producing the die work j inside the enterprise, em _ PO j Is a processing time margin Tm j Less than outside preparation time T pre The outside cooperation cost of the die part, ord _ PO, is urgent j The outside cooperation cost is met for the workpiece j under the condition of meeting the outside cooperation preparation time;
the hangover penalty is as follows:
Figure BDA0003959509500000033
wherein J b For all the work sets contained in lot b,
Figure BDA0003959509500000034
The hold-off amount of the workpiece j in the batch b,
Figure BDA0003959509500000035
A drag penalty factor for the workpiece j in batch b;
the objective function is as follows:
Figure BDA0003959509500000036
Extra-Cost is the Extra Cost in the roll cycle, and Extra-Cost is the total production Cost, max-Tw k The maximum batch hold penalty, bn, is the number of batches processed by the furnace in a production cycle.
Preferably, the solving of the rolling scheduling mathematical model by the M-GEP algorithm is as follows:
step S1: initializing individuals and populations in the M-GEP, taking the scheduling schemes of the outsource plan and the hot workshop as the individuals in the M-GEP, and taking the scheduling schemes of the multiple outsource plans and the multiple hot workshops as the populations in the M-GEP;
wherein each individual in the M-GEP algorithm represents a scheduling plan, wherein decision information about the scheduling plans of the outside cooperation plan and the hot shop is stored in a chromosome, wherein the chromosome comprises two types of genes: the gene coding of the batch determines the batch rule of the workpieces, and the structural gene coding determines the loading sequence of the batch;
step S2: decoding chromosomes of an individual, wherein the decoding comprises batch gene decoding and structural gene decoding;
wherein the batching genes are decoded to calculate the priority of the workpieces, an outsource task set is selected and the rest task sets are batched;
decoding the structural gene into adjusting batch loading sequence;
and step S3: obtaining the reciprocal of each individual in the target function as the fitness value of the individual;
and step S4: judging whether the current evolution algebra reaches the maximum evolution algebra, if so, outputting an individual with the highest fitness value as a production scheduling scheme of the external cooperation workpiece; if not, evolving the algebra +1, selecting the individual with the highest response value to enter the next generation of population, and selecting N-2 individuals from the father generation of population to enter the next generation of population by adopting a roulette algorithm, wherein N is the number of the father generation of population;
step S5: and (3) carrying out genetic operation on individual chromosomes in the population, wherein the genetic operation comprises gene variation, insertion and chromosome cross recombination, and repeating the steps S2-S5.
Preferably, the chromosome is composed of a head segment and a tail segment, wherein the head segment is composed of a function and a terminator together or only of a function, and the tail segment is composed of a terminator;
wherein the function symbol is four operation symbols of { +, -, - ×,% }, wherein% is protective division, when the divisor is 0, the return value is 1;
the terminator is the attribute of the workpiece;
the head fragment and the tail fragment have the following relationship: lt = lh +1,lt is the tail fragment length and lh is the head fragment length.
Preferably, the online scheduling process of the workpieces on the hot shop is as follows:
a1, when a heat treatment machine of a heat treatment workshop is switched into an idle state or a new workpiece reaches the idle heat treatment machine, recording the current time as tnow, and establishing a prediction time window [ tnow, tnow + LAW ] of the current workpiece task i;
step A2: forming a plurality of workpiece tasks which reach the predicted time window at the moment of tnow into a task set J NT And selecting any workpiece task in a task set JNT to construct an LAB scheme, wherein the total mass of any workpiece task after combination is smaller than the rated load of the heat treatment machine, and each workpiece task correspondingly constructs an LAB squareCase;
step A3: judging whether the task time length of the next workpiece task I +1 exceeds the time length of the predicted time window, if not, updating the current time tnow by using the time of the next workpiece I +1 reaching the heat treatment machine, if so, selecting the optimal LAB scheme of a single workpiece in the workpiece task from the constructed LAB schemes, and recording the collection of the optimal LAB schemes of all the single workpieces as a set I *
Step A4: in set I * And selecting an optimal LAB scheme, recording the optimal LAB scheme as a final LAB scheme, using the final LAB scheme as a computer-on scheme of the heat treatment machine, judging whether the final LAB scheme is operated on the computer at the current moment, if so, immediately performing heat treatment processing on the workpiece of the final LAB scheme by the heat treatment machine, and if not, immediately processing the workpiece of the final LAB scheme after the heat treatment machine enters an idle state.
Preferably, the LAB protocol constructed in step A2 is specifically as follows:
Figure BDA0003959509500000061
wherein gamma is a target weight coefficient, D j Is a task set J NT The delivery date of the middle work piece j,
Figure BDA0003959509500000062
Is a task set J NT Minimum lead time of all workpieces,
Figure BDA0003959509500000063
Is a task set J NT Maximum lead time, w, in all workpieces j Is a task set J NT The weight of the middle workpiece j,
Figure BDA0003959509500000064
Is a task set J NT Average weight of the workpiece,
Figure BDA0003959509500000065
Is a task set J NT Middle workpieceMaximum difference in weight from the average weight.
Preferably, the specific steps of obtaining the optimal LAB scheme for a single workpiece in step A3 are as follows:
substituting the LAB scheme into an evaluation function, and acquiring a workpiece LAB scheme with the minimum evaluation function value in the LAB scheme as an optimal LAB scheme of a single workpiece;
wherein the merit function is as follows:
Figure BDA0003959509500000066
wherein gamma is the target weight coefficient, I WTI Weighted pull-off index, min (I) for LAB version I WTI ) Is the smallest weighted lag index among all LAB solutions; max (I) WTI ) Is the largest weighted lag index in all LAB solutions; i is EEI Energy consumption index, min (I) for LAB protocol I EEI ) Is the smallest energy consumption index in all LAB schemes; max (I) EEI ) Is the largest weighted pull finger in all LAB schemes.
Preferably, the formula for obtaining the final LAB solution is as follows:
Figure BDA0003959509500000067
wherein I is set I * Middle LAB protocol, J i* Representation set I * Set of all workpieces in, J i* Representing non-collections I * The collection of all the workpieces in the group,
Figure BDA0003959509500000068
is shown in the set J i* The pull-off penalty factor for the jth workpiece,
Figure BDA0003959509500000069
is shown in the set J i* The amount of the pull-out period of the jth workpiece,
Figure BDA0003959509500000071
is shown in the set J i* The amount of pull-out for the jth workpiece in (j'),
Figure BDA0003959509500000072
for the heat treatment completion time of the workpiece j',
Figure BDA0003959509500000073
represents a set J i* Heat treatment delivery date, n, of the jth workpiece i* Represents a set J i* Number of middle workpieces, n i* Representing non-aggregate J i* The number of workpieces;
wherein
Figure BDA0003959509500000074
Figure BDA0003959509500000075
Represents the average heat treatment time of the workpiece in j', W represents the payload of the heat treatment machine,
Figure BDA0003959509500000076
represents a set J i* The weight of the middle work piece j',
Figure BDA0003959509500000077
represents a set J i* The number of batches that can be made up,
Figure BDA0003959509500000078
represents a set J i* Average processing time of the medium workpiece task;
and acquiring an LAB scheme with the minimum M (i) value as a final LAB scheme.
One of the above technical solutions has the following advantages or beneficial effects: the invention provides a proactive periodical rolling scheduling, which decomposes a dynamic scheduling problem into a series of deterministic subproblems according to a time axis in sequence, and formulates a workpiece external cooperation plan and a pre-scheduling scheme in a period, thereby balancing the production load of a short-term heat treatment workshop; therefore, a heat treatment production plan in a future period is obtained, a reliable outsourcing decision is made, and the proactive of the production plan is realized.
Drawings
Fig. 1 is a schematic structural diagram of an embodiment of the present invention.
FIG. 2 is a flowchart of batch gene decoding according to one embodiment of the present invention.
FIG. 3 is a flow chart of gene decoding according to an embodiment of the present invention.
FIG. 4 is a schematic diagram of gene mutation according to an embodiment of the present invention.
FIG. 5 is a schematic cross-regrouping diagram of one embodiment of the present invention.
FIG. 6 is a LAB scheme build flow diagram of one embodiment of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
In the description of the embodiments of the present invention, the terms "first", "second", and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or to implicitly indicate the number of technical features indicated. Thus, features defined as "first", "second", may explicitly or implicitly include one or more of the described features. In the description of the embodiments of the present invention, "a plurality" means two or more unless specifically limited otherwise.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless otherwise specified. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
As shown in fig. 1 to 6, a method for dynamic batch scheduling for on-time and energy-saving production of mold heat treatment comprises the following steps:
periodically acquiring information of a workpiece to be thermally treated;
according to the rolling period in the information of the workpiece to be heat-treated, taking part of the workpiece in the rolling period as an external cooperation workpiece;
and performing on-line dispatching on the rest workpieces on a hot workshop.
Aiming at the management and control problem of heat treatment production, the invention makes two researches: firstly, selecting an external cooperation workpiece, and secondly, carrying out real-time on-machine decision in a heat treatment workshop. The researches in the two aspects are not independent, the selection of the outsourced workpiece has the effect of balancing the production load of a workshop, and the decision of real-time computer operation of the workpiece is directly influenced. In view of the uncertainty and multiple optimization goals of the mold thermal processing production environment, a hierarchical scheduling mechanism is proposed herein. Firstly, acquiring part of workpieces needing outsourcing, interfering the workpiece processing amount of a rolling period, and then carrying out online scheduling on the workpieces based on the rest, wherein the workpieces in an outsourcing plan are not added with a task set of online scheduling, and the rest of the workpieces realize real-time control through online scheduling, so that batching and on-machine decisions of a heat treatment workshop are completed.
The invention provides a proactive periodical rolling scheduling, which decomposes a dynamic scheduling problem into a series of deterministic subproblems according to a time axis in sequence, and formulates a workpiece external cooperation plan and a pre-scheduling scheme in a period, thereby balancing the production load of a short-term heat treatment workshop; therefore, a heat treatment production plan in a future period is obtained, a reliable outsourcing decision is made, and the proactive of the production plan is realized.
Preferably, the step of obtaining the outsourced workpiece is as follows:
establishing a rolling scheduling mathematical model through the information of the workpiece to be heat-treated;
and solving the rolling scheduling mathematical model through an M-GEP algorithm to obtain a scheduling scheme of the outsourced workpieces.
The heat treatment production process of the die is shown in figure 1, and most workpieces dynamically reach a heat treatment workshop from an upstream process, are subjected to heat treatment and then are transferred out to a downstream processing process. The heat treatment workshop mainly comprises two parts, namely a task pool without capacity limitation and a heat treatment furnace. After dynamically reaching a heat treatment workshop, the workpieces firstly enter a task pool, are batched and put on the machine in batches when the time of putting on the machine arrives, and enter a downstream procedure after heat treatment; and (3) carrying out outside processing on a small part of workpieces after the upstream process is finished without entering a heat treatment workshop, returning to the downstream process after the outside processing is finished, and finishing the production of the die together.
Preferably, the specific steps of constructing the rolling scheduling mathematical model are as follows:
obtaining the production cost PI of the die workpiece j in the workpiece information to be heat-treated in the enterprise from the erp of the enterprise j Exterior cost PO j Calculating the total production cost in the rolling period and the extra cost for the outside cooperation;
when the external cooperation workpiece is selected, the extra cost, the pull-in penalty and the internal manufacturing cost of the external cooperation are required to be compromised to be the production cost, an optimization model (and the objective function) is established with the aim of minimizing the production cost, and an external cooperation plan and a pre-scheduling scheme are obtained through calculation;
the rolling scheduling decision selected by the outsourced workpieces on the upper layer can act on the on-machine decision of the rest workpieces on the lower layer, and the fundamental reason is that the upper layer and the lower layer both meet the delivery date and reduce the delay penalty as the optimization target. The upper-layer rolling scheduling mechanism considers the extra production cost brought by the outside cooperation plan while reducing the possibility of generating the drag period of the future workpiece; and the lower-layer online scheduling reduces the production energy consumption as much as possible while considering the reduction of the work-piece stall penalty. The two have a common target, so the upper layer rolling decision has a direct effect on the lower layer real-time decision.
The choice of the external cooperation workpiece aims at balancing the production load of the workshop. If the production load of the workshop is too large, the delivery delay of the order of the enterprise is often caused, and the enterprise bears corresponding delay punishment (paying default money). The effective outsourcing strategy is a method for reducing the production load of a workshop and also plays a main role in scheduling.
Setting the drag penalty of the workpiece according to the drag amount and the drag penalty coefficient of the workpiece;
with n o The production load of the workshop is reduced, so that the risk of order pull-off of the enterprise is reduced. But also increases Extra _ Cost, which causes the loss of production profit of enterprises. Therefore, the development of the outsourcing scheme is to balance the profit and the risk of order pull-off sufficiently.
Predicting a stall penalty according to the total production cost, and establishing an objective function of a minimum stall penalty coefficient and an extra cost coefficient;
obtaining a scheduling scheme of the outside cooperation workpiece by analyzing the objective function;
in actual production, a contradiction relation exists between the extra cost and the pull-out penalty, and the addition of the external cooperation workpiece can reduce the pull-out penalty on one hand and increase the extra cost on the other hand. The need for upper layer scheduling finds the optimum between the two, for which an objective function is established that minimizes the stall penalty factor and the extra cost factor.
The total production cost is as follows:
Figure BDA0003959509500000111
the additional cost in the scroll cycle is as follows:
Figure BDA0003959509500000112
where Nb is the number of workpieces in the rolling period, n o Number n of outer cooperation workpieces Em For number of emergency outsourcing workpieces, PI j For die work j in the enterpriseCost of internal production, em _ PO j Is a processing time margin Tm j Less than outside preparation time T pre The outside cooperation cost of the die part, ord _ PO, is urgent j The outside cooperation cost is met for the workpiece j under the condition of meeting the outside cooperation preparation time;
the hangover penalty is as follows:
Figure BDA0003959509500000113
wherein J b For all the work sets contained in lot b,
Figure BDA0003959509500000114
The hold-off amount of the workpiece j in the batch b,
Figure BDA0003959509500000115
A drag penalty factor for the workpiece j in batch b;
the objective function is as follows:
Figure BDA0003959509500000116
Extra-Cost is the Extra Cost in the roll cycle, and Extra-Cost is the total production Cost, max-Tw k The maximum batch hold penalty, bn, is the number of batches processed by the furnace in a production cycle.
Scheduling decision is carried out based on predicted data of the workshop production condition in the whole period, and low robustness is an important bottleneck faced by the decision mode. In order to improve the robustness of the decision, the optimistic arrival time and the pessimistic arrival time of the workpiece can be respectively calculated to obtain two objective functions, and two outsourcing schemes are formulated according to the two objective functions. Finally, two outsourced schemes are given to the management layer for decision.
Preferably, the solving of the rolling scheduling mathematical model by the M-GEP algorithm is as follows:
the traditional GEP algorithm flow is as follows:
step1: initializing individuals and populations;
step2: calculating an individual fitness value;
step3: and (5) iteration times are plus 1, and whether the maximum iteration times are reached is judged. If the optimal solution is reached, outputting the optimal solution, and finishing the algorithm; otherwise, continuing to execute Step4;
step4: selecting excellent individuals to form a next generation population according to the fitness value of the individuals;
step5: random genetic manipulation was performed on the selected superior individuals. The genetic manipulation comprises: chromosomal gene variation, insertion, chromosomal cross recombination. Returning to Step2.
And the M-GEP is a GEP algorithm which is used with the aim of minimizing the lag time cost coefficient and the additional cost coefficient.
Step S1: initializing individuals and populations in the M-GEP, using the scheduling schemes of the external cooperation plan and the hot workshop as the individuals in the M-GEP, and using the scheduling schemes of the multiple external cooperation plans and the multiple hot workshops as the populations in the M-GEP;
wherein each individual in the M-GEP algorithm represents a scheduling plan, wherein decision information about the scheduling plans of the outside cooperation plan and the hot shop is stored in a chromosome, wherein the chromosome comprises two types of genes: the gene coding of the batch determines the batch rule of the workpieces, and the structural gene coding determines the loading sequence of the batch;
preferably, the chromosome is composed of a head segment and a tail segment, wherein the head segment is composed of a function and a terminator together or only of a function, and the tail segment is composed of a terminator;
the function symbol is four operation symbols of { +, - ×,% }, wherein% is protective division, and when the divisor is 0, the return value is 1;
the terminator is the attribute of the workpiece;
for example: the terminator may be
Figure BDA0003959509500000131
Sequentially comprises the following steps: arrival time, processing time allowance, weight, delivery date,A lingering penalty coefficient, an internal production profit, and an external profit. Chromosomes require the use of multiple batching genes to batch workpieces from different workpiece families, respectively. After the batching is finished, initializing a batch loading scheme according to the batching finishing sequence, and then adjusting the batch loading sequence by the same structural gene, wherein the batching gene is shown as follows in one embodiment:
Figure BDA0003959509500000132
in the above-mentioned batch genes, the workpieces are from 3 workpiece groups in total, and thus the chromosome contains 3 batch genes. Wherein the length of the head segment is 4 and the length of the tail segment is 5. In the context of the present invention,
Figure BDA0003959509500000134
the head segment and the tail segment are separated for convenient understanding.
In addition, the structural genes in this example are shown below:
Figure BDA0003959509500000133
the head fragment and the tail fragment have the following relationship: lt = lh +1,lt is the tail fragment length and lh is the head fragment length.
Step S2: decoding chromosomes of an individual, wherein the decoding comprises batch gene decoding and structural gene decoding;
wherein the batching genes are decoded to calculate the priority of the workpieces, an outsource task set is selected and the rest task sets are batched;
decoding the structural gene into adjusting batch loading sequence;
as shown in fig. 2, the batch genes of the workpiece family are converted into K-expression trees one by one, and the K-expression trees are decoded into a priority calculation formula by a width-first traversal method, that is, in the order from top to bottom and from left to right. And substituting the attribute value of the workpiece into a priority calculation formula corresponding to the workpiece family to calculate the priority value of the workpiece. The workpieces in each workpiece group are arranged from large to small according to the priority, and the smaller the priority is, the lower the priority of the workpiece is, and sufficient machining time allowance is provided for carrying out external coordination.
As shown in FIG. 3, after all the workpieces are batched, the batching order is initialized according to the batching completion time. The lot completion time is defined as the time of arrival at the heat treatment shop of the latest arriving workpiece in the lot. Without loss of generality, the processing time margin of some batches is smaller although the batch completion time is later, so that the batch processing sequence can be adjusted through the structural genes, the production scheduling plan is optimized, and the outsourcing scheme is more robust. In addition, although the empty lot is allowed to exist in the structural gene, the empty lot can be actually regarded as a lot with a heat treatment time of 0, and therefore the calculated completion time of the heat treatment of the lot is the same as the actually calculated value.
And step S3: obtaining the reciprocal of each individual in the target function as the fitness value of the individual;
and step S4: judging whether the current evolution algebra reaches the maximum evolution algebra, if so, outputting an individual with the highest fitness value as a production scheduling scheme of the external cooperation workpiece; if not, evolving the algebra +1, selecting the individual with the highest response value to enter the next generation of population, and selecting N-2 individuals from the father generation of population to enter the next generation of population by adopting a roulette algorithm, wherein N is the number of the father generation of population;
step S5: and (3) carrying out genetic operation on individual chromosomes in the population, wherein the genetic operation comprises genetic variation, insertion and chromosome cross recombination, and repeating the steps S2-S5.
The genetic manipulation of the M-GEP specifically comprises gene variation, insertion and chromosome cross recombination;
and carrying out genetic operation on the group genes and the structural genes simultaneously by the M-GEP. The individual batch scheme, the outsource scheme and the batch computer-on scheme are adjusted through the parallel evolution of the two genes, so that the evolution capability of the M-GEP is greatly improved, and the convergence speed is accelerated.
1) Genetic variation:
as shown in fig. 4: in M-GEP, the group genes and the structural genes adopt variant genetic operators.
For the head fragment of the set of genes, the first element can only be mutated into a function, the other elements of the head can be mutated into a function or terminator, and the elements on the tail fragment can only be mutated into a terminator.
For the structural gene, the head segment can only be mutated into a functional character, and the tail segment does not do any operation.
2) Inserting strings:
in the invention, only the group genes adopt genetic operators. Since the head segment of the structural gene cannot have a terminator, the structural gene does not employ an insertion operator.
The insertion string has two operations: IS insertion string and RIS root insertion string.
IS insertion string: a gene fragment, which may be 2 or 3 or 4 in length, is selected and copied from the tail fragment of the bulk gene and inserted at any position after the first element of the head fragment. After the head segment is inserted into the gene segment, elements behind the insertion point are moved backwards, redundant segments of the head are deleted, and the length of the head is kept unchanged.
The RIS inserts the cluster: selecting and copying a gene segment from any position behind the first element of the head of the group gene, wherein the gene segment should contain a complete expression, inserting the gene segment in front of the first element of the head, deleting the redundant segment of the head, and keeping the length of the head unchanged. If no gene fragment containing the complete formula can be found, the procedure is not performed.
3) Chromosome cross recombination:
as shown in fig. 5: both the batch genes and the structural genes adopt chromosome single-point crossing genetic operators.
Chromosome cross recombination is the most destructive genetic operator. Its destructive power is reflected in the possibility of generating entirely new chromosome structures by simple crossover operations.
Preferably, the online scheduling process of the workpieces on the hot shop is as follows:
as shown in fig. 6: a1, when a heat treatment machine of a heat treatment workshop is switched into an idle state or a new workpiece reaches the idle heat treatment machine, recording the current time as tnow, and establishing a prediction time window [ tnow, tnow + LAW ] of the current workpiece task i (a plurality of workpieces need to perform heat treatment in one workpiece task);
step A2: forming a plurality of workpiece tasks which reach the predicted time window at the moment of tnow into a task set J NT Selecting any workpiece task in a task set JNT to construct an LAB scheme, wherein the total mass of any workpiece task after combination is smaller than the rated load of a heat treatment machine, and each workpiece task correspondingly constructs an LAB scheme;
step A3: judging whether the task time length of the next workpiece task I +1 exceeds the time length of the predicted time window, if not, updating the current time tnow by using the time of the next workpiece I +1 reaching the heat treatment machine, if so, selecting the optimal LAB scheme of a single workpiece in the workpiece task from the constructed LAB schemes, and recording the collection of the optimal LAB schemes of all the single workpieces as a set I *
Step A4: in set I * And selecting an optimal LAB scheme, recording the optimal LAB scheme as a final LAB scheme, using the final LAB scheme as a computer-on scheme of the heat treatment machine, judging whether the final LAB scheme is operated on the computer at the current moment, if so, immediately performing heat treatment processing on the workpiece of the final LAB scheme by the heat treatment machine, and if not, immediately processing the workpiece of the final LAB scheme after the heat treatment machine enters an idle state.
Preferably, the LAB protocol constructed in step A2 is specifically as follows:
Figure BDA0003959509500000161
where γ is the target weight coefficient, D j Is a task set J NT The delivery date of the middle work piece j,
Figure BDA0003959509500000162
Is a task set J NT Minimum lead time of all workpieces,
Figure BDA0003959509500000163
Is a task set J NT Maximum lead time, w, in all workpieces j Is a task set J NT The weight of the middle workpiece j,
Figure BDA0003959509500000164
Is a task set J NT Average weight of the workpiece,
Figure BDA0003959509500000165
Is a task set J NT The maximum difference between the average weight and the median workpiece weight.
Preferably, the specific steps of obtaining the optimal LAB scheme for a single workpiece in step A3 are as follows:
substituting the LAB scheme into an evaluation function, and acquiring a workpiece LAB scheme with the minimum evaluation function value in the LAB scheme as an optimal LAB scheme of a single workpiece;
wherein the merit function is as follows:
Figure BDA0003959509500000171
wherein I WTI Weighted pull-off index, min (I) for LAB version I WTI ) Is the smallest weighted lag index among all LAB solutions; max (I) WTI ) Is the largest weighted lag index in all LAB solutions; i is EEI Energy consumption index, min (I) for LAB protocol I EEI ) Is the smallest energy consumption index in all LAB schemes; max (I) EEI ) Is the largest weighted pull finger in all LAB schemes.
Preferably, the formula for obtaining the final LAB solution is as follows:
Figure BDA0003959509500000172
wherein I is set I * Middle LAB protocol, J i* Representation set I * All work set in, J' i* Representing non-collections I * The collection of all the workpieces in the group,
Figure BDA0003959509500000173
is shown in the set J i* The pull-off penalty factor for the jth workpiece,
Figure BDA0003959509500000174
is shown in the set J i* The amount of the pull-out period of the jth workpiece,
Figure BDA0003959509500000175
is represented by the aggregate J' i* The amount of pull-out for the jth workpiece in (j'),
Figure BDA0003959509500000176
for the time to complete the heat treatment of the workpiece j',
Figure BDA0003959509500000177
represents a set J' i* Heat treatment delivery date, n, of the jth workpiece i* Represents a set J i* Number of workpieces, n' i* Denotes non-collective J' i* The number of workpieces;
wherein
Figure BDA0003959509500000178
Figure BDA0003959509500000179
Represents the average heat treatment time of the workpiece in j', W represents the payload of the heat treatment machine,
Figure BDA00039595095000001710
represents a set J' i* The weight of the middle work piece j',
Figure BDA00039595095000001711
represents a set J' i* The number of batches that can be made up,
Figure BDA00039595095000001712
represents a set J' i* Average processing time of the medium workpiece task;
and acquiring an LAB scheme with the minimum M (i) value as a final LAB scheme.
In the description of the present specification, reference to the description of "one embodiment," "some embodiments," "an illustrative embodiment," "an example," "a specific example," or "some examples" or the like means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and not to be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made in the above embodiments by those of ordinary skill in the art within the scope of the present invention.

Claims (9)

1. A dynamic batch scheduling method for on-time and energy-saving production of die heat treatment is characterized by comprising the following steps:
periodically acquiring information of a workpiece to be thermally treated;
according to the rolling period in the information of the workpiece to be heat-treated, taking part of the workpiece in the rolling period as an external cooperation workpiece;
and performing on-line dispatching on the rest workpieces on a hot workshop.
2. The method for dynamic batch scheduling for on-time and energy-saving production of die heat treatment according to claim 1, wherein the step of obtaining the outsourced workpiece is as follows:
establishing a rolling scheduling mathematical model through the information of the workpiece to be heat-treated;
and solving the rolling scheduling mathematical model through an M-GEP algorithm to obtain a scheduling scheme of the outsourced workpieces.
3. The method for dynamically scheduling the on-time and energy-saving production of the heat treatment of the mold according to claim 2, wherein the specific steps for constructing the mathematical model of the rolling scheduling are as follows:
obtaining the production cost PI of the die workpiece j in the workpiece information to be heat-treated in the enterprise from the erp of the enterprise j Exterior cost PO j Calculating the total production cost in the rolling period and the extra cost for the outside cooperation;
setting the drag penalty of the workpiece according to the drag amount and the drag penalty coefficient of the workpiece;
predicting a stall penalty according to the total production cost, and establishing an objective function of a minimum stall penalty coefficient and an extra cost coefficient;
obtaining a scheduling scheme of the outside cooperation workpiece by analyzing the target function;
the total production cost is as follows:
Figure FDA0003959509490000011
the additional cost in the scroll cycle is as follows:
Figure FDA0003959509490000021
where Nb is the number of workpieces in the rolling period, n o Number n of outer cooperation workpieces Em For number of emergency outsourcing workpieces, PI j Cost of production of the die work j inside the enterprise, em _ PO j Tm as a margin of processing time j Less than outside preparation time T pre The outside cooperation cost of the die part, ord _ PO, is urgent j The outside cooperation cost is met for the workpiece j under the condition of meeting the outside cooperation preparation time;
the hangover penalty is as follows:
Figure FDA0003959509490000022
wherein J b For all the work sets contained in lot b,
Figure FDA0003959509490000024
The hold-off amount of the workpiece j in the batch b,
Figure FDA0003959509490000025
A drag penalty factor for the workpiece j in batch b;
the objective function is as follows:
Figure FDA0003959509490000023
Extra-Cost is the Extra Cost in the roll cycle, and Extra-Cost is the total production Cost, max-Tw k The maximum batch hold penalty, bn, is the number of batches processed by the furnace in a production cycle.
4. The method for dynamically scheduling the batch for the on-time and energy-saving production of the die heat treatment according to claim 2, wherein the mathematical model of the rolling scheduling is solved by an M-GEP algorithm as follows:
step S1: initializing individuals and populations in the M-GEP, using the scheduling schemes of the external cooperation plan and the hot workshop as the individuals in the M-GEP, and using the scheduling schemes of the multiple external cooperation plans and the multiple hot workshops as the populations in the M-GEP;
wherein each individual in the M-GEP algorithm represents a scheduling plan, wherein decision information about the scheduling plans of the outside cooperation plan and the hot shop is stored in a chromosome, wherein the chromosome comprises two types of genes: the gene coding of the batch determines the batch rule of the workpieces, and the structural gene coding determines the loading sequence of the batch;
step S2: decoding chromosomes of an individual, wherein the decoding comprises batch gene decoding and structural gene decoding;
wherein the batching genes are decoded to calculate the priority of the workpieces, an outsource task set is selected and the rest task sets are batched;
decoding the structural gene into adjusting batch loading sequence;
and step S3: obtaining the reciprocal of each individual in the target function as the fitness value of the individual;
and step S4: judging whether the current evolution algebra reaches the maximum evolution algebra, if so, outputting an individual with the highest fitness value as a production scheduling scheme of the external cooperation workpiece; if not, evolving the algebra +1, selecting the individual with the highest response value to enter the next generation of population, and selecting N-2 individuals from the father generation of population to enter the next generation of population by adopting a roulette algorithm, wherein N is the number of the father generation of population;
step S5: and (3) carrying out genetic operation on individual chromosomes in the population, wherein the genetic operation comprises genetic variation, insertion and chromosome cross recombination, and repeating the steps S2-S5.
5. The dynamic batch scheduling method for the on-time and energy-saving production of the mold-oriented heat treatment, according to claim 4, wherein the chromosome is composed of a head segment and a tail segment, wherein the head segment is composed of a function symbol and a terminator together or only of a function symbol, and the tail segment is composed of a terminator;
the function symbol is four operation symbols of { +, - ×,% }, wherein% is protective division, and when the divisor is 0, the return value is 1;
the terminator is the attribute of the workpiece;
the head fragment and the tail fragment have the following relationship: lt = lh +1,lt is the tail fragment length and lh is the head fragment length.
6. The method for dynamically scheduling the on-time and energy-saving production of the heat treatment of the die as claimed in claim 1, wherein the on-line scheduling process of the workpieces on the hot workshop is as follows:
a1, when a heat treatment machine of a heat treatment workshop is switched into an idle state or a new workpiece reaches the idle heat treatment machine, recording the current time as tnow, and establishing a prediction time window [ tnow, tnow + LAW ] of the current workpiece task i;
step A2: forming a plurality of workpiece tasks which reach the predicted time window at the moment of tnow into a task set J NT Selecting any workpiece task in a task set JNT to construct an LAB scheme, wherein the total mass of any workpiece task after combination is smaller than the rated load of the heat treatment machine, and each workpiece task correspondingly constructs an LAB scheme;
step A3: judging whether the task time length of the next workpiece task I +1 exceeds the time length of the predicted time window, if not, updating the current time tnow by using the time of the next workpiece I +1 reaching the heat treatment machine, if so, selecting the optimal LAB scheme of a single workpiece in the workpiece task from the constructed LAB schemes, and recording the collection of the optimal LAB schemes of all the single workpieces as a set I *
Step A4: in set I * And selecting an optimal LAB scheme, recording the optimal LAB scheme as a final LAB scheme, using the final LAB scheme as a computer-on scheme of the heat treatment machine, judging whether the final LAB scheme is operated on the computer at the current moment, if so, immediately performing heat treatment processing on the workpiece of the final LAB scheme by the heat treatment machine, and if not, immediately processing the workpiece of the final LAB scheme after the heat treatment machine enters an idle state.
7. The method for dynamic batch scheduling for on-time and energy-saving production of mold heat treatment according to claim 6, wherein the LAB scheme constructed in step A2 is specifically as follows:
Figure FDA0003959509490000041
where γ is the target weight coefficient, D j As a set of tasks J NT The delivery date of the middle work piece j,
Figure FDA0003959509490000042
Is a task set J NT Minimum lead time of all workpieces,
Figure FDA0003959509490000043
Is a task set J NT Maximum lead time, w, in all workpieces j Is a task set J NT The weight of the middle workpiece j,
Figure FDA0003959509490000051
Is a task set J NT Average weight of workpiece,
Figure FDA0003959509490000052
Is a task set J NT The maximum difference between the average weight and the median workpiece weight.
8. The method for dynamically scheduling the on-time and energy-saving production of the die heat treatment according to claim 6, wherein the specific steps for obtaining the optimal LAB scheme of a single workpiece in the step A3 are as follows:
substituting the LAB scheme into an evaluation function, and acquiring a workpiece LAB scheme with the minimum evaluation function value in the LAB scheme as an optimal LAB scheme of a single workpiece;
wherein the merit function is as follows:
Figure FDA0003959509490000053
wherein I WTI Weighted pull-off index, min (I) for LAB version I WTI ) Is the smallest weighted lag index among all LAB solutions; max (I) WTI ) Is the largest weighted lag index in all LAB solutions; i is EEI Energy consumption index, min (I) for LAB protocol I EEI ) Is the minimum energy consumption index in all LAB schemes; max (I) EEI ) Is the largest weighted pull finger among all LAB schemes.
9. The method for dynamic batch scheduling for on-time and energy-saving production of mold heat treatment according to claim 6, wherein the formula for obtaining the final LAB scheme is as follows:
Figure FDA0003959509490000054
wherein I is set I * Middle LAB protocol, J i* Representation set I * Set of all workpieces in, J i* Represents a non-aggregate I * The collection of all the workpieces in the group,
Figure FDA0003959509490000055
is shown in the set J i* The pull-off penalty factor for the jth workpiece,
Figure FDA0003959509490000056
is shown in the set J i* The amount of the pull-out period of the jth workpiece,
Figure FDA0003959509490000057
is shown in the set J i* The amount of pull-out for the jth workpiece in (j'),
Figure FDA0003959509490000067
for the time to complete the heat treatment of the workpiece j',
Figure FDA0003959509490000061
represents a set J i* Heat treatment delivery date, n, of the jth workpiece i* Represents a set J i* Number of middle workpieces, n i* Representing non-aggregate J i* The number of middle workpieces;
wherein
Figure FDA0003959509490000062
Figure FDA0003959509490000063
Represents the average heat treatment time of the workpiece in j', W represents the payload of the heat treatment machine,
Figure FDA0003959509490000064
represents a set J i* The weight of the middle work piece j',
Figure FDA0003959509490000065
represents a set J i* The number of batches that can be made up,
Figure FDA0003959509490000066
represents a set J i* Average processing time of the medium workpiece task;
and acquiring an LAB scheme with the minimum M (i) value as a final LAB scheme.
CN202211474876.0A 2022-11-23 2022-11-23 Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment Pending CN115808909A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211474876.0A CN115808909A (en) 2022-11-23 2022-11-23 Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211474876.0A CN115808909A (en) 2022-11-23 2022-11-23 Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment

Publications (1)

Publication Number Publication Date
CN115808909A true CN115808909A (en) 2023-03-17

Family

ID=85483917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211474876.0A Pending CN115808909A (en) 2022-11-23 2022-11-23 Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment

Country Status (1)

Country Link
CN (1) CN115808909A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116402424A (en) * 2023-04-07 2023-07-07 广州拓威天海国际物流有限公司 Logistics configuration optimization adjustment method and system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106339774A (en) * 2016-08-18 2017-01-18 广东工业大学 Mould thermal treatment workshop dynamic batch scheduling method
CN112859761A (en) * 2020-12-06 2021-05-28 北京工业大学 Distributed forging flow shop energy-saving scheduling method considering centralized heat treatment
CN112947345A (en) * 2021-03-09 2021-06-11 河海大学 Deterministic reentrant sensor workshop-oriented dynamic batch scheduling intelligent optimization method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106339774A (en) * 2016-08-18 2017-01-18 广东工业大学 Mould thermal treatment workshop dynamic batch scheduling method
CN112859761A (en) * 2020-12-06 2021-05-28 北京工业大学 Distributed forging flow shop energy-saving scheduling method considering centralized heat treatment
CN112947345A (en) * 2021-03-09 2021-06-11 河海大学 Deterministic reentrant sensor workshop-oriented dynamic batch scheduling intelligent optimization method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116402424A (en) * 2023-04-07 2023-07-07 广州拓威天海国际物流有限公司 Logistics configuration optimization adjustment method and system
CN116402424B (en) * 2023-04-07 2024-01-23 广州拓威天海国际物流有限公司 Logistics configuration optimization adjustment method and system

Similar Documents

Publication Publication Date Title
WO2022000924A1 (en) Double-resource die job shop scheduling optimization method based on ammas-ga nested algorithm
CN111507641B (en) Batch processing equipment scheduling method and device
Sang et al. Intelligent factory many-objective distributed flexible job shop collaborative scheduling method
Zolfaghari et al. Comparative study of simulated annealing, genetic algorithms and tabu search for solving binary and comprehensive machine-grouping problems
CN112051825A (en) Multi-target production scheduling method considering employee operation capacity in automobile trial-manufacturing workshop
CN112668901B (en) Steel mill production scheduling method and system considering energy consumption
CN115808909A (en) Dynamic batch scheduling method for on-time and energy-saving production of die heat treatment
CN113822525B (en) Flexible job shop multi-target scheduling method and system based on improved genetic algorithm
CN113805545A (en) Flexible flow shop combined scheduling rule generation method considering batch processing
CN109858780A (en) A kind of Steelmaking-Continuous Casting Production Scheduling optimization method
CN110263970B (en) Hybrid discrete artificial bee colony algorithm for solving automatic rail trolley feeding scheduling problem
CN113139720B (en) Casting shop production scheduling optimization method with learning effect
Tian et al. Robust optimization of the continuous annealing process based on a novel Multi-Objective Dragonfly Algorithm
Tian et al. Dynamic operation optimization based on improved dynamic multi-objective dragonfly algorithm in continuous annealing process.
CN110728404A (en) Aluminum alloy part surface integrity prediction and optimization system
CN114021934A (en) Method for solving workshop energy-saving scheduling problem based on improved SPEA2
CN111665799A (en) Time constraint type parallel machine energy-saving scheduling method based on collaborative algorithm
CN114298567B (en) Method and system for scheduling casting schedule and dynamic decision-making of casting time of continuous casting machine
Sun et al. Scheduling of steel-making and continuous casting system using the surrogate subgradient algorithm for Lagrangian relaxation
CN114237166A (en) Method for solving multi-rotating-speed energy-saving scheduling problem based on improved SPEA2 algorithm
CN114879618A (en) Knowledge-driven workshop intelligent scheduling method
CN111445050B (en) Furnace assembling method for producing special alloy aluminum ingot
CN110688758B (en) Forging resistance furnace production optimization method based on SPEA2 algorithm
Hu A study of resource planning for precast production
CN113780871A (en) Multi-target low-carbon flexible job shop scheduling method

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