<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JSS</journal-id><journal-title-group><journal-title>Open Journal of Social Sciences</journal-title></journal-title-group><issn pub-type="epub">2327-5952</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jss.2016.47008</article-id><article-id pub-id-type="publisher-id">JSS-69477</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Business&amp;Economics</subject><subject> Social Sciences&amp;Humanities</subject></subj-group></article-categories><title-group><article-title>
 
 
  Study of Optimization Assigned on Location Selection of an Automated Stereoscopic Warehouse Based on Genetic Algorithm
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tongjuan</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiaoxia</surname><given-names>Ma</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiaohui</surname><given-names>Zhan</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Beijing Wuzi University, Beijing, China</addr-line></aff><pub-date pub-type="epub"><day>18</day><month>07</month><year>2016</year></pub-date><volume>04</volume><issue>07</issue><fpage>52</fpage><lpage>58</lpage><history><date date-type="received"><day>13</day>	<month>May</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>18</month>	<year>July</year>	</date><date date-type="accepted"><day>25</day>	<month>July</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   In twenty-first century, automated stereoscopic warehouse has attracted many attentions of the enterprises because it has high working efficiency. The operation efficiency and management benefit of an automated stereoscopic warehouse are affected directly by the order picking efficient, it is an important symbol of service level about the automated stereoscopic warehouse. This paper from the perspective of the warehouse management to discuss the issue, the optimization assigned of location selection is regarded as the ultimate goal. A multi-objective mathematical model is established by setting “job efficiency” and “turnover rate” as objective functions, “warehouse space layout” as the constraint condition, and the genetic algorithm is used for coding this mathematical model which is applied to a pharmaceutical warehouse and solved by using Matlab software, result shows that the picking efficiency of normal temperature zone and the turnover of cargo area are have been greatly improved. 
 
