• Views 587
• Citations 4
• ePub 16
• PDF 298
`Mathematical Problems in EngineeringVolume 2014, Article ID 657170, 16 pageshttp://dx.doi.org/10.1155/2014/657170`
Research Article

## A Quasiphysical and Dynamic Adjustment Approach for Packing the Orthogonal Unequal Rectangles in a Circle with a Mass Balance: Satellite Payload Packing

1School of Information and Engineering, Xiangtan University, Xiangtan, Hunan 411105, China
2School of Mathematics and Computer Science, Xiangtan University, Xiangtan, Hunan 411105, China
3Department of Aeronautics, Xiamen University, Xiamen, Fujian 361005, China

Received 2 June 2014; Revised 28 August 2014; Accepted 22 September 2014; Published 2 December 2014

Copyright © 2014 Ziqiang Li et al. 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

Packing orthogonal unequal rectangles in a circle with a mass balance (BCOURP) is a typical combinational optimization problem with the NP-hard nature. This paper proposes an effective quasiphysical and dynamic adjustment approach (QPDAA). Two embedded degree functions between two orthogonal rectangles and between an orthogonal rectangle and the container are defined, respectively, and the extruded potential energy function and extruded resultant force formula are constructed based on them. By an elimination of the extruded resultant force, the dynamic rectangle adjustment, and an iteration of the translation, the potential energy and static imbalance of the system can be quickly decreased to minima. The continuity and monotony of two embedded degree functions are proved to ensure the compactness of the optimal solution. Numerical experiments show that the proposed QPDAA is superior to existing approaches in performance.

#### 1. Introduction

2D rectangle packing problems are derived from the industry and antiaircraft field [13]. They occur in logistics packing, plate cutting, the layout design of the very large scale integration (VLSI), and satellite modules. They can be divided into unconstrained rectangle packing problems [1] and constrained ones [3]. Both are NP-hard problems and are difficult to be solved. However they have attracted much attention and some packing approaches for different containers have been reported in literatures.

For the 2D rectangle container, the packing approaches mainly include graph theories [47], branch-and-bound methods [810], dynamic planning [11], heuristics [1215], artificial intelligent [16], evolutionary approaches [17], and hybrid approaches [1820]. Regarding the strip container, main packing approaches are branch-and-bound methods [21], heuristics [22, 23], and evolutionary approaches [24, 25]. On the 2D polygon or 3D polyhedron container, the existing packing approaches have heuristics [26], evolutionary approaches [2729], and integer programming [30]. Some scholars are interested in the packing problem of the convex region and have proposed heuristics [26] and branch-and-bound approaches [31].

The layout design problem of the satellite module described in [32] is an important packing problem, which can be transformed into the problem of packing 2D orthogonal unequal rectangles within a circular container with the mass balance (BCOURP). In 1999, Feng et al. [33] built a mathematical model of this problem and analyzed the isomorphism and equivalent intrinsic properties among its layout schemes by using the graph theory and group theory and proposed a theoretical global optimization algorithm. In 2007, Xu et al. [34] defined embedded degree functions between two rectangles and between the rectangle and circular container and presented a compaction algorithm with the particle swarm local search (CA-PSLS). Their idea is that a feasible solution with a smaller envelope radius obtained through the gradient method is taken as an elite individual and the optimal solution is obtained by the PSO iteration. In 2010, Xu et al. [35] suggested a heuristic algorithm ordered by GA (GA-HA, see its algorithm steps in Appendix A and Figure 1). Its key technology is the positioning strategy of constructing the feasible solution. By combining it with GA, the computational efficiency and solution quality are improved.

Figure 1: The sketch map for available positions of the rectangle.

