Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 317697, 6 pages

http://dx.doi.org/10.1155/2014/317697

Research Article

Integer Programming Formulations for Approximate Packing Circles in a Rectangular Container

Igor Litvinchev1,2 and Edith Lucero Ozuna Espinosa²

¹Complex Systems Department, Computing Center, Russian Academy of Sciences, 119333 Moscow, Russia

²Faculty of Mechanical and Electrical Engineering, Nuevo Leon State University, 66450 Monterrey, NL, Mexico

Received 23 December 2013; Accepted 15 February 2014; Published 25 March 2014

Academic Editor: Vyacheslav Kalashnikov

Copyright © 2014 Igor Litvinchev and Edith Lucero Ozuna Espinosa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A problem of packing a limited number of unequal circles in a fixed size rectangular container is considered. The aim is to maximize the (weighted) number of circles placed into the container or minimize the waste. This problem has numerous applications in logistics, including production and packing for the textile, apparel, naval, automobile, aerospace, and food industries. Frequently the problem is formulated as a nonconvex continuous optimization problem which is solved by heuristic techniques combined with local search procedures. New formulations are proposed for approximate solution of packing problem. The container is approximated by a regular grid and the nodes of the grid are considered as potential positions for assigning centers of the circles. The packing problem is then stated as a large scale linear 0-1 optimization problem. The binary variables represent the assignment of centers to the nodes of the grid. Nesting circles inside one another is also considered. The resulting binary problem is then solved by commercial software. Numerical results are presented to demonstrate the efficiency of the proposed approach and compared with known results.

1. Introduction

Packing problems generally consist of packing of a certain number of items of known dimensions into one or more large objects or containers so as to minimize a certain objective (e.g., the unused part of the objects or waste). The shape of items and containers may vary from a circle, a square, a rectangle, and so forth.

This problem has been applied in different areas, such as the coverage of a geographical area with cell transmitters, storage of a cylindrical drums into containers or stocking them into an open area, packaging bottles or cans into the smallest box, planting trees in a given region so as to maximize the forest density and the distance between the trees, and so forth [1–3]. One can find other applications in the motor cycle industry, circular cutting, communication networks, facility location, and dashboard layout [4–8].

In this paper we address the problem of packing a set of circular items in a rectangular container. There are two principal types of objectives that have been used in the literature: (a) regard the circles (not necessary equal) as being of fixed size and the container as being of variable size and (b) regard the circles and the container as being of fixed size and minimize “waste.”

Examples of the first approach include the following [9].

(i)For the square container minimize the length of the side and hence minimize the perimeter and area of the square.(ii)Minimize the perimeter of the rectangle.(iii)Minimize the area of the rectangle.(iv)Considering one dimension of the rectangle as fixed, minimize the other dimension. Problems of this type are often referred to as strip packing problems (or as circular open dimension problems).

For the second approach various definitions of the waste can be used. The waste can be defined in relation to circles not packed (e.g., the number of unpacked circles or the perimeter/area of unpacked circles) or introducing a value associated with each circle that is packed (e.g., area of the circles packed), and so forth.

Many variants of packing circular objects in the plane have been formulated as nonconvex (continuous) optimization problems with decision variables being coordinates of the centres. The nonconvexity is mainly provided by no overlapping conditions between circles. These conditions typically state that the Euclidean distance separating the centres of the circles is greater than a sum of their radii. The nonconvex problems can be tackled by available nonlinear programming (NLP) solvers; however, most NLP solvers fail to identify global optima. Thus, the nonconvex formulation of circular packing problem requires algorithms which mix local searches with heuristic procedures in order to widely explore the search space. It is impossible to give a detailed overview on the existing solution strategies and numerical results within the framework of a single short paper. We will refer the reader to review papers presenting the scope of techniques and applications for the circle packing problem (see, e.g., [1, 2, 10–13] and the references therein).

