Three-Dimensional Container Loading: A Simulated Annealing .

3y ago
20 Views
2 Downloads
913.80 KB
15 Pages
Last View : 16d ago
Last Download : 3m ago
Upload by : Elisha Lemon
Transcription

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.comThree-Dimensional Container Loading: A Simulated AnnealingApproach1Hanan Mostaghimi,2Bryan St. Amour and1#Walid Abdul-KaderFaculty of Engineering, University of Windsor, Windsor (Ontario), Canada12Department of Computer Science, University of Windsor, Windsor (Ontario), Canada#ORCID: 0000-0002-2926-5059to a compact, which results in a potential loss of the marketwhen compared to other competitor's vehicles (Hifi, 2002;Bortfeldt et al., 2003; and Wang et al., 2008).AbstractHigh utilization of cargo volume is an essential factor in thesuccess of modern enterprises in the market. Althoughmathematical models have been presented for container loadingproblems in the literature, there is still a lack of studies thatconsider practical constraints. In this paper, a Mixed IntegerLinear Programming is developed for the problem of packing asubset of rectangular boxes inside a container such that the totalvalue of the packed boxes is maximized while some realisticconstraints, such as vertical stability, are considered. Thepacking is orthogonal, and the boxes can be freely rotated intoany of the six orientations. Moreover, a sequence triple-basedsolution methodology is proposed, simulated annealing is usedas modeling technique, and the situation where some boxes arepreplaced in the container is investigated. These preplacedboxes represent potential obstacles. Numerical experiments areconducted for containers with and without obstacles. Theresults show that the simulated annealing approach issuccessful and can handle large number of packing instances.The basic form of the container loading problem is packing thebest subset of rectangular boxes (called cargo) into a largeobject (called container) to maximize the total value of theloaded boxes. The boxes should not be overlapped and shouldlie entirely in the container. According to the topologyintroduced by Wascher et al. (2007), the containers can beeither homogeneous or heterogeneous. If the boxes placed inthe given container are identical, it is called homogeneous;however, if various types of boxes are loaded, the container isconsidered as heterogeneous. Besides the non-overlappingconstraints, some other practical constraints should beconsidered in the real-world container loading problem, such ascargo vertical stability, preplaced boxes, and box rotation(Junqueira et al., 2012; Bortfeldt et al., 2012). However, notmany papers have considered these practical constraints in theirproposed models.The problem addressed in this research and as per the topologypresented by Dyckhoff (1990), belongs to 3/B/O/F (3: threedimensional, B/O: one object/bin and items selection, F: fewitems of different types), while Wascher et al. (2007) classify itas the three-dimensional single orthogonal knapsack problem.The given problem considers the packing of rectangular itemsinto a container to maximize the total value of the packed itemsby minimizing the amount of lost space. The value of boxes isassumed to be equal to their volume. The rotation of the boxesis considered as well. The multidimensional knapsack Problem(MKP) is a NP-hard optimization problem that can be shownby reduction from the one-dimensional packing problem(Egeblad and Pisinger, 2009). Although technologicalknowledge has been enhanced, solving real knapsack problemsis still a challenge. Due to NP-hardness of the packing problem,only heuristic methods, and a few exact algorithms have beenpresented.Keywords: knapsack, packing sequence, rotation, obstacles,simulated annealingINTRODUCTIONLogistics has recently played an important role in the successof modern enterprises. Packing boxes inside a container is anessential material handling activity in manufacturing andtransportation industries. It is also a key function for operatingsupply chain efficiently. The efficient use of transportationdevices, like containers and palettes, leads to significant costsaving. Moreover, high utilization of transportation devicesreduces the traffic of goods and protects natural resources.Therefore, optimal loading of a container decreases theshipping cost and increases the stability of the load. Containerloading problem has practical values, and it can be applied tovarious fields. Loading cars, trucks, trains, or ships can be alsoconsidered as a container loading or a three-dimensionalpacking problem. Furthermore, cargo volume is an importantfactor used by motor vehicle companies to market theirproducts as sub-compact, compact, midsize, or full size. Ifinaccurate, the cargo volume found may downgrade the vehiclefrom its real size. For example, a mid-size may be downgradedIn this paper, a mixed integer linear model and a simulatedannealing algorithm are developed to address a iderations, such as vertical stability and preplaced(obstacles) constraints, are tackled. These practical constraintsand box rotation contribute significantly to a study of a more1290

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.comproposed for packing a batch of objects. Goncalves et al. (2012)propose a multi-population biased random-key geneticalgorithm for the single container loading problem. Maximalspace representation is used to manage the container free space.The authors consider stability and orientation constraints;however, they do not develop a mathematical model for thegiven problem. Peng et al. (2009) present a hybrid simulatedannealing algorithm for three-dimensional container loadingproblem. Firstly, a heuristic algorithm is used for encodingfeasible packing solutions, and then the simulated annealingalgorithm is applied to search in the encoding neighborhood.Hongtao et al. (2012) address a three-dimensional singlecontainer loading problem by using a multi-stage search basedsimulated annealing algorithm.Wei et al. (2012) use areference length approach to address the three-dimensionalstrip packing problem. In another paper, Wei et al. (2015)address the problem of multiple container loading costminimization problem by using a new approach that combinescolumn generation technique with a goal driven search.realistic 3D knapsack problem as proposed in this researchwork. The aim of this proposed research is to provide managersand decision-makers with an adequate modeling tool that helpsmake shipping goods more efficient and eco-friendly as fewertrips can be made if higher container utilization is achieved.Also, auto-makers can market the class size of their vehiclesmore accurately.LITERATURE REVIEWThe focus of most of the container loading or three-dimensionalcutting and packing problems is on the rectangular bins. Feketeand Schepers (2004) propose a new method for obtaining lowerbound classes for higher-dimensional packing problem. Themajor objective of this paper is to define good criteria forremoving a candidate set of boxes. Dual feasible function is away to build conservative scales. The computational results aremainly limited to the two-dimensional packing problem. Hifi(2004) proposes a dynamic algorithm and an exact depth-firstsearch to solve the three-dimensional cutting problem.Orientation and guillotine constraints are considered. Optimalsolutions are obtained for a significant number of instances, butnot all of them. Althaus et al. (2007) consider the trunk packingproblem where the box dimensions are as per the SAE J1100standard. They propose two discretized methods. First, thespace to be packed is discretized. Then, an approximationapproach is considered using linear inequalities. The spacediscretization causes insufficient representation of the boxes.Additionally, the runtime of the enumerative algorithms isexponential.Other research works related to design automation focus onthree-dimensional placement of circuit elements by exploringthe layout of the integrated circuits (Obenaus and Szymansky,2003), and Cheng et al. (2005) address floor planning for 3-DVLSI Design. While this technique is less known to containerloading practices, it carries some similarities and it is moreefficient in terms of search time than other methods such aspartitioning placement.Models that provide information on optimal objective functionvalue and bounds help to assess the solution quality of heuristicalgorithms. Although modeling three-dimensional knapsackproblems considers practical constraints, it is still at itsbeginnings. Junqueira et al. (2012) present mixed integer linearprogramming models for the container loading problem.Vertical and horizontal stability of the cargo, as well as cargoload bearing strength, are considered in the proposed model.However, the models are only able to handle moderate sizedproblems. Table 1 compares some relevant papers and models,and shows their similarities and differences.Although considerable advancement has been made in thedevelopment of exact algorithms, heuristic algorithms still playan important role in solving three-dimensional knapsackproblems. Only heuristic methods can provide reasonablesolutions within an acceptable running time for instances ofreal-world sized problems. Pisinger (2002) develops a wallbuilding based heuristic. Both homogenous and heterogeneousinstances are considered. Moreover, several ranking rules arestudied to select the best layers’ depths. Bortfeldt et al. (2003)propose a parallel tabu search approach for a single containerloading problem and give little consideration to heterogeneousinstances. Wang et al. (2008) also present a heuristic methodfor a heterogeneous container loading problem and developeda dynamic space decomposition approach based on the tertiarytree structure. Egeblad and Pisinger (2009) propose a simulatedannealing based methodology for the two and threedimensional knapsack problems, and a three-dimensionalknapsack model is presented. The authors present an iterativeheuristic approach for the knapsack problem that is based onthe sequence triple representation. Also, Yamazaki et al. (2000)apply a variety of packing sequences including sequence-triplein their 3-D packing solution approach. To control the heuristicmethod, simulated annealing is used. However, rotation ofboxes and preplaced boxes (obstacles) are not considered in thethree-dimensional model and experiments. Wu et al. (2010)consider the three-dimensional bin packing problem withvariable bin heights. A mixed integer programming model isproposed, and they also present the case when more than onetype of bin is used. A genetic algorithm-based heuristic isPer the literature, not all papers consider box rotation since itincreases the search space significantly. Bin stabilityconstraints have likewise been just considered in a few papers.To the best of our knowledge, preplaced boxes (obstacles) havenot been studied in three-dimensional knapsack problems,although it is so essential for such problems since it is oftenrequired to place certain boxes in certain positions. Theseconstraints can also be used in the case of having a nonrectangular container. Therefore, it is important to study morepractical constraints in the knapsack problem. This proposedresearch work aims to contribute to the literature so that a 3Dknapsack problem can be tackled where box rotation isconsidered to help finding more practical packingconfigurations. Furthermore, preplaced boxes (bin with someobstacles) and vertical stability constraints are considered andaddressed. This is useful, especially when considering trunkloading for auto size classification as indicated earlier in theIntroduction section.1291

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.comthat the boxes can be freely rotated in six different orientations.However, it is possible to relax this constraint and fix a box ina specific orientation. The boxes need not to be packed inlayers, and the bottom of each box must be supported by thetop of other boxes or the bin floor. In addition, some boxeswhose left-bottom-behind (LBB) corner should be placed in aspecific position are considered as preplaced boxes orobstacles. The value of each box is equal to its volume. It isassumed that the dimensions of all boxes and the bin areintegers; thus, the placement is to be done in integer steps. LetC be a rectangular bin with width W, height H and depth D.The origin of the Cartesian coordinate system is located at theLBB corner of the container, wi, hi, and di are respectively, thewidth, height and depth of box type i. and (xi, yi, zi) representthe coordinates of the LBB corner of the box.Goncalveset al.(2012) Wu et al.(2010) Egeblad &Pisinger(2009)Just 2DWang et al.(2008) Bortfeldt etal. (2003) Pisinger(2002) ProposedResearch Mathematical ModelPreplaced Boxes Heuristic ApproachJunqueiraet al.(2012)Vertical StabilityPaperBox RotationTable 1: Models Parameters ComparisonAble tosolveLargeProblems A mixed integer programming formulation is presented for thegiven problem. Some real-world knapsack problem constraintsare considered in the model, and to the best of our knowledge,they have not been addressed previously. These constraints arevertical stability and preplaced boxes (obstacles). Since thethree-dimensional knapsack problem is NP-hard, it is difficultto solve. Additionally, the flexibility of the orientation of boxesincreases the search space significantly so that the difficulty offinding the optimal solution is increased as well. Some exactalgorithms and heuristic methods are proposed in the publishedliterature. As exact algorithms require more time to find asolution, heuristic approaches are more popular and can beeffective alternatives to finding an optimal or near optimalsolution. The proposed three-dimensional solutionmethodology is based on sequence triple representation, whichis defined below under Placement Algorithm. Simulatedannealing is used as the meta-heuristic method. As the numberof box types (or box dimensions and variety) is finite, the useof simulated annealing is favoured by its efficiency inneighborhood search.PROBLEM DEFINITIONMATHEMATICAL FORMULATIONIn this study, the three-dimensional knapsack problem isconsidered where there is one bin with fixed size and a set ofboxes, and each box has an associated size. The aim is to findan efficient solution methodology to pack rectangular boxes ina single bin so that the total value of the packed boxes ismaximized, or equivalently the empty spaces left areminimized. The boxes are assumed to be stronglyheterogeneous, which means that there is relatively manydifferent types of boxes and a small number of boxes for eachbox type (Wascher et al., 2007). Moreover, the packing isconsidered feasible if each box lies entirely in the bin and thepacked boxes do not overlap. The edges of all boxes must beparallel to the edges of the bin (orthogonal packing). The boxesare assumed to be of rectangular shape; however, the bin canbe considered either of rectangular or nonrectangular shape. Inthe case of having preplaced boxes (obstacles), the bin isassumed to be non-rectangular.A mixed-integer programming model of the 3D-knapsackproblem is formulated in this section. The mathematical modelis based on work done by Egeblad and Pisinger (2009) and Wuet al. (2010). Some modifications are made to their models andinclude new constraints addressing vertical stability andobstacles, which were not considered in any previouslypublished works. Constraints (1) – (4) are adapted fromEgeblad and Pisinger (2009); however, the authors did notconsider box orientation in their model. The binary positionvariables, which show the orientation of the boxes, areintegrated in constraints (5) – (17). This makes the model morecomprehensive. They are described below in this section.The following are the main assumptions considered for the mixinteger linear model:1.2.3.4.The boxes are strongly heterogeneous,The boxes must be located orthogonally,The boxes can freely rotate,The box and bin dimensions are assumed to be nonnegative integer,5. The value of a box is equal to its volume, andSome practical considerations that play an important role inmodeling more realistic knapsack problems, such as boxrotation and bin stability, are presented. The algorithm assumes1292

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.com 6. The X, Y, and Z axes of the bin are shown in the followingfigure.Parameters:(wi, hi, di): width, height, and depth of box iY(W, H, D): width, height, and depth of the bin(r, s, k): Left-Bottom Behind (LBB) coordinates of thepreplaced boxes (obstacles)XZ(a, b, c, d): binary orientation parameters of the preplaced boxesPi: value of box iFigure 1: The X, Y, and Z axes of the binNotationM: large numberThe model notation for the parameters and variables are asfollows:Pb: set of preplaced boxes Variables:(xi,yi,zi): Left-Bottom Behind (LBB) coordinates of box iXwi:Zwi:Yhi:Zdi:rij, lij:oij, uij:bij, fij:si:π‘¦π‘–π‘—π‘Ž :1if width of box i is parallel to the container’s X0otherwise1if width of box i is parallel to the container’s Z0otherwise1if height of box i is parallel to the container’s Y0otherwise1if depth of box i is parallel to the container’s Z0Otherwise1if box i is to the right (rij) or to left (lij) of box j0otherwise1if box i is over (oij) or under (uij) box j0otherwise1if box i is behind (bij) or in front (fij) of box j0otherwise1if box i is packed0otherwise1if π‘₯𝑗 π‘₯𝑖 (x coordinate of the LBB corner of box j is greater than or equal to x coordinate of the LBB corner ofbox i)π‘₯π‘–π‘—π‘Ž :0otherwise1if π‘₯𝑗 π‘₯𝑖′ (x coordinate of the LBB corner of box j is less than x coordinate of the Right-Bottom-Behind (RBB)0otherwisecorner of box i)1293

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.com𝑦𝑖𝑗𝑏 :π‘₯𝑖𝑗𝑏 :1if 𝑧𝑗 𝑧𝑖 (z coordinate of the LBB corner of box j is greater than or equal to z coordinate of the LBB corner of0box i)otherwise1if 𝑧𝑗 𝑧𝑖′ (z coordinate of the Left-Bottom-Front (LBB) corner of box j is less than z coordinate of the LeftBottom-Front (LBF) corner of box i)0𝑦𝑖𝑗𝑐 :1otherwiseif π‘₯𝑗′ π‘₯𝑖 (x coordinate of Right-Bottom-Behind (RBB) corner of box j is greater than x coordinate of the LBBcorner of box i)π‘₯𝑖𝑗𝑐 :𝑦𝑖𝑗𝑑 :π‘₯𝑖𝑗𝑑 :π‘§π‘–π‘—π‘Ž :𝑧𝑖𝑗𝑏 :𝑧𝑖𝑗𝑐 :𝑧𝑖𝑗𝑑 :𝐢𝑠1𝑖𝑗 :𝐢𝑠2𝑖𝑗 :𝐢𝑠3𝑖𝑗 :𝐢𝑠4𝑖𝑗 :0otherwise1if π‘₯𝑗′ π‘₯𝑖′ (x coordinate of RBB corner of box j is less than or equal to x coordinate of the RBB corner of box i)0otherwise1if 𝑧𝑗′ 𝑧𝑖 (z coordinate of the LBB corner of box j is greater than z coordinate of the LBF corner of box i)0otherwise1if 𝑧𝑗′ 𝑧𝑖′ (z coordinate of the LBF corner of box j is less than or equal to z coordinate of the LBF corner of box i)0otherwise1if π‘₯𝑖 π‘₯𝑗 π‘₯𝑖′ (both π‘₯π‘–π‘—π‘Ž and π‘¦π‘–π‘—π‘Ž are equal to one)0otherwise1if 𝑧𝑖 𝑧𝑗 𝑧𝑖′ (both π‘₯𝑖𝑗𝑏 and 𝑦𝑖𝑗𝑏 are equal to one)0otherwise1if π‘₯𝑖 π‘₯𝑗′ π‘₯𝑖′ (both π‘₯𝑖𝑗𝑐 and 𝑦𝑖𝑗𝑐 are equal to one)0otherwise1if 𝑧𝑖 𝑧𝑗′ 𝑧𝑖′ (both π‘₯𝑖𝑗𝑑 and 𝑦𝑖𝑗𝑑 are equal to one)0otherwise1if π‘₯𝑖 π‘₯𝑗 π‘₯𝑖′ and 𝑧𝑖 𝑧𝑗 𝑧𝑖′ (both π‘§π‘–π‘—π‘Ž and 𝑧𝑖𝑗𝑏 are equal to one)0otherwise1if π‘₯𝑖 π‘₯𝑗 π‘₯𝑖′ and 𝑧𝑖 𝑧𝑗′ 𝑧𝑖′ (both π‘§π‘–π‘—π‘Ž and 𝑧𝑖𝑗𝑑 are equal to one)0otherwise1if π‘₯𝑖 π‘₯𝑗′ π‘₯𝑖′ and 𝑧𝑖 𝑧𝑗 𝑧𝑖′ (both 𝑧𝑖𝑗𝑐 and 𝑧𝑖𝑗𝑏 are equal to one)0otherwise1if π‘₯𝑖 π‘₯𝑗′ π‘₯𝑖′ and 𝑧𝑖 𝑧𝑗′ 𝑧𝑖′ (both 𝑧𝑖𝑗𝑐 and 𝑧𝑖𝑗𝑑 are equal to one)0otherwise1294