Generally, there exist strong points and deficiencies for each type of approaches.(i)For approaches based on the graph theory, there exists combinatorial explosion when the adjacent topological relation is transformed into the layout diagraph without a size limit for the large-scale layout problem. This is because only two limited relations called the vicinity and distance can be used in pruning branch.(ii)The heuristic method can be used to quickly construct a feasible solution. But it is generally not easy to devise a good heuristic strategy, unless the designer makes a long time painstaking effort and has good luck.(iii)Stochastic algorithms have the global search ability, but there exists the bottleneck of time-costing overlapping area calculation for them [33]. By combining the heuristic method with the stochastic algorithm, their respective advantages can be exerted to the utmost. Based on this mechanism, CA-PSLS and GA-HA are consecutively proposed to solve this problem. According to No Free Lunch Theory [36], how to obtain the knowledge from the problem itself and its area and fuse it into the heuristic and stochastic search mechanism is a way of designing a high performance approach for this problem.

Huang et al. [3740] presented a quasiphysical and quasihuman heuristic algorithm and its variants for the circle packing problem. They obtained excellent results. For BCOURP, Xu et al. [34] proposed CA-PSLS based on embedded degrees between circumcircles of two rectangles and between the container and rectangle’s circumcircle. But due to the discontinuity of the two embedded degree functions, it is difficult to obtain a high quality solution by using CA-PSLS. That is, constructing the continuous rectangular embedded function and exploring a better optimized mechanism are necessary for this problem. Therefore, in this paper, we consider two definitions of monotonous and continuous embedded degrees between two orthogonal rectangles and between an orthogonal rectangle and the container and suggest a dynamic adjustment strategy. We merge them into the proposed QPDAA to improve the solution quality of BCOURP. Numerical experiments will test effectiveness of the considered QPDAA.

The remainder of this paper is organized as follows. The problem statement and mathematical model are in Section 2. The compact and feasible solution strategy and dynamic adjustment strategy are given in Sections 3 and 4, respectively. This algorithm is presented in Section 5. Section 6 is experiments and analysis. The conclusion is shown in Section 7. The final part is acknowledgment.

#### 2. Problem Statement and Mathematical Model

Consider the following two related definitions where and is the number of rectangles.

Definition 1. As shown in Figure 2, let the origin of the Cartesian coordinate system be the center of the container. Let denote the th () rectangle and let (), , , , and be its center, length, width, mass, and direction angle between its long side and the positive direction of the -axis, respectively. Then a layout scheme of rectangles () () can be denoted by ().

Figure 2: The definition of a rectangle.

Definition 2. For a layout scheme , if or 90°, then is an orthogonal rectangle packing scheme, () is an orthogonal rectangle, and the packing is the orthogonal rectangle packing (see Figure 3).

Figure 3: The orthogonal rectangle packing scheme.

Herein this paper considers only orthogonal rectangle packing schemes.

Suppose that the center of each rectangle coincides with its mass center. Let be a rectangle set ; then the mathematical model of this problem can be described as follows. Find a solution and satisfies Formulas (1)–(4). Consider

In Formula (1), denotes the radius of the enveloping circle of the scheme whose circular center is at . Formula (2) indicates that there is no overlap region between two rectangles and . Formula (3) indicates that all rectangles are contained in the container. In Formulas (2) and (3), denotes the interior region of the rectangle . Formula (4) means that the static imbalance of the solution is less than its threshold , where .

#### 3. Compact and Feasible Solution Strategy

Based on the potential energy function of the embedded degree between two circles, Huang et al. [3840] proposed the quasiphysical strategy and its variants for the circle packing problem. Inspired by the quasiphysical idea, we suggest a compact and feasible strategy for BCOURP.

##### 3.1. Embedded Degree Function and Related Properties

Xu et al. [34] defined the embedded degrees between two rectangles and between the rectangle and container by Definitions 3 and 4, respectively.

Definition 3. Let and ( and ) denote the radii of the circumscribed circles of two rectangles and , respectively, and (see Figure 4(a)) the embedded degree between them (see Figure 4(a)); then can be calculated by

Figure 4: Two embedded degree definitions in [34].

Definition 4. Let and () denote the radii of the container and circumscribed circle of the rectangle , respectively, and (see Figure 4(b)) the embedded degree between the rectangle and container; then can be calculated by