In this paper we propose a new formulation for approximate solution of packing problems based on using a regular grid to approximate the container. The nodes of the grid are considered as potential positions for assigning centers of the circles. The packing problem is then stated as a large scale linear 0-1 optimization problem. The binary variables represent the assignment of centers to the nodes of the grid. The resulting binary problem is then solved by commercial software. To the best of our knowledge, the idea to use a grid was first implemented by Beasley [14] in the context of cutting problems. Recently this approach was applied in [15, 16] for packing problems. This work is a continuation of [16].

2. The Model

Suppose we have nonidentical circles of known radius , . Let at most circles be available for packing, and at least of them have to be packed. Denote by the node points of a regular grid covering the rectangular container. Let be the grid points lying on the boundary of the container. Denote by the Euclidean distance between points and of the grid. Define binary variables if centre of a circle is assigned to the point ; , otherwise.

In what follows we will distinguish two cases of circle packing, depending on whether nesting circles inside one another is permitted or not. To the best of our knowledge, nesting problem was first mentioned in [6] in the context of packing pipes of different diameters into a shipping container. Compared to the standard packing, packing with nesting is much less investigated.

Consider first the problem without nesting. In order to let the circle assigned to the point be nonoverlapping with other circles being packed, it is necessary that for , , such that . For fixed , let . Let be the cardinality of . Then the problem of maximizing the area covered by the circles can be stated as follows:

Constraints (2) ensure that the number of circles packed is between and ; constraints (3) ensure that at most one centre is assigned to any grid point; constraints (4) ensure that the point cannot be a centre of the circle if the distance from to the boundary is less than ; pairwise constraints (5) guarantee that there is no overlapping between the circles; constraints (6) represent the binary nature of variables.

Similar to plant location problems [17] we can state nonoverlapping conditions in a more compact form. Summing up pairwise constraints (5) over we get

Note that constraints similar to (7) were used in [15] for packing equal circles.

Proposition 1. Constraints (5) and (6) are equivalent to constraints (6) and (7).

Proof. If (5) are fulfilled, then obviously (7) hold by construction. If (7) are fulfilled, then for we have and hence for as in (5). If in (7), then holds for all .

Thus the problem (1)–(6) is equivalent to the problem (1)–(4), (6), and (7). To compare two equivalent formulations, let

Proposition 2. .

Proof. Since constraints of are obtained by summing up some constraints of , then . To show that we need to find a point in that is not in .

This point can be constructed as follows. By the definition, and hence if , then . Choose and such that , . Set to zero all the variables except , . Obviously all constraints (5) corresponding to zero variables are fulfilled. Define , to fulfil the two remaining constraints as equalities:

with , . The corresponding solution is

with

Thus this point violates corresponding constraint (5) in and hence as desired.

As follows from Proposition 2, the pairwise formulation (1)–(6) is stronger than the compact one [17].

By the definition, and hence if , then . Thus a half of the constraints in (5) are redundant:

The redundant constraints can be eliminated giving a reduced pairwise nonoverlapping formulation. The overall set of the reduced constraints is independent of whether we will eliminate the first constraint above or the second. However, the way of eliminating redundant constraints will affect the corresponding compact formulation obtained by summing up the reduced constraints.

To consider nesting circles inside one another, we only need to modify the nonoverlapping constraints. In order to let the circle assigned to the point be nonoverlapping with other circles being packed (including circles places inside this circle), it is necessary that for , , such that . Note that the later condition is always fulfilled for (), such that only smaller circles can be placed inside a given circle. For fixed , let . Then the problem of packing circles with nesting can be stated as follows:

Here weighting coefficients may be associated with the area of circles and/or represent the relative importance of subsets of the container.

3. Computational Results

A rectangular uniform grid was used in the numerical experiments, such that all grid points are defined by the grid points on its edges. Let be a horizontal dimension (length) and let be a vertical dimension (width) of the container; let be a number of the equidistant grid points on the horizontal edge of the container, while is a number of the equidistant grid points on its vertical edge. Hence the grid has node points ().