</p></abstract><kwd-group><kwd>Automated Stereoscopic Warehouse</kwd><kwd> Optimization Assigned on Location Selection</kwd><kwd> Genetic Algorithm</kwd><kwd> The MATLAB Simulation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Comparing with the ordinary warehouse, the automated warehouse has some advantages, such as low error rate of operation, high space utilization, and so on. The key factor of automated warehouse favored by enterprises is the high working efficiency. There are many factors that influence operation efficiency of the automated warehouse, which can be roughly divided into: mechanical equipment and warehouse management. Optimization assigned on location selection, belongs to the warehouse management problem, which is the focus of this paper.</p></sec><sec id="s2"><title>2. Problem Description</title><sec id="s2_1"><title>2.1. Object of Study</title><p>Three key entities are selected as the research objects to analyze the problem of location optimization accurately, namely “cargo”, “cargo location” and “handling equipment”, the mathematical model is established to analyze the equilibrium relationship among cargo, cargo location and handling equipment, and finally we realize the goal of optimization assigned on location selection.</p></sec><sec id="s2_2"><title>2.2. Fundamental Assumption</title><p>To make the mathematical model standardized, on the basis of objective reality, to the best solution oriented, the model assumptions are as follows:</p><p>(1) The storage unit is a pallet or container, goods stored in a single pallet or container is the same kind;</p><p>(2) All the size of the pallet or container are the same, and match with the cargo box, (the box size is the same), each cargo box, just can accommodate a pallet (or container).</p><p>(3) Shelves display are rectangular, the cargo box is a cube shape;</p><p>(4) In-out warehouse’s frequency of all types of cargo in the system will not change too much in a short time;</p><p>(5) The initial position of a stacker is located on the position which shelves closed to the entrance and exit;</p><p>(6) In the horizontal and vertical operations, the stacker speed value is a constant;</p><p>(7) The stacker can only save or take a cargo unit at a time;</p><p>(8) The goods can be directly complete stock pick without putting into warehouse;</p></sec><sec id="s2_3"><title>2.3. Warehouse Layout Description</title><p>Automated stereoscopic warehouse has a total of A row of the shelf, each row has B column and C layer, one of the closest to the ground as the first layer, the nearest to entrance and exit as the first column, there is a stacker machine in each roadway, and a row of shelves on both sides of the roadway, a container (or pallet) on the shelf sides can be made accessing operation by a stacker.</p></sec></sec><sec id="s3"><title>3. Model Building</title><sec id="s3_1"><title>3.1. The Objective Function and Constraint Conditions</title><p>In the process of optimization assigned on location selection , this article selects two optimizing principles that are “nearest picking” principle based on in-out warehouse frequency and “first in first out” principle based on turnover rate, this paper establishes the specific optimization objectives which are “minimizing the total operation time” and “the total storage time of chosen goods is maximum”, combining with the constraint: warehouse space limitation, the “multi-objective mathematical model of optimization assigned on location selection” is built.</p><sec id="s3_1_1"><title>3.1.1. Objective Function S<sub>1</sub></title><p>The mathematical expression of objective function S<sub>1</sub>considering minimum operation time is as shown below:</p><disp-formula id="scirp.69477-formula83"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x4.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula84"><graphic  xlink:href="http://html.scirp.org/file/69477x5.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula85"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x6.png"  xlink:type="simple"/></disp-formula><p>S<sub>1</sub>: expresses the sum picking time of goods waiting for out the warehouse;</p><p>T<sub>xyz</sub>: expresses the outbound operating time of goods located at position (x, y, z);</p><p>q<sub>ixyz</sub>: 0 - 1 decision variables;</p><p>i: the ith kind goods to be taken in the outbound order;</p><p>n<sub>i</sub>: the ex-warehouse quantity of the ith kind goods;</p><p>N: the type number of goods in the outbound order;</p><p>(x, y, z): the position coordinates of goods waiting for picking;</p><p>(A, B, C): expresses the maximum number of row, column and layer of the shelves, respectively, those have a total of A row, B column, and C layer;</p><p>w: the length of the cube cargo box;</p><p>l: the roadway width;</p><p>V<sub>1</sub>: the running speed of conveyor belt;</p><p>V<sub>2</sub>: the mobile speed of a stacker in horizontal direction (Y direction);</p><p>V<sub>3</sub>: the mobile speed of a stacker invertical;</p></sec><sec id="s3_1_2"><title>3.1.2. The Objective Function S<sub>2</sub></title><p>The mathematical expression of the objective function S<sub>2</sub> based on the principle as follows:</p><disp-formula id="scirp.69477-formula86"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x7.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula87"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x8.png"  xlink:type="simple"/></disp-formula><p>S<sub>2</sub>: indicates the sum storage time of selected goods;</p><p>T: indicates current time point of picking operation;</p><p>T<sub>XYZ</sub>: indicates the storage time point of goods located at position (x, y, z).</p></sec><sec id="s3_1_3"><title>3.1.3. Constraint Condition</title><p>1 ≤ x ≤ A; 1 ≤ y ≤ L; 1 ≤ z ≤ C, x, y and Z are all positive integers.</p></sec></sec><sec id="s3_2"><title>3.2. Multi-Objective Optimization Model</title><p>The multi-objective optimization model is shown.</p><disp-formula id="scirp.69477-formula88"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x9.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula89"><graphic  xlink:href="http://html.scirp.org/file/69477x10.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4"><title>4. Design of Model Algorithm Based on Genetic Algorithm</title><p>The optimization assigned on location selection s a problem that which locations cargo should be taken out to complete stock pick, ensure the highest operation efficiency and meet the rule “first in first out”.</p><sec id="s4_1"><title>4.1. Chromosome Coding Design</title><p>Coding problem, which is how to code the solution of problem becomes chromosome [<xref ref-type="bibr" rid="scirp.69477-ref1">1</xref>]. This article encodes the chromosome by adopting the binary coding method. Specifically, the chromosome gene site refers to the selected outbound goods, if a certain location was selected out of the warehouse, the corresponding gene value is 1, otherwise, the value is 0. The binary coding is a conclusive presentation, lacking specific description of the cargo information; therefore, it needs to use auxiliary string to describe the cargo information.</p></sec><sec id="s4_2"><title>4.2. Fitness Function Setting and Genetic Operator Designing</title><sec id="s4_2_1"><title>4.2.1. Fitness Function Setting</title><p>In general, fitness function of genetic algorithm is the maximum problem, when it is not, it should be transformed at first [<xref ref-type="bibr" rid="scirp.69477-ref2">2</xref>].</p><p>The fitness function is designed to the value which the objective function S<sub>2</sub>add the reciprocal of the objective function value S1 plus 1. The fitness function F (x, y, z) is expressed as:</p><disp-formula id="scirp.69477-formula90"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x11.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula91"><graphic  xlink:href="http://html.scirp.org/file/69477x12.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4_3"><title>4.3. Genetic Operator Design</title><p>(1) Selection operator design</p><p>It is selected the operator selection in this paper by using “the roulette wheel selection method”, Assuming the population number is N, and the fitness function value of individual i is Fi, the corresponding selection probability is P<sub>i</sub>:</p><disp-formula id="scirp.69477-formula92"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x13.png"  xlink:type="simple"/></disp-formula><p>Each round will produce a uniform random number between [0, 1] as a pointer to determine selected individual. The individual can be made the crossover operation after being selected.</p><p>(2) Crossover operator design</p><p>This paper chooses the way of “single point crossover”, namely selecting two individuals from the group to make interlace operation after selection operation and randomly generate a crossing point which the two individuals interchange parts genetic code, thus forming two new individuals.</p><p>(3) Mutation operator design</p><p>Mutation operation is to change the chromosome genevalue of each individual in the current group according to initial setting of the mutation probability [<xref ref-type="bibr" rid="scirp.69477-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.69477-ref4">4</xref>]. Based on the binary coding, this paper makes the gene value of chromosome code string to inverse operation and achieves the effect of variation, thus forming a new individual.</p></sec><sec id="s4_4"><title>4.4. Design of Control Parameters and Algorithm Termination Conditions</title><p>In advance, some parameters are set in the algorithm run, “the control parameters”. In this article, the control parameters of the model are: population size M, termination evolution algebra T, crossover probability Pc, and mutation probability P<sub>m</sub>. Among them, the population size refers to the individuals number contained in the group, the range is from 20 to 100; the termination evolution algebra of genetic algorithm, the range is 100 - 500; crossover probability P<sub>C</sub>, the range is from 0.4 to 0.99; mutation probability P<sub>M</sub>, the range is from 0.0001 to 0.1.</p></sec></sec><sec id="s5"><title>5. Case Verification and Result Analysis</title><sec id="s5_1"><title>5.1. Case Background</title><p>The storage area of an automated warehouse of a pharmaceutical company in Daxing Beijing Economic and Technological Development Zone including normal temperature area, easy odor area, cold storage area, special medicine area and so on.</p></sec><sec id="s5_2"><title>5.2. Case Description</title><p>Two rows shelves area of the 1th roadway in normal temperature zone is selected as the analysis area, the region has 2*10*8 cargo locations.</p><p>The basic constant parameters of the handling equipment and shelves in the automated warehouse are as follow:</p><p>v<sub>1</sub> = 2 m/s, v<sub>2</sub> = 1 m/s, v<sub>3</sub> = 1 m/s.</p><p>A = 2, B = 10, C = 8, w = 1 m, l = 1.2 m.</p><p>The time how long it takes on putting in storage and other information about ten locations of drug 4 is seen in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>This article, firstly, solves assignment location model of one class cargo, gets the drugs chosen location assignment scheme and solves the rest class drugs waiting for out-warehouse by category, ultimately, acquiring picking operation list.</p></sec><sec id="s5_3"><title>5.3. Model Solving Process</title><sec id="s5_3_1"><title>5.3.1. Model Deformation</title><p>This article considers the two goals are equally important and give equal weight, changes the second objective function into the reciprocal. The transformed single objective model is shown below:</p><disp-formula id="scirp.69477-formula93"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x14.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula94"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69477x15.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69477-formula95"><graphic  xlink:href="http://html.scirp.org/file/69477x16.png"  xlink:type="simple"/></disp-formula><p>1 ≤ x ≤ 2, 1 ≤ y ≤ 10, 1 ≤ z ≤ 8, ∑q<sub>ixyz</sub> = n<sub>i</sub>.</p></sec><sec id="s5_3_2"><title>5.3.2. Solution Process</title><p>Step 1: constructing chromosomes by using a binary encoding method;</p><p>Step 2: setting control parameters: population size NIND=20; crossover probability pc=0.8; mutation probability pm = 0.1; maximum genetic algebra MAXGEN=200;</p><p>Step 3: generating the initial population: Gen=1, randomly generating initial group Pop (0), including 160 individuals;</p><p>Step 4: if Gen &gt; MAXGen, the algorithm meets the termination conditions, it stops, otherwise, it should be continue;</p><p>Step 5: calculating the objective function value according to the population of all chromosome information;</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Drug 4 store state information</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Drug type</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th></tr></thead><tr><td align="center" valign="middle" >Location</td><td align="center" valign="middle" >1 6 1</td><td align="center" valign="middle" >1 7 1</td><td align="center" valign="middle" >1 8 1</td><td align="center" valign="middle" >1 9 1</td><td align="center" valign="middle" >1 10 1</td></tr><tr><td align="center" valign="middle" >In-warehouse time</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-5</td><td align="center" valign="middle" >2015-3-5</td></tr><tr><td align="center" valign="middle" >Drug type</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Location</td><td align="center" valign="middle" >2 6 1</td><td align="center" valign="middle" >2 7 1</td><td align="center" valign="middle" >2 8 1</td><td align="center" valign="middle" >2 9 1</td><td align="center" valign="middle" >2 10 1</td></tr><tr><td align="center" valign="middle" >In-warehouse time</td><td align="center" valign="middle" >2015-2-28</td><td align="center" valign="middle" >2015-2-28</td><td align="center" valign="middle" >2015-3-10</td><td align="center" valign="middle" >2015-3-10</td><td align="center" valign="middle" >2015-3-10</td></tr></tbody></table></table-wrap><p>Step 6: calculating the fitness value of each individual;</p><p>Step 7: selecting individuals according to the fitness value size through the roulette method, then getting the next generation population by crossover and mutation operation;</p><p>Step 8: Gen = Gen + 1;</p><p>Step 9: go to Step 4.</p></sec></sec><sec id="s5_4"><title>5.4. Results Analysis</title><p>After running in the MATLAB software, the change curve of the objective function S can be derived, see <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>The curve change trend in the figure is showed that the objective function value is decreasing in the iteration process; the total operating time tends to be minimized, while the total storage time tends to be maximized. After 160 iterations, the fluctuation in the objective function value decreases and tends to be stable, and its value is S = 20.39. Therefore, the model is feasible.</p><p>After the end of the iteration, assignment result of the location selection is shown in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>It can be seen from the results of the table above, the selected operation location is located on the six cargo locations of column 6, 7, 8 in the first row and the second row.</p><p>From the above analysis, the assignment model of location selection can be solved by MATLAB software and obtain ideal location assignment strategy, meanwhile the model meets the requirements which the operation efficiency and turnover rate, so it has the feasibility.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The objective function change curve</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69477x17.png"/></fig><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> The distribution result would be</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Drug type</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th><th align="center" valign="middle" >4</th></tr></thead><tr><td align="center" valign="middle" >Drug location</td><td align="center" valign="middle" >1 6 1</td><td align="center" valign="middle" >1 7 1</td><td align="center" valign="middle" >1 8 1</td><td align="center" valign="middle" >1 9 1</td><td align="center" valign="middle" >1 10 1</td></tr><tr><td align="center" valign="middle" >In-warehouse time</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-2</td><td align="center" valign="middle" >2015-3-5</td><td align="center" valign="middle" >2015-3-5</td></tr><tr><td align="center" valign="middle" >Selection decision</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Drug type</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Drug location</td><td align="center" valign="middle" >2 6 1</td><td align="center" valign="middle" >2 7 1</td><td align="center" valign="middle" >2 8 1</td><td align="center" valign="middle" >2 9 1</td><td align="center" valign="middle" >2 10 1</td></tr><tr><td align="center" valign="middle" >In-warehouse time</td><td align="center" valign="middle" >2015-2-28</td><td align="center" valign="middle" >2015-2-28</td><td align="center" valign="middle" >2015-3-10</td><td align="center" valign="middle" >2015-3-10</td><td align="center" valign="middle" >2015-3-10</td></tr><tr><td align="center" valign="middle" >Selection decision</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td></tr></tbody></table></table-wrap></sec></sec><sec id="s6"><title>6. Conclusions</title><p>In his article, using the multi-objective genetic algorithm theory and MATLAB software tool.</p><p>Selecting the pharmaceutical automation stereoscopic warehouse as research object, solving the chosen location assignment problem by optimization and simulation analysis, we confirm the built model is feasible. The model and the algorithm compiler have scientificity and certain reference value for the study of others related picking optimization problem, meanwhile it can provide a good decision-making basis for managers. In addition, the different conditions hypothesis of the model simplification will have different effects on the solution. This issue has not been discussed in this paper, which will be the focus of further research.</p></sec><sec id="s7"><title>Acknowledgements</title><p>This work was supported by Beijing outstanding talent training project (2014000020124G093) and Beijing Intelligent Logistics System Collaborative Innovation Centre.</p></sec><sec id="s8"><title>Cite this paper</title><p>Tongjuan Liu,Xiaoxia Ma,Xiaohui Zhan, (2016) Study of Optimization Assigned on Location Selection of an Automated Stereoscopic Warehouse Based on Genetic Algorithm. Open Journal of Social Sciences,04,52-58. doi: 10.4236/jss.2016.47008</p></sec></body><back><ref-list><title>References</title><ref id="scirp.69477-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Liu, E.G., Wang, G.C. and Zhu, X.S. (2014) MATLAB and Mathematics Experiment. National Defence Industry Press, Beijing.</mixed-citation></ref><ref id="scirp.69477-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Pan, Q. (2013) Multi-Source Information Fusion Theory and Application. Tsinghua University Press, Beijing.</mixed-citation></ref><ref id="scirp.69477-ref3"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Huang</surname><given-names> B.L. </given-names></name>,<etal>et al</etal>. (<year>2011</year>)<article-title>Application of Adaptive Genetic Algorithm in Intelligent Generating Test Paper</article-title><source> Computer Engineering</source><volume> 37</volume>,<fpage> 161</fpage>-<lpage>163</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.69477-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Tao, X.M., Liu, Y., Fu, D.D., et al. (2011) Hybrid Clonal Selection Algorithm for Multi Objective Optimization. Com- puter Simulation, 28, 199-203.</mixed-citation></ref></ref-list></back></article>