Both of the above two embedded degree functions are discontinuous in the critical state from overlapping to separating, as has been discussed by Stoyan and Yaskov [41]. For example, as shown in Figure 5, by moving one rectangle along a direction, two rectangles with the overlap area of (state 1) are changed into and (state 2). By Formula (5), we know that for state 1, but for state 2. So, in Definition 3 is discontinuous where ( and ). Similarly, () in Definition 4 is also discontinuous (see Figure 6). Owing to their discontinuity, it is difficult to select an appropriate step length to obtain the feasible and compact layout scheme for the gradient iteration of CA-PSLS. Inspired by [41], Definitions 5 and 6 are given for the considered QPDAA.

Figure 5: The embedded degree between two rectangles in the critical state.
Figure 6: The embedded degree between the rectangle and container in the critical state.

Definition 5. For two rectangles and ( or , and ) (shown in Figure 7), let and be the radii of their circumscribed circles, respectively, and let denote their embedded degree; then can be calculated by In Formula (7), , . Here, if , , , , and (see Figure 7(a)); otherwise, , , , and (see Figure 7(b)).

Figure 7: Two cases of the embedded degree definition between two rectangles.

In Formula (7), the embedded degree between two rectangles is the moving distance of the rectangle from an overlap state with the stationary to the separation state along the direction from the center to the center . If and are two squares, and the center of is on the diagonal line of and enough close to its center (i.e., and ), the moving distance of from the initial state (shown in Figure 8(a)) to the separate state (shown in Figure 8(b)) along the direction of their diagonal lines is about . Thus in the initial state, their embedded degree is close to the maximal value . In addition, when and/or , and/or . That is, .

Figure 8: The geometric interpretation of definition of the embedded degree between two rectangles.

Definition 6. For the rectangle ( or , ) and the circle container as shown in Figure 7, let denote the embedded degree between the rectangle and container; then can be calculated by
The geometric interpretation of Definition 6 is that when the farthest vertex of the rectangle () from the coordinate origin is within the container, their embedded degree ; otherwise it is the length of the straight line segment pointed by in Figure 9.

Figure 9: The schematic diagram overlapped between the orthogonal rectangle and container.

For embedded degree functions in Definitions 3 and 5, their geometric figures are two curved semi-cone surfaces shown in Figures 10(a) and 10(b), respectively, where it is obvious that there is a gap between the semi-cone surface and plane in Figure 10(a) but there is no gap between them in Figure 10(b). It is not difficult to assert that the difference of geometric figures of two embedded degree functions in Definitions 4 and 6 is the same as the above one. After describing Lemma 7, we propose properties of two embedded degree functions in Definitions 5 and 6, respectively.

Figure 10: The curved semi-cone surface of the embedded degree of two orthogonal rectangles.

Lemma 7. If a binary function is continuous for each variable in a domain, respectively, and is monotonous for the variable or , then the function is continuous in the domain.

Property 1. , let and ( and or ) be two rectangles, and the domain . Then in Definition 5 is a continuous binary function in the domain .

Proof. and , set . From Definition 5, we know that the function is continuous on both and . Here, we prove that it is continuous on the domain .
For , Simultaneously, This is because
According to Lemma 7, the binary function is continuous on the domain . Therefore, is continuous on the domain .

Property 2. and , for the container with the radius and rectangle ( or , ), set and ; then the binary function in Definition 6 is continuous on the domain .

##### 3.2. Extruded Force and Energy Function

In order to quickly decrease the overlapping area of the rectangle packing system, we define the extruded forces between two rectangles and between the rectangle and container.

Definition 8. Let and be two rectangles with the embedded degree ( and ). Then the extruded force between and is calculated by

Definition 9. Let () and denote the rectangle and container. Then the extruded force between and can be calculated by Formula (13), whose direction is the direction from the center of the container to the furthermost rectangular vertex (see Figure 6):
By experiments, we can know that .
So, the extruded resultant force of () in the rectangle packing scheme can be calculated by

Definition 10. Let and ( and ) denote extruded potential energies of with respect to and the container, respectively. Then can be calculated by Formula (15), where and denote two embedded degrees between two rectangles and and between the rectangles and container:

Definition 11. The total extruded potential energy of can be calculated by
Let be the area of the rectangle ; then and () are its absolute and relative extruded potential energies, respectively.