All optimization problems were solved by the system CPLEX 12.5. The runs were executed on a DELL Power Edge T410, Intel Xeon 2.53 Ghz and 16 Gb RAM.

First, we compare different formulations for the case of packing equal circles. The same set of 9 instances as in [15, Table 3] was used for the experiment. The number of circles available for packing was not limited. The following formulations were considered: pairwise formulation (1)–(6) (complete), reduced formulation (1)–(6) without redundant constraints (half), compact formulation (1)–(4), (6), and (7) (compact), and compact formulation obtained from (1)–(6) without redundant constraints (compact half). The results of the numerical experiment are presented in Table 1. Here the first four columns present circle radius, dimension of the quadratic cell used to construct the grid (), the number of grid points (), and the number of circles packed. For all formulations the same number of circles packed was achieved. The last four columns give CPU time for different formulations (in seconds). For all problem instances mipgap = 0 was set for running CPLEX. Asterisks in the column indicate that computation was interrupted after the computational time exceeded 15-hour CPU time.

Table 1: Results of the numerical experiments.

Circle radius

Problem dimension
Circle number
Complete
Half
Compact
Compact half

0.625
0.15625
1403
10
6.4
125.5
144.8
110.6

0.5625
0.0703125
2449
13
50.4
647.0
5034.5
1890.6

0.5
0.125
697
18
2.6
3.2
5.2
24.1

0.4375
0.0546875
3666
21
849.2
310.8
7459.9
3690.5

0.375
0.046875
1425
32
50.0
21.4
3873.5
847.7

0.3125
0.078125
2139
45
403.8
183.9
6514.3
4451.8

0.275
0.06875
2880
61
1032.4
415.3
*
6985.4

0.25
0.0625
3649
74
1234.8
535.5
14645.9
7840.6

0.1875
0.046875
6897
140
1427.3
725.9
*
8765.2

As we can see from Table 1, CPU time for the complete formulation is typically much lower than for the compact formulation, especially for large problems. Note that according to Proposition 2, complete (pairwise) formulation is stronger than the compact one. Eliminating redundant constraints typically (but not always) reduces CPU time, although for complete formulation eliminating redundancy does not change the continuous relaxation. This effect is well known for users of CPLEX since simply rearranging constraints may result in increase/decrease of computational time, depending on the path selected by branch and bound algorithm. However, we may conclude that eliminating redundant constraints is useful for large problems.

Figures 1, 2, 3, and 4 present packing pictures, for instance, from Table 1 with radii 0.375, 0.3125, 0.4375, and 0.275, respectively.

Figure 1: .

Figure 2: .

Figure 3: .

Figure 4: .

Figures 5 and 6 present the packing pictures obtained for the container with and , , and using pairwise formulation (1)–(6). We used for Figure 5 and for Figure 6. The solution presented in Figure 5 was obtained in 53.239 sec with mipgap = 0, while for Figure 6 the computation was interrupted after the computational time exceeded 12-hour CPU time with mipgap = 8% achieved.

Figure 5

Figure 6

Figures 7 and 8 present packing pictures obtained by the formulation (13) for the quadratic container with and four radii , , , and with nesting permitted. Figure 7 corresponds to the grid with , while was used for Figure 8. We see how packing approximation changes with the parameters of the grid. In Figure 7 the following circles were packed: 741 circles , 52 , 20 , and 3 . In Figure 8 we have 932 circles , 77 , 13 , and 5  and there is still room for more circles. For the case of Figures 7 and 8 the computations were interrupted after the computational time exceeded 12-hour CPU time with mipgap = 67% and mipgap = 83%, correspondingly.

Figure 7

Figure 8

4. Conclusions