International Journal of Applied Engineering Research ISSN 0973-4562 Volume 12, Number 7 (2017) pp. 1290-1304 Research India Publications. http://www.ripublication.comπ‘₯𝑖′ π‘₯𝑖 𝑀𝑖 𝑋𝑀𝑖 β„Žπ‘– (𝑍𝑀𝑖 π‘Œβ„Žπ‘– 𝑍𝑑𝑖 ) 𝑑𝑖 (1 𝑋𝑀𝑖 𝑍𝑀𝑖 π‘Œ

tree structure. Egeblad and Pisinger (2009) propose a simulated annealing based methodology for the two and three-dimensional knapsack problems, and a three-dimensional knapsack model is presented. The authors present an iterative heuristic approach for the knapsack problem that is based on the sequence triple representation.

Related Documents:

loading stations for loading up to 1000TPH feature the popular and rugged MD30-EV, MV30-EV Series loading spouts and the MA30-EV Series which includes a filter module, air purging system and clean air fan. From a single stand alone loading station with a sepa-rate loading spout and dust collector to an integral Vaculoader, spout positioner

container container container container container networking storage registry security logs & metrics container orchestration & cluster management (kubernetes) fedora / centos / red hat enterprise linux container runtime & packaging (docker) atomic host infrastructure automation & cockpit