##### 3.3. Compact and Feasible Solution Strategy

By predetermining envelope radius of this problem, the extruded force and direction can be calculated by Formula (13). The extruded force of the envelope circle makes each rectangle close to the center of the container in the direction. For each rectangle, the extruded forces of all others with respect to the rectangle drive away themselves to relieve the pressure in respective direction and their extruded resultant force is calculated by Formula (12). In this paper, the above two steps are used to decrease the overlapping area of the packing scheme and make it compact.

Considering the problem of the low efficiency and possible local optimum (e.g., large static imbalance) of the iteration of two steps in Section 3.3, we propose Property 3 and dynamic adjustment strategy for QPDAA.

##### 4.1. Related Property

For optimizing the static imbalance of the layout scheme, we introduce Property 3.

Property 3. Assume is the mass center of an orthogonal packing scheme of this problem and and is another scheme obtained by interchanging centers of two rectangles with a mass and ( and ) with a mass in . If , , centers and satisfy Formula (17), then , where (): In Formula (17), .

Proof . Consider
Property 3 indicates that, for selecting rectangle , we can find a rectangle in the sector area with an angle (shown as in Figure 11) and interchange them to obtain whose static imbalance is less than . Here, the angle satisfies

Figure 11: The geometric area of Formula (17).

Let ; we consider the following dynamic adjustment strategy.

###### 4.2.1. Rectangle-Interchanging

According to Property 3, the static imbalance of the packing scheme can be decreased by interchanging positions of two rectangles.

Two rectangles, with a smaller mass and with a larger mass , are selected from . If and satisfy Formula (17), then we update by , , .

###### 4.2.2. Rotation and Off-Trap

(i) A rectangle with a larger pain degree is found out from and is rotated 90° round its center counterclockwise direction to relieve its pain. (ii) A rectangle with a larger pain degree is found out from and is moved to such a place in the container where its pain degree is smaller. We know that the role of (ii) is similar to a construction of the nonisomorphic layout pattern [42, 43].

###### 4.2.3. Center Translation

If its mass center , then .

#### 5. The Proposed Algorithm

Through an organic combination of the compact and feasible solution strategy and dynamic adjustment strategy, we present QPDAA for BCOURP.

Let and be the predetermined value and an allowable maximum of the envelope radius of the solution, respectively; is the number of rectangles. , , and denote the mass center, envelope radius, and extruded potential energy of the packing scheme , respectively. () is the extruded potential energy of the rectangle ; is the step length; is the maximum translation times. Then steps of the proposed QPDAA are shown in Algorithm 1.

Algorithm 1

#### 6. Experiments and Analysis

##### 6.1. Experiments

The proposed QPDAA is coded in VC++ 6.0 and carried on a Pentium 3 GHZ PC with 512 MB memory. CA-PSLS [34], GA-HA [35] are coded in VC++ 6.0 and two algorithms are carried on a Pentium 1.83 GHZ with 512 MB memory; IGA [44] is carried on an IBM 586 166 MHz.

Experiment 1. Five examples are taken from [33, 35] and are used in testing the performance of the proposed QPDAA. Data of all examples are shown in Table 1. For the proposed algorithm, we take , , , , , and , respectively. For Examples , we take = 11.4, 14.3, 17.5, 22.3, and 115.5 and take , , , , and , respectively. Running the proposed QPDAA 30 times for each example (the success rate is 100%), we show its average running time, average envelope radius, standard variance of the radius, and the maximum and minimum envelope radii in Table 2; for each example, its layout scheme diagraph is shown in Figure 12; the other data in Table 2 is taken from [3335]. The optimal layout schemes of the proposed QPDAA are shown in Tables 3, 4, 5, and 6 for 5 examples.

Table 1: Parameters of rectangles for five layout examples.
Table 2: Performance comparisons of four algorithms.
Table 3: The layout schemes of the proposed QPDAA for Examples  1 and  2.
Table 4: The layout schemes of the proposed QPDAA for Example  3.
Table 5: The layout scheme of the proposed QPDAA for Example  4.
Table 6: The layout scheme of the proposed QPDAA for Example  5.
Figure 12: The layout diagraphs of five examples for the proposed QPDAA.