Different integer formulations were proposed for approximate solution of the circle packing problem. We demonstrate that the pairwise formulation is stronger than the compact one obtained by summing up the nonoverlapping constraints. The presented approach can be easily generalized to the three- (and more) dimensional case and to different shapes of the container, including irregulars. We also proposed a formulation permitting nesting circles inside one another. This problem was mentioned in [7] in the context of packing pipes of different diameters into a shipping container and has not received much attention so far. An interesting direction for the future research is to study the use of Lagrangian relaxation and corresponding heuristics [3] to cope with large dimension of arising problems. Another area for future study is the use of valid inequalities for strengthening the formulations, as well as combining continuous and integer formulations in the solution process. Some complements in these directions are in course [16].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was partially supported by Grants from RFBR, Russia (12_01_00893_a), and CONACYT, Mexico (167019). The authors would like to thank anonymous referees for their constructive comments and suggestions.

References

E. G. Birgin and J. M. Gentil, “New and improved results for packing identical unitary radius circles within triangles, rectangles and strips,” Computers and Operations Research, vol. 37, no. 7, pp. 1318–1327, 2010.

M. Hifi and R. M'Hallah, “A literature review on circle and sphere packing problems: models and methodologies,” Advances in Operations Research, vol. 2009, Article ID 150624, 22 pages, 2009.

I. Litvinchev, S. Rangel, and J. Saucedo, “A lagrangian bound for many-to-many assignment problems,” Journal of Combinatorial Optimization, vol. 19, no. 3, pp. 241–257, 2010.

E. Baltacioǧlu, J. T. Moore, and R. R. Hill Jr., “The distributor's three-dimensional pallet-packing problem: a human intelligence-based heuristic approach,” International Journal of Operational Research, vol. 1, no. 3, pp. 249–266, 2006.

E. G. Birgin, J. M. Martínez, and D. P. Ronconi, “Optimizing the packing of cylinders into a rectangular container: a nonlinear approach,” European Journal of Operational Research, vol. 160, no. 1, pp. 19–33, 2005.

I. Castillo, F. J. Kampas, and J. D. Pintér, “Solving circle packing problems by global optimization: numerical results and industrial applications,” European Journal of Operational Research, vol. 191, no. 3, pp. 786–802, 2008.

H. J. Fraser and J. A. George, “Integrated container loading software for pulp and paper industry,” European Journal of Operational Research, vol. 77, no. 3, pp. 466–474, 1994.

J. A. George, “Multiple container packing: a case study of pipe packing,” Journal of the Operational Research Society, vol. 47, no. 9, pp. 1098–1109, 1996.

C. O. López and J. E. Beasley, “A heuristic for the circle packing problem with a variety of containers,” European Journal of Operational Research, vol. 214, pp. 512–525, 2011.

H. Akeb and M. Hifi, “Solving the circular open dimension problem using separate beams and look-ahead strategies,” Computers & Operations Research, vol. 40, no. 5, pp. 1243–1255, 2013.

M. H. Correia, J. F. Oliveira, and J. S. Ferreira, “Cylinder packing by simulated annealing,” Pesquisa Operacional, vol. 20, pp. 269–286, 2000.

C. O. Lopez and J. E. Beasley, “Packing unequal circles using formulation space search,” Computers & Operations Research, vol. 40, pp. 1276–1288, 2013.

Y. G. Stoyan and G. N. Yaskov, “Packing congruent spheres into a multi-connected polyhedral domain,” International Transactions in Operational Research, vol. 20, no. 1, pp. 79–99, 2013.

J. E. Beasley, “An exact two-dimensional non-guillotine cutting tree search procedure,” Operations Research, vol. 33, no. 1, pp. 49–64, 1985.

S. I. Galiev and M. S. Lisafina, “Linear models for the approximate solution of the problem of packing equal circles into a given domain,” European Journal of Operational Research, vol. 230, no. 3, pp. 505–514, 2013.

I. Litvinchev and E. L. Ozuna, “Packing circles in a rectangular container,” in Proceedings of the International Congress on Logistics and Supply Chain, pp. 24–25, Queretaro, Mexico, October 2013.

L. A. Wolsey, Integer Programming, Wiley, New York, NY, USA, 1999.