CN110909930B - Goods position distribution method of mobile goods shelf storage system for refrigeration house - Google Patents
Goods position distribution method of mobile goods shelf storage system for refrigeration house Download PDFInfo
- Publication number
- CN110909930B CN110909930B CN201911138483.0A CN201911138483A CN110909930B CN 110909930 B CN110909930 B CN 110909930B CN 201911138483 A CN201911138483 A CN 201911138483A CN 110909930 B CN110909930 B CN 110909930B
- Authority
- CN
- China
- Prior art keywords
- goods
- picking
- objective function
- value
- item
- 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
Links
- 238000003860 storage Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000009826 distribution Methods 0.000 title claims abstract description 21
- 238000005057 refrigeration Methods 0.000 title claims abstract description 14
- 241000196324 Embryophyta Species 0.000 claims abstract description 55
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 29
- 238000005457 optimization Methods 0.000 claims abstract description 20
- 238000009792 diffusion process Methods 0.000 claims abstract description 11
- 230000005484 gravity Effects 0.000 claims abstract description 6
- 230000002068 genetic effect Effects 0.000 claims abstract description 5
- 238000011156 evaluation Methods 0.000 claims description 21
- 238000003780 insertion Methods 0.000 claims description 15
- 230000037431 insertion Effects 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001174 ascending effect Effects 0.000 claims description 6
- 240000008791 Antiaris toxicaria Species 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000004083 survival effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000000151 deposition Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
- G06Q10/087—Inventory or stock management, e.g. order filling, procurement or balancing against orders
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Game Theory and Decision Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Educational Administration (AREA)
- Evolutionary Computation (AREA)
- Finance (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Accounting & Taxation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Warehouses Or Storage Devices (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A goods position distribution method of a mobile goods shelf storage system facing a refrigeration house distributes items with strong correlation to the same sorting roadway, reduces the possibility of opening the sorting roadway for many times, takes the similarity coefficient of order items as the basis of the correlation size, comprehensively considers the item sorting frequency and the goods shelf gravity center, establishes a multi-target goods position distribution optimization model, then solves by adopting an improved invasive weed algorithm to obtain the optimal storage position of goods, generates part of initial populations by adopting a greedy algorithm, then sets a reasonable space diffusion operator, and finally introduces the evolution reversal operation of a genetic algorithm. The invention has stronger overall search and local search capability, obvious optimization effect and effectively improved warehouse picking efficiency and shelf stability.
Description
Technical Field
The invention belongs to the field of storage management, and particularly relates to a goods position distribution method of a mobile goods shelf storage system for a refrigeration house.
Background
In recent years, with the rapid development of the cold-chain logistics industry, the refrigeration storage is concerned by more and more logistics enterprises. The problems of warehouse energy consumption, investment cost and efficiency are always pain points in the refrigeration house, so that the selection of a storage system with compact storage and taking space and timely service time becomes a new direction for the development of the refrigeration house. As a new compact storage system, the movable shelf storage system only needs to leave a sorting roadway for the storage trolley to operate, the storage trolley is moved to distribute out of the roadway through the shelf, and then enters the roadway to complete the storage and discharge of goods through the storage trolley, and the storage system is simple in structure, high in space utilization rate and low in cost, and is widely applied to various large cold storages at home and abroad. The goods position distribution problem is the key problem of the mobile goods shelf storage system facing to the refrigeration house, and the problem that how to optimize the operation efficiency of the mobile goods shelf storage system through a reasonable goods position strategy and how to optimize the stability of goods shelves becomes urgent to solve is directly related to whether the warehouse can run efficiently and stably.
Commonly used cargo space allocation strategies include: positioning storage, random storage, nearby position storage, full turnover rate storage and classified storage. However, in the prior art, the distribution of the goods space of the mobile shelf warehouse system facing the cold storage is less researched, and the specific storage position of each item in the warehouse system is not researched temporarily. Aiming at the logistics characteristics of various storage varieties, high timeliness requirement, high cost, complex technical requirement and the like of the logistics of the refrigeration house, a reasonable goods location distribution strategy is adopted, the order response speed of the refrigeration house can be improved, the cost of the refrigeration house is reduced, and the stability of a goods shelf is improved.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a goods position distribution method of a mobile goods shelf storage system facing a cold storage, so as to optimize a storage mode of the storage and reduce the management cost of the storage.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a goods position distribution method of a mobile goods shelf storage system facing a cold storage comprises the following steps:
calculating an objective function f for improving the sorting efficiency1,
The movement of the storage trolley from the I/O position to the cargo space (x, y, z) is calculated, assuming that the pallet needs to be moved, the horizontal running time is:
wherein y represents the coordinate of the cargo space (x, y, z) along the y-direction, w represents the cargo space width, vyRepresenting the speed of the access trolley in the y direction; t is trIndicating rack movement to open picking laneRoad time, k denotes the goods to be picked at the picking position (x, y, z) in the picking lane k, d denotes the depth of the goods grid, l denotes the width of the picking lane, vxThe speed of the storage trolley in the x direction is shown, as can be seen from fig. 2, the goods on the first row of shelves need to be picked from the lane 1, the goods on the second row and the third row of shelves are both picked from the lane 2, … …, and so on;
and then calculating the running time of the access trolley moving from the I/O position to the cargo space (x, y, z) in the vertical direction as follows:
wherein z represents the coordinate of the cargo space (x, y, z) along the z direction, h is the cargo grid height, vzThe speed of the access trolley in the vertical direction is shown;
continuing to calculate the return of the access cart from the cargo space (x, y, z) to the I/O location, the return time on the horizontal plane is:
the run time in the vertical direction on the return of the recalculation is:
further, since the movement of the storage trolley in the horizontal plane and the movement in the vertical direction are simultaneous, the time taken to pick the goods is the maximum value of the running times in the two directions, and the time t taken to pick the goods at the goods location (x, y, z) is calculatedxyzComprises the following steps:
wherein, tpThe time for loading and unloading goods for the storage trolley;
the time taken to pick the item at the cargo space (x, y, z) is then:
txyz·pxyz (6)
wherein p isxyzIndicating the picking frequency of items deposited at the cargo space (x, y, z).
Further, the main objective of improving the picking efficiency is to minimize the total picking time of the goods, and the objective function expression is as follows:
wherein f is1Representing the total picking time of the goods;
calculating an objective function f for improving shelf stability2,
Wherein f is2Indicates the height of the center of gravity of the goods shelf as a whole, mxyzRepresenting the weight of an item stored at the cargo space (x, y, z);
calculating and improving item correlation objective function f of same picking roadway3,
The calculation formula of the similarity coefficient proposed by Russel and Rao is as follows:
wherein a represents the order quantity containing item i and item j simultaneously; b represents the order quantity containing only item i; c represents the order quantity containing only item j; d represents the order quantity that neither item i nor item j contains;
the method has the advantages that items with strong correlation are stored on the shelf of the same sorting roadway, so that the times of opening the sorting roadway can be effectively reduced, furthermore, the reciprocal of the sum of the similarity coefficients of the items of the same sorting roadway converted from an objective function is as small as possible, and the expression of the objective function is as follows:
wherein f is3Expressing the reciprocal of the sum of item similarity coefficients of the same picking lane, K is the number of the picking lane, i and j represent item numbers, g represents the number of the items, rikItem i is stored on the rack of the sorting lane k when the value 1 indicates that item i is not stored, otherwise rik=0,rjkThe same process is carried out;
further constraints of the model are as follows:
wherein x is more than or equal to 1 and less than or equal to a represents the limitation of the number of rows of the goods shelves; y is more than or equal to 1 and less than or equal to b, which represents the limitation of the number of rows of the goods shelves; z is more than or equal to 1 and less than or equal to c represents the limitation of the number of the shelf layers;indicating that each item can only occupy one cargo space;
and 2, constructing an evaluation function.
Processing the three objective functions by an ideal point method to construct an evaluation function, and firstly, finding each objective function fiOptimum value of fi *And taking it as an ideal point;
based on the distance between each objective function value and the ideal point, an evaluation function of each objective function is constructed:
Fi=(fi-fi *)2 (12)
wherein, FiAn evaluation function value representing an ith objective function; f. ofiA function value representing an ith objective function; f. ofi *An optimal value representing the ith objective function;
further, on the basis of the above formula, a weight coefficient lambda is introducediThe sum of which is 1, the multiple target is optimizedTransformation of the chemo-function into the merit-function:
wherein f is a multi-objective optimization evaluation function.
Wherein f represents the objective function value of the current weed, fmaxAnd fminRepresenting the maximum and minimum values of the objective function representing the weeds in the current population respectively,represents rounding down;
wherein, σ represents a standard deviation value corresponding to the current algebra; itermaxRepresenting the maximum number of iterations; iter represents the current algebra; sigmainitRepresents an initial value of standard deviation; sigmafinalRepresents the final value of the standard deviation;
Further, the process of step 4 is as follows:
step 4.1, firstly, the time t required for picking the goods at the goods position (x, y, z) is calculatedxyzCreating a picking time matrix from which the time t required to pick a good at any location can be knownxyz;
4.2, using a greedy algorithm to generate 40% of initial weed populations, randomly generating the rest 60% of the weed populations, and using two greedy strategies in the initial weeds generated by the greedy algorithm, wherein the first method is to preferentially select a goods space with short picking time; the second is to preferentially select a cargo space with a small number of layers.
Still further, the process of step 4.2 is as follows:
step 4.2.1, preferentially selecting the goods location method with short picking time: taking out the serial number of each goods position and the corresponding picking time to create a matrix, wherein the first column is the goods position serial number, the second column is the corresponding picking time, then the goods position serial number and the picking time of each row are still corresponding, then the goods position serial number and the picking time are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one;
step 4.2.2, preferentially selecting the method for the goods location with the small layer number: taking out each goods position number and the number of layers where the goods position number is located to create a matrix, wherein the first column is the goods position number, the second column is the corresponding layer, then the goods position number and the number of layers of each row are still corresponding, then the goods position number and the number of layers are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one.
Further, the process of step 7 is as follows:
and 7.1, adopting a space diffusion operator based on sliding insertion, firstly randomly generating two random positions j and k, then taking the cargo space numbers from the position j to the position k as a queue, taking the number of the position j as the head of the queue, taking the number of the position k as the tail of the queue, then sequentially deleting the numbers from the tail of the queue, and inserting the numbers into the head of the queue until d times of insertion are performed. A schematic view of the sliding insertion is shown in fig. 5.
And 7.2, the dispersion process of the seeds follows normal distribution with the mean value of 0 and the standard deviation of sigma, and the size of sigma determines the search range of the seeds. Inspired by this, the number of execution times d of sliding insertion is made to obey N (0, σ)2) The absolute value of the random number alpha is taken and then is rounded upAs the number of execution times d of the slide insertion.
The invention has the following beneficial effects: aiming at the problem that the goods location distribution of a mobile goods shelf storage system facing a refrigeration house is less researched, particularly, the specific storage location of each item is not considered in the previous research, the invention considers the problem that the time for opening a picking roadway of a mobile heavy goods shelf is too long, proposes to distribute the items with strong correlation to the same picking roadway, reduces the possibility of opening the picking roadway for multiple times, takes the similarity coefficient of order items as the basis of the correlation size, comprehensively considers the picking frequency and the gravity center of the goods shelf, establishes a multi-target goods location distribution optimization model, then solves the optimal storage location of the goods by adopting an improved invasive weed algorithm, generates part of initial population by adopting a greedy algorithm, then sets a reasonable space diffusion operator, finally introduces the evolution operation of a genetic algorithm, has strong global search and local search capabilities, and has obvious optimization effect, effectively improve the warehouse picking efficiency and the shelf stability.
Drawings
FIG. 1 is a schematic diagram of a mobile rack warehousing system, wherein 1 is a mobile rack, 2 is a mobile track, 3 is an I/O, and 4 is an access cart.
FIG. 2 is a top view of the mobile rack storage system.
Figure 3 is a flow chart of the IIWO algorithm.
Fig. 4 is a diagram of an example of encoding.
Fig. 5 is a drawing showing a slide insertion.
FIG. 6 is a graph of simulation results for improving picking efficiency.
FIG. 7 is a graph of simulation results for improving shelf stability.
Fig. 8 is a diagram of a simulation result of improving item relevance of the same picking lane.
FIG. 9 is a diagram of a simulation result of the multi-objective evaluation function.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1 to 9, a goods location allocation method of a mobile goods shelf storage system facing a cold storage includes the following steps:
calculating an objective function f for improving the sorting efficiency1,
The movement of the storage trolley from the I/O position to the cargo space (x, y, z) is calculated, assuming that the pallet needs to be moved, the horizontal running time is:
wherein y represents the coordinate of the cargo space (x, y, z) along the y-direction, w represents the cargo space width, vyRepresenting the speed of the access trolley in the y direction; t is trIndicating the time of the shelf moving to open the picking roadway, k indicating the goods to be picked at the goods position (x, y, z) on the picking roadway k, d indicating the depth of the goods grid, l indicating the width of the picking roadway, vxThe speed of the storage trolley in the x direction is shown, as can be seen from fig. 2, the goods on the first row of shelves need to be picked from the lane 1, the goods on the second row and the third row of shelves are both picked from the lane 2, … …, and so on;
and then calculating the running time of the access trolley moving from the I/O position to the cargo space (x, y, z) in the vertical direction as follows:
wherein z represents the coordinate of the cargo space (x, y, z) along the z direction, h is the cargo grid height, vzThe speed of the access trolley in the vertical direction is shown;
continuing to calculate the return of the access cart from the cargo space (x, y, z) to the I/O location, the return time on the horizontal plane is:
the run time in the vertical direction on the return of the recalculation is:
further, since the movement of the storage trolley in the horizontal plane and the movement in the vertical direction are simultaneous, the time taken for picking the goods is the maximum value of the running time in the two directions, and the time t taken for picking the goods at the goods location (x, y, z) is calculatedxyzComprises the following steps:
wherein, tpThe time for loading and unloading goods for the storage trolley.
The time taken to pick the item at the cargo space (x, y, z) is then:
txyz·pxyz (6)
wherein p isxyzIndicating the picking frequency of items deposited at the cargo space (x, y, z).
Further, the main objective of improving the picking efficiency is to minimize the total picking time of the goods, and the objective function expression is as follows:
wherein f is1Representing the total pick time of the goods.
Calculating an objective function f for improving shelf stability2,
Wherein f is2Indicates the height of the center of gravity of the goods shelf as a whole, mxyzIndicating storage in goodsThe weight of the item at bit (x, y, z);
calculating and improving item correlation objective function f of same picking roadway3,
The calculation formula of the similarity coefficient proposed by Russel and Rao is as follows:
wherein a represents the order quantity containing item i and item j simultaneously; b represents the order quantity containing only item i; c represents the order quantity containing only item j; d represents the order quantity that neither item i nor item j contains.
The method has the advantages that items with strong correlation are stored on the shelf of the same sorting roadway, so that the times of opening the sorting roadway can be effectively reduced, furthermore, the reciprocal of the sum of the similarity coefficients of the items of the same sorting roadway converted from an objective function is as small as possible, and the expression of the objective function is as follows:
wherein f is3Expressing the reciprocal of the sum of item similarity coefficients of the same picking lane, K is the number of the picking lane, i and j represent item numbers, g represents the number of the items, rikItem i is stored on the rack of the sorting lane k when the value 1 indicates that item i is not stored, otherwise rik=0,rjkThe same is true.
Further constraints of the model are as follows:
wherein x is more than or equal to 1 and less than or equal to a represents the limitation of the number of rows of the goods shelves; y is more than or equal to 1 and less than or equal to b, which represents the limitation of the number of rows of the goods shelves; z is more than or equal to 1 and less than or equal to c represents the limitation of the number of the shelf layers;representing each itemOnly one cargo space is occupied.
And 2, constructing an evaluation function.
The invention selects an ideal point method to process the three objective functions and constructs an evaluation function. First, each objective function f is foundiOptimum value of fi *And takes it as an ideal point.
Further, based on the distance between each objective function value and the ideal point, an evaluation function of each objective function is constructed:
Fi=(fi-fi *)2 (12)
wherein, FiAn evaluation function value representing an ith objective function; f. ofiA function value representing an ith objective function; f. ofi *Representing the optimal value of the ith objective function.
Further, on the basis of the above formula, a weight coefficient lambda is introducediAnd if the sum is 1, converting the multi-objective optimization function into an evaluation function:
wherein f is a multi-objective optimization evaluation function.
Step 4.1, firstly, the time t required for picking the goods at the goods position (x, y, z) is calculatedxyzCreating a picking time matrix from which the time t required to pick a good at any location can be knownxyz;
And 4.2, generating 40% of initial weed population by using a greedy algorithm, and randomly generating the remaining 60% of weed population. Two greedy strategies are adopted in initial weeds generated by a greedy algorithm, wherein the first method is to preferentially select a goods location with short picking time; the second is to preferentially select a cargo space with a small number of layers.
Step 4.2.1, preferentially selecting the goods location method with short picking time: taking out the serial number of each goods position and the corresponding picking time to create a matrix, wherein the first column is the goods position serial number, the second column is the corresponding picking time, then the goods position serial number and the picking time of each row are still corresponding, then the goods position serial number and the picking time are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one;
step 4.2.2, preferentially selecting the method for the goods location with the small layer number: taking out each goods position number and the number of layers where the goods position number is located to create a matrix, wherein the first column is the goods position number, the second column is the corresponding layer, then the goods position number and the number of layers of each row are still corresponding, then the goods position number and the number of layers are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one;
Wherein f represents the objective function value of the current weed, fmaxAnd fminRepresenting the maximum and minimum values of the objective function representing the weeds in the current population respectively,indicating a rounding down.
wherein, σ represents a standard deviation value corresponding to the current algebra; itermaxRepresenting the maximum number of iterations; iter represents the current algebra; sigmainitRepresents an initial value of standard deviation; sigmafinalThe final value of standard deviation is indicated.
And 7.1, adopting a space diffusion operator based on sliding insertion, firstly randomly generating two random positions j and k, then taking the cargo space numbers from the position j to the position k as a queue, taking the number of the position j as the head of the queue, taking the number of the position k as the tail of the queue, then sequentially deleting the numbers from the tail of the queue, and inserting the numbers into the head of the queue until d times of insertion are performed. A schematic view of the sliding insertion is shown in fig. 5.
And 7.2, the dispersion process of the seeds follows normal distribution with the mean value of 0 and the standard deviation of sigma, and the size of sigma determines the search range of the seeds. Inspired by this, the number of execution times d of sliding insertion is made to obey N (0, σ)2) The absolute value of the random number alpha is taken and then is rounded upAs a slideNumber of times d of execution of live insertion.
The actual refrigeration house movable shelf warehousing system of a certain enterprise is used as a research object, and is programmed and simulated in matlab. The basic simulation parameters of the mobile shelf warehousing system researched by the invention are shown in table 1.
Table 1 shows the simulation basic parameters of the mobile shelf warehousing system
Simulation parameters | Value taking |
Width w of the grid | 1.3m |
Height h of the grid | 1.4m |
Depth d of the grid | 1.1m |
Total number of |
6 |
Total number of |
10 |
Number of |
4 |
X-direction speed v of storage trolley on horizontal planex | 2m/s |
Y-direction speed v of storage trolley on horizontal planey | 2m/s |
Speed v of storage trolley in vertical directionz | 1m/s |
The time for loading and unloading goods of the access trolley is tp | 5s |
Width l of picking roadway | 4.3m |
Speed v of movement of the goods shelfr | 4m/min |
The time for opening the picking roadway by moving the goods shelf is tr | 64.5s |
TABLE 1
In the optimization task of the distribution of certain goods space in the warehouse, 200 goods items are required to be stored in the warehouse, the warehousing system has 240 goods spaces, each goods space can only store one goods item, the sorting frequency and the weight of each warehousing goods item are known, and the similarity coefficient between every two goods items is calculated by the formula (9) based on the historical orders. In order to enable the warehouse to operate efficiently and stably, the improved invasive weed algorithm is adopted to carry out simulation analysis on the optimization of the allocation of the goods space of the warehouse.
IIWO algorithm parameter settings are shown in table 2. Table 2 shows IIWO algorithm parameters.
TABLE 2
The sorting efficiency simulation result is improved:
fig. 6 is an iterative curve of the modified invasive weed algorithm, with the ideal point of objective function 1 being 582.43, which is 0.7% better than 586.5 at the beginning.
Improving the simulation result of the center of gravity of the goods shelf:
fig. 7 is an iterative curve of the improved invasive weed algorithm, the ideal point of the objective function 2 is 2.59, and the optimization effect is 14.5% compared with 3.03 at the beginning.
And (3) improving the correlation simulation result of the same picking roadway item:
fig. 8 is an iterative curve of the modified invasive weed algorithm, the ideal point of the objective function 3 is 0.11, and the optimization effect is 17.3% compared with the initial 0.133.
And (3) multi-target evaluation function simulation results:
the enterprise focuses on improving shelf stability and item relevance of the same picking roadway by comprehensively considering three objective functions, and the importance of improving the picking efficiency is the lowest because the shelf moving and opening the picking roadway is long, and the improvement of the picking efficiency depends on the item relevance of the same picking roadway to a great extent, namely (lambda)1,λ2,λ3) And (5) taking values of (0.2, 0.4 and 0.4), and substituting the values into a multi-objective optimization evaluation function to carry out simulation solution.
FIG. 9 is an iterative curve of the improved invasive weed algorithm, the optimal value of the multi-objective optimization evaluation function is 0.24, and the optimization effect is 87.4% compared with the initial 1.91.
The simulation results of the four objective functions can be verified, and the method provided by the invention can better solve the problem of goods allocation optimization.
Claims (4)
1. A goods position distribution method of a mobile goods shelf warehousing system facing a refrigeration house is characterized by comprising the following steps:
step 1, establishing a multi-target goods space allocation optimization model for the purpose of improving the picking efficiency, improving the shelf stability and improving the relevance of the items of the same picking roadway;
calculating an objective function f for improving the sorting efficiency1,
The movement of the storage trolley from the I/O position to the cargo space (x, y, z) is calculated, assuming that the pallet needs to be moved, the horizontal running time is:
wherein y represents the coordinate of the cargo space (x, y, z) in the y-direction, w represents the cargo compartment width, vyRepresenting the speed of the access trolley in the y direction; t is trIndicating the time of the shelf moving to open the picking roadway, k indicating the goods to be picked at the goods position (x, y, z) on the picking roadway k, d indicating the depth of the goods grid, l indicating the width of the picking roadway, vxRepresenting the speed of the access trolley in the x direction;
and then calculating the running time of the access trolley moving from the I/O position to the cargo space (x, y, z) in the vertical direction as follows:
wherein z represents the coordinate of the cargo space (x, y, z) along the z direction, h is the cargo grid height, vzThe speed of the access trolley in the vertical direction is shown;
continuing to calculate the return of the access cart from the cargo space (x, y, z) to the I/O location, the return time on the horizontal plane is:
the run time in the vertical direction on the return of the recalculation is:
further, since the movement of the storage trolley in the horizontal plane and the movement in the vertical direction are simultaneous, the time taken to pick the goods is the maximum value of the running times in the two directions, and the time t taken to pick the goods at the goods location (x, y, z) is calculatedxyzComprises the following steps:
wherein, tpThe time for loading and unloading goods for the storage trolley;
the time taken to pick the item at the cargo space (x, y, z) is then:
txyz·pxyz (6)
wherein p isxyzRepresenting a picking frequency of items deposited at the cargo space (x, y, z);
further, the goal of improving the picking efficiency is to minimize the total picking time of the goods, and the expression of the objective function is as follows:
wherein f is1Representing the total picking time of the goods;
calculating an objective function f for improving shelf stability2,
Wherein f is2Indicates the height of the center of gravity of the goods shelf as a whole, mxyzRepresenting the weight of an item stored at the cargo space (x, y, z);
calculating and improving item correlation objective function f of same picking roadway3,
The calculation formula of the similarity coefficient proposed by Russel and Rao is as follows:
wherein A represents the order quantity simultaneously containing item i and item j; b represents the order quantity containing only item i; c represents the order quantity containing only item j; d represents the order quantity that neither item i nor item j contains;
the method has the advantages that items with strong correlation are stored on the shelf of the same sorting roadway, so that the times of opening the sorting roadway can be effectively reduced, furthermore, the reciprocal of the sum of the similarity coefficients of the items of the same sorting roadway converted from an objective function is as small as possible, and the expression of the objective function is as follows:
wherein f is3Expressing the reciprocal of the sum of item similarity coefficients of the same picking lane, K is the number of the picking lane, i and j represent item numbers, g represents the number of the items, rikItem i is stored on the rack of the sorting lane k when the value 1 indicates that item i is not stored, otherwise rik=0;rjkDepositing item j as 1On the shelf of lane k, otherwise rjk=0;
Further constraints of the model are as follows:
wherein x is more than or equal to 1 and less than or equal to a represents the limitation of the number of rows of the goods shelves; y is more than or equal to 1 and less than or equal to b, which represents the limitation of the number of rows of the goods shelves; z is more than or equal to 1 and less than or equal to c represents the limitation of the number of the shelf layers;indicating that each item can only occupy one cargo space;
step 2, establishing an evaluation function
Processing the three objective functions by an ideal point method to construct an evaluation function, and firstly, finding each objective function fiOptimum value of fi *And taking it as an ideal point;
based on the distance between each objective function value and the ideal point, an evaluation function of each objective function is constructed:
Fi=(fi-fi *)2 (12)
wherein, FiAn evaluation function value representing an ith objective function; f. ofiA function value representing an ith objective function; f. ofi *An optimal value representing the ith objective function;
further, on the basis of the above formula, a weight coefficient lambda is introducediAnd if the sum is 1, converting the multi-objective optimization function into an evaluation function:
wherein f is a multi-objective optimization evaluation function;
step 3, coding design: numbering the items and the goods positions simultaneously, adopting a natural number arrangement coding mode, wherein the coding length depends on the number of the items, when the number of the items is N, one code consists of N unrepeated natural numbers, and each natural number corresponds to a goods position number;
step 4, generating an initial population: initializing algorithm related parameters: initial number of weeds N0Maximum number of weed population NmaxMaximum number of iterations itermaxMaximum number of seeds S that can be produced per weedmaxAnd minimum value SminNonlinear modulation index n, initial value of standard deviation σ of weeds in spatial diffusioninitSum standard deviation final value σfinal;
Step 5, recording an objective function value of each weed, wherein the objective function value comprises f1,f2,f3And f, then recording the optimal weed individuals and the optimal solution;
step 6, a weed propagation stage: calculating the number of seeds generated by each weed according to the objective function value of each weed by using the following formula;
wherein f represents the objective function value of the current weed, fmaxAnd fminRepresenting the maximum and minimum values of the objective function representing the weeds in the current population respectively,represents rounding down;
step 7, a spatial diffusion stage: weed seeds are scattered around the parent weeds according to normal distribution with the mean value of 0 and the standard deviation of sigma, and the sigma can also be distributed from the initial value sigma as the iteration number is increasedinitDecrease to the final value σfinalSpecifically, the standard deviation is calculated as follows when a certain generation is reached:
wherein, σ represents a standard deviation value corresponding to the current algebra; itermaxRepresenting the maximum number of iterations; iter represents the current algebra; sigmainitRepresents an initial value of standard deviation; sigmafinalRepresents the final value of the standard deviation;
step 8, competition survival stage: after the parent individuals and the offspring individuals in the population are combined into a new population, in order to improve the local search capability of the algorithm, an evolution reversal operator of the genetic algorithm is introduced, two positions are randomly selected, the numbers of the two positions are interchanged, if the objective function value is reduced, the individuals are accepted, otherwise, the evolution reversal is invalid, then the new population is sorted according to the objective function value, excellent individuals with small objective function values are left, the disadvantaged individuals with large objective function values are eliminated, and the number of the individuals cannot exceed the maximum number N of the populationmax;
Step 9, completing the iteration once, and judging whether the maximum iteration number iter is reachedmaxAnd if so, outputting the optimal weed individuals, otherwise, returning to the step 6.
2. The method for allocating the cargo space of the mobile goods shelf warehouse system facing the cold storage according to claim 1, wherein the process of the step 4 is as follows:
step 4.1, firstly, the time t required for picking the goods at the goods position (x, y, z) is calculatedxyzCreating a picking time matrix from which the time t required to pick a good at any location can be knownxyz;
4.2, using a greedy algorithm to generate 40% of initial weed populations, randomly generating the rest 60% of the weed populations, and using two greedy strategies in the initial weeds generated by the greedy algorithm, wherein the first method is to preferentially select a goods space with short picking time; the second is to preferentially select a cargo space with a small number of layers.
3. The method for allocating the cargo space of the mobile goods shelf warehouse system facing the cold storage according to claim 2, wherein the process of the step 4.2 is as follows:
step 4.2.1, preferentially selecting the goods location method with short picking time: taking out the serial number of each goods position and the corresponding picking time to create a matrix, wherein the first column is the goods position serial number, the second column is the corresponding picking time, then the goods position serial number and the picking time of each row are still corresponding, then the goods position serial number and the picking time are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one;
step 4.2.2, preferentially selecting the method for the goods location with the small layer number: taking out each goods position number and the number of layers where the goods position number is located to create a matrix, wherein the first column is the goods position number, the second column is the corresponding layer, then the goods position number and the number of layers of each row are still corresponding, then the goods position number and the number of layers are arranged according to the ascending order of the second column, and at the moment, the generated partial population codes correspond to the first column one by one.
4. The goods location distribution method of the mobile goods shelf storage system facing the cold storage according to any one of claims 1 to 3, wherein the process of the step 7 is as follows:
step 7.1, adopting a space diffusion operator based on sliding insertion, firstly randomly generating two random positions j and k, then taking the goods space numbers from the position j to the position k as a queue, taking the number of the position j as the head of the queue, taking the number of the position k as the tail of the queue, then sequentially deleting the numbers from the tail of the queue, and inserting the numbers into the head of the queue until d times of insertion are executed;
step 7.2, the dispersion process of the seeds obeys normal distribution with the mean value of 0 and the standard deviation of sigma, the size of sigma determines the search range of the seeds, and the execution times d of the sliding insertion obeys N (0, sigma)2) The absolute value of the random number alpha is taken and then is rounded upAs the number of execution times d of the slide insertion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138483.0A CN110909930B (en) | 2019-11-20 | 2019-11-20 | Goods position distribution method of mobile goods shelf storage system for refrigeration house |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138483.0A CN110909930B (en) | 2019-11-20 | 2019-11-20 | Goods position distribution method of mobile goods shelf storage system for refrigeration house |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110909930A CN110909930A (en) | 2020-03-24 |
CN110909930B true CN110909930B (en) | 2022-05-03 |
Family
ID=69818154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911138483.0A Active CN110909930B (en) | 2019-11-20 | 2019-11-20 | Goods position distribution method of mobile goods shelf storage system for refrigeration house |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909930B (en) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111553063B (en) * | 2020-04-20 | 2022-03-08 | 广州地铁设计研究院股份有限公司 | Scheduling method for solving resource-limited project by invasive weed algorithm |
CN113554199B (en) * | 2020-04-23 | 2024-04-23 | 华晨宝马汽车有限公司 | Computer-implemented method and apparatus for optimizing inventory of materials |
CN111582787A (en) * | 2020-04-29 | 2020-08-25 | 上海久耶供应链管理有限公司 | Freezing warehouse scale type picking system and method thereof |
CN111709681B (en) * | 2020-06-03 | 2023-07-14 | 大连九州创智科技有限公司 | Discrete storage position selection method |
CN111860837B (en) * | 2020-07-20 | 2024-06-18 | 上海汽车集团股份有限公司 | Method and device for processing boxing problem and computer readable storage medium |
CN112100861B (en) * | 2020-09-22 | 2024-05-14 | 河南中烟工业有限责任公司 | Cigarette production material cargo space distribution method based on invasive weed optimization algorithm |
CN112700183B (en) * | 2020-12-08 | 2023-10-17 | 中电九天智能科技有限公司 | Automatic sorting and carrying method based on AGV |
CN112633729B (en) * | 2020-12-29 | 2022-06-10 | 杭州电子科技大学 | Multi-compartment material vehicle cargo space optimization method based on human factors and Epsilon greedy algorithm |
CN112950109B (en) * | 2021-01-28 | 2022-05-17 | 浙江大学 | Complex network-based associated article storage position optimization method |
CN113064392B (en) * | 2021-03-22 | 2023-09-08 | 聊城大学 | Discrete optimization method based on matrix workshop AGV scheduling |
JP7386485B2 (en) * | 2021-03-25 | 2023-11-27 | ダイキン工業株式会社 | Information processing equipment and programs |
CN113371383B (en) * | 2021-06-25 | 2023-01-10 | 深圳市库宝软件有限公司 | Goods shelf scheduling method, device, equipment, warehousing system and storage medium |
CN113371382B (en) * | 2021-06-25 | 2023-01-10 | 深圳市库宝软件有限公司 | Robot path planning method, device, equipment, warehousing system and storage medium |
CN113371381B (en) * | 2021-06-25 | 2022-12-13 | 深圳市库宝软件有限公司 | Shelf scheduling method, device, equipment, system, medium and program product |
CN113371380B (en) * | 2021-06-25 | 2022-11-22 | 深圳市库宝软件有限公司 | Path generation method, device, equipment, storage medium and program product |
CN113570025B (en) * | 2021-07-16 | 2022-11-18 | 东南大学 | E-commerce storage center goods space distribution method based on discrete particle swarm optimization |
CN114648272B (en) * | 2022-04-01 | 2023-07-21 | 上海聚货通电子商务有限公司 | Commodity layout adjustment method and system based on goods picking thermodynamic diagram |
CN115439069B (en) * | 2022-09-24 | 2023-08-01 | 北京融安特智能科技股份有限公司 | Intelligent inventory method, device, equipment and storage medium for unmanned archive warehouse |
CN116579721B (en) * | 2023-07-14 | 2023-09-19 | 中油管道物资装备有限公司 | Warehouse goods position optimization method and device, electronic equipment and storage medium |
CN118037187A (en) * | 2023-11-02 | 2024-05-14 | 北京极智嘉科技股份有限公司 | Warehouse management method, device, equipment and readable storage medium |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106779153B (en) * | 2016-11-15 | 2021-08-03 | 浙江工业大学 | Intelligent stereoscopic warehouse goods space allocation optimization method |
JP6323787B1 (en) * | 2017-04-10 | 2018-05-16 | 株式会社オープンロジ | Warehouse management server and warehouse management method |
CN107480922B (en) * | 2017-07-07 | 2021-03-16 | 西安建筑科技大学 | Method for establishing goods position distribution scheduling model under two-end type same-rail double-vehicle running mode |
CN107808215B (en) * | 2017-10-23 | 2021-11-19 | 南昌大学 | Goods allocation optimization method applied to Flying-V type non-traditional layout warehouse |
CN109597304B (en) * | 2018-11-30 | 2022-02-11 | 北京工业大学 | Intelligent partitioned storage method for mold library based on artificial bee colony algorithm |
-
2019
- 2019-11-20 CN CN201911138483.0A patent/CN110909930B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110909930A (en) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110909930B (en) | Goods position distribution method of mobile goods shelf storage system for refrigeration house | |
CN107480922B (en) | Method for establishing goods position distribution scheduling model under two-end type same-rail double-vehicle running mode | |
CN109886478B (en) | Goods space optimization method for finished wine automatic stereoscopic warehouse | |
CN108550007B (en) | Goods space optimization method and system for automatic stereoscopic warehouse of pharmaceutical enterprise | |
CN113222293B (en) | Intelligent stereoscopic warehouse optimal scheduling method | |
CN106779153B (en) | Intelligent stereoscopic warehouse goods space allocation optimization method | |
CN107808215B (en) | Goods allocation optimization method applied to Flying-V type non-traditional layout warehouse | |
CN109597304B (en) | Intelligent partitioned storage method for mold library based on artificial bee colony algorithm | |
CN103559396B (en) | Based on the automatic dispensary stock's allocation optimization method improving chaos particle cluster algorithm | |
CN114417696B (en) | Automatic stereoscopic warehouse cargo space distribution optimization method based on genetic algorithm | |
CN113222410B (en) | Method for establishing cargo space distribution model in bidirectional layout mode | |
CN109081030B (en) | Configuration optimization method of primary and secondary shuttle type intensive warehousing system | |
CN101968860A (en) | Order sorting method and system | |
CN107967586A (en) | A kind of power grid goods and materials storage optimization method | |
CN109583660B (en) | Method for realizing dynamic goods picking strategy | |
CN111815233B (en) | Goods position optimization method based on total logistics amount and energy consumption | |
CN111626516B (en) | Order ordering optimization method of double-deep four-way shuttle system considering cargo pouring strategy | |
CN105858043B (en) | The warehousing system Optimization Scheduling that a kind of lift is combined with shuttle | |
Mahajan et al. | A retrieval sequencing heuristic for miniload end-of-aisle automated storage/retrieval systems | |
CN112100861B (en) | Cigarette production material cargo space distribution method based on invasive weed optimization algorithm | |
CN112580852A (en) | Intensive automatic stereoscopic warehouse goods space optimization method for electric power materials | |
CN112561448B (en) | Multi-task stacker scheduling algorithm based on improved genetic algorithm | |
CN117408604A (en) | Large-scale stereoscopic warehouse cargo space distribution method based on intelligent calculation | |
CN116596440A (en) | Automatic stereoscopic warehouse-in and warehouse-out intelligent scheduling method | |
CN116342039A (en) | Optimizing method for goods distribution and sorting of stereoscopic warehouse |
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 |