Experiment 2. For testing the effects of and on the minimal radius and running time with the proposed QPDAA, we take another set of and (five examples and other parameters are the same as those of  Experiment 1) and run the proposed QPDAA procedure 30 times for each example. The minimal radii and running times are given in Table 7. Their layout schemes and layout diagraphs are shown in Tables 8, 9, 10, and 11 and Figure 13. It can be found from Table 7 that changing values of and we can obtain the packing scheme with a smaller envelop radius, but it costs more time for each example. So, values and in Experiment 1 can be applied to make a tradeoff between the computational effectiveness and solution quality.

Table 7: The effect of parameters and on the optimal radius and running time for the proposed QPDAA.
Table 8: The optimal layout schemes of Examples  1 and  2 for the proposed QPDAA.
Table 9: The optimal layout schemes of Example  3 for the proposed QPDAA.
Table 10: The optimal layout scheme of Example  4 for the proposed QPDAA.
Table 11: The optimal layout scheme of Example  5 for the proposed QPDAA.
Figure 13: The optimal layout scheme diagraphs of five examples for the proposed QPDAA.

In order to further test the effectiveness of the proposed QPDAA, we consider Experiment 3.

Experiment 3. Numbers of rectangles (generated randomly) of three examples are 50, 55, and 60 and their lengths and widths are between 20 and 40. For the proposed QPDAA, we take , , , , , , = 125.9, 133.8, 137.2, and = 126.8, 134.6, 138.8 for Examples , respectively. For GA + HA, the population size, mutation probability, and max number of the iteration are 30, 0.125, and 50, respectively. By running HA + GA and the proposed QPDAA 30 times for three examples, respectively, their optimal envelop radii and average times are shown in Table 12, respectively. The optimal packing scheme diagrams of GA + HA and the proposed QPDAA are shown in Figures 14(a)14(c) and Figures 14(d)14(f). We can know from Table 12 that both the solution quality and computational quality of the proposed QPDAA are obviously higher than those of GA + HA.

Table 12: The effect of and on the optimal radius and running time for the proposed QPDAA.
Figure 14: Packing scheme diagrams of GA + HA and the proposed QPDAA for Experiment 3.

Note that, in Experiment 3, the procedure of GA + HA is coded by author and is carried on a Pentium 3 GHZ with 512 MB memory.

##### 6.2. Analysis

From data of Tables 2, 7, and 12, we know that the solution quality of the proposed QPDAA algorithm is higher than those of CA-PSLS and GA-HA. Compared with those of CA-PSLS, embedded degree functions of the proposed QPDAA can make the layout scheme more compact. Due to the fixed candidate positions of GA + HA, it is difficult to find the best position for some rectangles close to the marginal region of the container. The deficiency of a mechanism of decreasing the static balance in the process of packing rectangles also limits the solution quality of GA + HA. The experimental results illustrate the effectiveness of the proposed QPDAA.

The computational efficiency of the proposed QPDAA is one magnitude higher than those of CA-PSLS. There are two reasons. (i) Owing to orthogonal packing, searching the optimal solution in the 2D solution of the proposed QPDAA is easier than that in the 3D solution space of CA-PSLS. (ii) In order to improve the solution quality of CA-PSLS, PSO is used to optimize the feasible solution obtained through the gradient method based on two discontinuous embedded degree functions, but the proposed QPDAA need not do the PSO optimization without reducing its solution quality. Except for Example 1, the computational efficiency of the proposed QPDAA is higher than that of GA + HA. And with the increase of the size of the packing problem, the advantage of the proposed QPDAA is more obvious. This is because computational complexities of the extruded resultant force and potential energy are in this paper, but the computational complexity of noninterference judgment is . In addition, for GA + HA, with the increasing of the number of rectangles, the number of candidate positions of each rectangle increases dramatically. These reasons lead to the computational efficiency of the proposed QPDAA higher than that of GA + HA for the BCOURP with a large size.