container container container container container networking storage registry security logs & metrics container orchestration & cluster management (kubernetes) fedora / centos / red hat enterprise linux container runtime & packaging (docker) atomic host infrastructure automation & cockpit

We offer three types of loading spouts to suit your needs: 1. Tanker loading spout (type NZU) 2. Truck/Wagon/Open pile loading spout (type NZO) 3. Combined loading spout (type NZK) HENNLICH ENGINEERING is the leading manufacturer of telescopic loading spouts - we have been producing them since 1996. We create customized solutions for your .

Dimensions: Main loading door is 8' wide x 10' high (2.4m x 3m). Loading door opens into loading area. This loading area then opens into a hallway, followed by the backstage shop. From main loading door to stage is 60' (18.2m) straight on except for a slight 15 angle between the loading area and backstage shop.

2008 komatsu wa430-6 wheel loader 2013 envirotank 73000 litre skid mounted steel fuel tank 2001 jindo 48 ft high cube container 2003 cimc 40 ft container 2005 quingdao 40 ft container 1996 changzou 20 ft container 2002 evergreen 20 ft container jindo 20 ft storage container 2001 alta-fab wellsite 2010sentag 12 ft x 60 ft 3 unit skid mounted .

Oracle Container Runtime for Docker 19.03 1-2 Oracle Container Runtime for Docker 18.09 1-3 Oracle Container Runtime for Docker 18.03 1-3 Oracle Container Runtime for Docker 17.06 1-4 Docker 17.03 1-5 Docker 1.12 1-6 2 Installing Oracle Container Runtime for Docker Setting Up the Unbreakable Enterprise Kernel 2-1

B.A.G. CORP. SUPER SACK CONTAINER CATALOG S U P E R S A C K C O N T A I N E R D E S I G N S S U P E R S A C K C O N T A I N E R D E S I G N S 7 Spread Strap container Tubular Super Sack container Four-Panel Super Sack container Hardwall container in open position. Barrel Bag container LINER OPTIONS ARE AVAILABLE FOR ALL OF OUR FIBCS .