#### 7. Conclusions

Taking the layout design of a satellite module as the application background, we have proposed the QPDAA for the BCOURP problem in this paper. Two continuous embedded functions between orthogonal rectangles and between the rectangle and container are constructed to overcome the weakness of embedded functions in [34]. And the suggestion of the extruded resultant force formula and the potential energy function of the rectangle packing system based on the proposed embedded functions make solving the BCOURP problem simple and effective as solving the circle packing problem [3740]. The proposed dynamic adjustment strategy can quickly decrease the static imbalance of the packing scheme and make the iteration skip the local optimum. The experiment results show that the proposed QPDAA is superior to existing algorithms in performance for the BCOURP problem, especially for the BCOURP problem with the large size. The next work is to extend the above algorithm into solving the 3D satellite module payload packing problem.

#### A. HA + GA [35]

Input the length, width, and mass of the rectangle () in turn, initialize the number of the maximal iteration times, and generate their placing sequence set .

Step 1. Set .

Step 2. Set the center of at and its long side is parallel to -axis, , .

Step 3. For , calculate centers and direction angles of 16 candidate positions (see Figure 1) of the rectangle with respect to the rectangle . From its candidate positions eliminate unfeasible ones and calculate the optimal one (i.e., compared with other feasible candidate positions, it makes the packing scheme of the first rectangles have less envelop radius).

Step 4. If then update the current packing scheme, , and go to Step 3; otherwise, go to Step 5.

Step 5. If then update the optimal packing scheme and use GA to generate a new placing sequence set and go to Step 2; otherwise, go to Step 6.

Step 6. Output the optimal packing scheme and envelop radius; algorithm ends.

#### B. Results of Experiment 2

See Tables 8, 9, 10, and 11 and Figure 13.

#### Conflict of Interests

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

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 61272294) and Research Foundation of Education Bureau of Hunan Province, China (Grant no. 11A120), and the Construct Program of the Key Discipline in Hunan Province. And the authors are also grateful to the anonymous referees for their advice and reviews for this paper.

#### References

1. K. He, W. Huang, and Y. Jin, “An efficient deterministic heuristic for two-dimensional rectangular packing,” Computers & Operations Research, vol. 39, no. 7, pp. 1355–1363, 2012.
2. N. Lesh, J. Marks, A. McMahon, and M. Mitzenmacher, “Exhaustive approaches to 2D rectangular perfect packings,” Information Processing Letters, vol. 90, no. 1, pp. 7–14, 2004.
3. E.-M. Feng and X.-L. Wang, “An optimization model for the layout of a group of cuboids in a satellite module and its global optimization algorithm,” Operation Research Transaction, vol. 5, no. 3, pp. 71–77, 2001 (Chinese).
4. J. Roth and R. Hasimshony, “Comparison of existing three-room apartment plans with computer-generated layouts,” Planning and Design, vol. 14, no. 2, pp. 149–161, 1987.
5. J. Leung, “A new graph-theoretic heuristic for facility layout,” Management Science, vol. 38, no. 4, pp. 594–606, 1992.
6. J. M. V. de Carvalho, “Exact solution of bin-packing problems using column generation and branch-and-bound,” Annals of Operations Research, vol. 86, pp. 629–659, 1999.
7. R. Macedo, C. Alves, and J. M. V. de Carvalho, “Arc-flow model for the two-dimensional guillotine cutting stock problem,” Computers and Operations Research, vol. 37, no. 6, pp. 991–1001, 2010.
8. Y.-D. Cui, C. L. Zhang, and Y. Zhao, “A continued fractions and branch-and-bound algorithm for generating cutting patterns with equal rectangles,” Journal of Computer-Aided Design & Computer Graphics, vol. 16, no. 2, pp. 252–256, 2004.
9. F. Clautiaux, A. Jouglet, J. Carlier, and A. Moukrim, “A new constraint programming approach for the orthogonal packing problem,” Computers & Operations Research, vol. 35, no. 3, pp. 944–959, 2008.
10. K. Yoon, S. Ahn, and M. Kang, “An improved best-first branch-and-bound algorithm for constrained two-dimensional guillotine cutting problems,” International Journal of Production Research, vol. 51, no. 6, pp. 1680–1693, 2013.
11. E. G. Birgin, R. D. Lobato, and R. Morabito, “Generating unconstrained two-dimensional non-guillotine cutting patterns by a recursive partitioning algorithm,” Journal of the Operational Research Society, vol. 63, no. 2, pp. 183–200, 2012.
12. Y.-L. Wu, W. Huang, S.-C. Lau, C. Wong, and G. H. Young, “An effective quasi-human based heuristic for solving the rectangle packing problem,” European Journal of Operational Research, vol. 141, no. 2, pp. 341–358, 2002.
13. D.-F. Zhang, S.-H. Han, and W.-G. Ye, “Bricklaying heuristic algorithm for the orthogonal rectangular packing problem,” Chinese Journal of Computers, vol. 31, no. 3, pp. 509–514, 2008.
14. R. E. Korf, M. D. Moffitt, and M. E. Pollack, “Optimal rectangle packing,” Annals of Operations Research, vol. 179, no. 1, pp. 261–295, 2010.
15. C. Charalambous and K. Fleszar, “A constructive bin-oriented heuristic for the two-dimensional bin packing problem with guillotine cuts,” Computers and Operations Research, vol. 38, no. 10, pp. 1443–1451, 2011.
16. S. Polyakovsky and R. M'Hallah, “An agent-based approach to the two-dimensional guillotine bin packing problem,” European Journal of Operational Research, vol. 192, no. 3, pp. 767–781, 2009.
17. S. Khebbache1, C. Prins, and A. Yalaoui, “Iterated local search algorithm for the constrained two-dimensional non-guillotine cutting problem,” Journal of Industrial and Systems Engineering, vol. 2, no. 3, pp. 164–179, 2008.
18. J. F. Gonçalves and M. G. C. Resende, “A parallel multi-population genetic algorithm for a constrained two-dimensional orthogonal packing problem,” Journal of Combinatorial Optimization, vol. 22, no. 2, pp. 180–201, 2011.
19. E. G. Birgin, R. D. Lobato, and R. Morabito, “An effective recursive partitioning approach for the packing of identical rectangles in a rectangle,” Journal of the Operational Research Society, vol. 61, no. 2, pp. 306–320, 2010.
20. M. Dolatabadi, A. Lodi, and M. Monaci, “Exact algorithms for the two-dimensional guillotine knapsack,” Computers and Operations Research, vol. 39, no. 1, pp. 48–53, 2012.
21. S. Martello, M. Monaci, and D. Vigo, “An exact approach to the strip-packing problem,” INFORMS Journal on Computing, vol. 15, no. 3, pp. 310–319, 2003.
22. B. S. Baker, J. Coffman, and R. . Rivest, “Orthogonal packings in two dimensions,” SIAM Journal on Computing, vol. 9, no. 4, pp. 846–855, 1980.
23. D. Liu and H. Teng, “An improved BL-algorithm for genetic algorithm of the orthogonal packing of rectangles,” European Journal of Operational Research, vol. 112, no. 2, pp. 413–420, 1999.
24. L. H. W. Yeung and W. K. S. Tang, “Strip-packing using hybrid genetic approach,” Engineering Applications of Artificial Intelligence, vol. 17, no. 2, pp. 169–177, 2004.
25. D. Ye, X. Han, and G. Zhang, “A note on online strip packing,” Journal of Combinatorial Optimization, vol. 17, no. 4, pp. 417–423, 2009.
26. A. Cassioli and M. Locatelli, “A heuristic approach for packing identical rectangles in convex regions,” Computers and Operations Research, vol. 38, no. 9, pp. 1342–1350, 2011.
27. J. Błazewicz, P. Hawryluk, and R. Walkowiak, “Using a tabu search approach for solving the two-dimensional irregular cutting problem,” Annals of Operations Research, vol. 41, no. 4, pp. 313–325, 1993.
28. V. Petridis, S. Kazarlis, and A. Bakirtzis, “Varying fitness functions in genetic algorithm constrained optimization: the cutting stock and unit commitment problems,” IEEE Transactions on Systems, Man, and Cybernetics. Part B: Cybernetics, vol. 28, no. 5, pp. 629–640, 1998.
29. Y. Chen, M. Tang, R. Tong, and J. Dong, “Packing of polygons using genetic simulated annealing algorithm,” Journal of Computer-Aided Design and Computer Graphics, vol. 15, no. 5, pp. 598–609, 2003 (Chinese).
30. R. Andrade and E. G. Birgin, “Symmetry-breaking constraints for packing identical rectangles within polyhedra,” Optimization Letters, vol. 7, no. 2, pp. 375–405, 2013.
31. E. G. Birgin and R. D. Lobato, “Orthogonal packing of identical rectangles within isotropic convex regions,” Computers and Industrial Engineering, vol. 59, no. 4, pp. 595–602, 2010.
32. H.-F. Teng, S. L. Sun, W. H. Ge, and W. X. Zhong, “Layout optimization for the dishes installed on a rotating table—the packing problem with equilibrium behavioural constraints,” Science in China A: Mathematics, Physics, Astronomy, vol. 37, no. 10, pp. 1272–1280, 1994.
33. E. Feng, X. Wang, X. Wang, and H. Teng, “An algorithm of global optimization for solving layout problems,” European Journal of Operational Research, vol. 114, no. 2, pp. 430–436, 1999.
34. Y.-C. Xu, R.-B. Xiao, and M. Amos, “Particle swarm algorithm for weighted rectangle placement,” in Proceedings of the 3rd International Conference on Natural Computation (ICNC '07), China IEEE Press, Haikou, China, August 2007.
35. Y.-C. Xu, F.-M. Dong, Y. Liu, and R.-B. Xiao, “Genetic algorithm for rectangle layout optimization with equilibrium constraints,” Pattern Recognition and Artificial Intelligence, vol. 23, no. 6, pp. 794–801, 2010 (Chinese).
36. D. H. Wolpert and W. G. Macready, “No free lunch theorems for optimization,” IEEE Transactions on Evolutionary Computation, vol. 1, no. 1, pp. 67–82, 1997.
37. W.-Q. Huang and R.-C. Xu, “Two Quasi-human strategies for the circle packing problem,” Science China (Series E), vol. 29, no. 4, pp. 347–353, 1999.
38. W.-Q. Huang and M. Chen, “Note on: an improved algorithm for the packing of unequal circles within a larger containing circle,” Computers and Industrial Engineering, vol. 50, no. 3, pp. 338–344, 2006.
39. W.-Q. Huang and T. Ye, “Quasi-physical algorithm for the equal circle packing problem,” Journal System Science & Math, vol. 50, no. 3, pp. 993–1001, 2008.
40. W.-Q. Huang and T. Ye, “Quasi-physical global optimization method for solving the equal circle packing problem,” Science China: Information Sciences, vol. 54, no. 7, pp. 1333–1339, 2011.
41. Y. G. Stoyan and G. N. Yaskov, “Mathematical model and solution method of optimization problem of placement of rectangles and circles taking into account special constraints,” International Transactions in Operational Research, vol. 5, no. 1, pp. 45–57, 1998.
42. H.-F. Teng, Z.-Q. Li, Y.-J. Shi, and Y.-S. Wang, “An approach to constructing isomorphic or non-isomorphic layout pattern,” Chinese Journal of Computers, vol. 29, no. 6, pp. 985–991, 2006.
43. H.-F. Teng, Y. Chen, W. Zeng, Y.-J. Shi, and Q.-H. Hu, “A dual-system variable-grain cooperative coevolutionary algorithm: satellite-module layout design,” IEEE Transactions on Evolutionary Computation, vol. 14, no. 3, pp. 438–455, 2010.
44. E.-M. Feng, Z.-H. Gong, C.-Y. Liu, and Z. Xu, “Improved GA for satellite module layout problem with performance constraints,” Journal of Dalian University of Technology, vol. 45, no. 3, pp. 457–463, 2005.