Available online at www.sciencedirect.com
ScienceDirect Procedia CIRP 33 (2015) 87 – 92
9th CIRP Conference on Intelligent Computation in Manufacturing Engineering - CIRP ICME '14
A Novel Bi-level Continuous Formulation for the Cellular Manufacturing System Facility Layout Problem Maral Zafar Allahyaria, Ahmed Azabb*
a,b
Product Life Cycle Research Lab, Dept. of Industrial & Manufacturing Systems Engineering, University of Windsor, 401 Sunset Ave, Windsor, N9B 3P4, Canada
* Corresponding author. Tel.: +1-519-253-3000 EXT: 4958; fax: +1 (519) 973-7062. E-mail address:
[email protected]
Abstract Cellular Manufacturing Systems (CMSs) play an important role in today’s small- to mid-size production enterprises. The facility layout problem for such manufacturing system is a group-one that includes both intercellular and intracellular relationships. A novel continuous formulation has been developed for the problem to model manufacturing shops with vertical and horizontal aisles and to eliminate any possible overlap between machine tools as well that between cells. The overall approach adopted is a bi-level one; initially an upper-level leader facility layout problem is being solved for each cell at a time; then, a lower-level follower facility layout problem (FLP) is being solved after at the shop level determining the overall layout of the shop, where the position of the different cells are to be determined. A case study from the local machining industry has been utilized to verify the model. The Authors. Authors. Published Publishedby byElsevier ElsevierB.V. B.V. This is an open access article under the CC BY-NC-ND license © 2014 The Peer-review under responsibility of the International Scientific Committee of “9th CIRP ICME Conference". (http://creativecommons.org/licenses/by-nc-nd/4.0/). Selection and peer-review under responsibility of the International Scientific Committee of “9th CIRP ICME Conference” Keywords: Cellular Manufacturing Systems; Unequal Facility Layout Problem; Mathematical Modelling; Overlap Elimination; Aisle Modelling
1. Introduction The facility layout problem can be defined as an optimization problem that tries to make layouts more efficient by taking into account various interactions between facilities and material handling systems while designing layouts [1]. There is a lot of interest in optimizing FLP (Facility Layout Problem); however, the majority of which have attempted to solve FLP by using a discrete approach [2,3,4]. This approach attempts to assigning n facilities with the same shape and size to n distinct predetermined locations. The main critique this widely used modeling scheme has received is the lack of geometric constraints to account for the possible unequal size of the different facilities and the lack of knowledge of the possible locations of the different facilities (cells and machines). Discrete representations are not suited to represent the exact position of facilities in the plant site, the orientation of facilities, pick-up and drop-off points or clearance between facilities [5]. A continuous mathematical representation has been introduced, where the problem becomes the optimal
allocation of the two-dimensional coordinates of facilities within a shop. This has been often addressed as Mixed Integer Programming Problems [6]. The continuous approach tries to find optimal arrangement of facilities within the planar site, and yet not overlap each other. Cellular manufacturing system (CMS) layout has recently begun to receive heightened attention worldwide. An effective CMS implementation help any company improve machine utilization and quality; it also makes reduction in setup time, work-in-process inventory, material handling cost, part makespan, and expediting costs [7,8,9]. In this paper, the layout problem of CMS is being tackled. A Bilevel Mathematical Programming approach has been adopted, in which two-tier mathematical formulations have been developed to solve the layout problem at two different hierarchical levels. First, the layout of the different cells of the manufacturing system is to be determined. Second, and once the coordinate of each cell is determined, the recursive layout of each cell, in which a cell’s machine tools’ layout are to be determined, is performed. Grouping of each manufacturing
2212-8271 © 2014 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Selection and peer-review under responsibility of the International Scientific Committee of “9th CIRP ICME Conference” doi:10.1016/j.procir.2015.06.017
88
Maral Zafar Allahyari and Ahmed Azab / Procedia CIRP 33 (2015) 87 – 92
cell machine tools is to be done in a prior step using a metaheuristic; however, this part of the overall methodology is not reported in this paper. 2. Critical Literature Review The block facility layout problem that was originally formulated by Armour and Buffa (1963) [10] is concerned with finding the most efficient arrangement of m indivisible departments with unequal area requirements within facility [11]. As defined in the literature, the objective of the block layout design problem is to minimize the material handling costs inside a facility by considering the following two sets of constraints: (a) department and floor area requirements and (b) department locational restrictions; i.e., departments cannot overlap, must be placed within the facility, and some must be fixed to a location or cannot be placed in specific regions [12,13,11,14,15,16] . Cellular layout has being categorized as one of the special cases in their layout problem classification. There are increasing interests in solving the block layout problem by taking continuous approach. Heragu and Kusiak (1991) [17] presented two new models of the facility layout problem, namely linear continuous with absolute values in the objective function and constraints and linear mixed integer. A coevolutionary approach to the numerical optimization of large facility layouts is introduced by Dunker et al., (2003) [18]. Their work was based on a mixed integer model for the layout constraints and objectives. It could be concluded, a few publications seem to deal with adopting the continuous approach for CMS. Wu et al., [19,20] developed a hierarchical genetic algorithm (HGA) to concurrently solve the CF (Cellular Facility) and GL (Group Layout) problems. However, they did not consider the overlap-elimination. There are recent studies that have adopted a continuous approach [21,22]. There are some common drawbacks to the constraints such as those ones used as non-overlapping constraints for machines; and also the constraints used to force machines to stay within shop floor boundaries. Finally, Arkat et al., assumed that the machines have equal square area and cells are rectangle. However, in the real world these are poor assumption. To conclude, very few researches seem to deal with UAFLPs with a continuous representation for CMS. Moreover, the majority of the models have been too complex to be attempted using mathematical programming and Operations Research solution methods, and hence are solved using heuristics and metaheuristics, which defeat the purpose of developing a mathematical model in the first place, one could argue. In this paper, a Bilevel a novel continuous mathematical formulation has been developed by taking into account the possibility of having both vertical and horizontal aisles in the shop, as well as constraints to eliminate overlap between machine tools, as well as that between manufacturing cells. The model has been later on being solved using commercial OR algebraic modeling language and solvers. A case study finally has been employed to demonstrate the working of the model and for validation.
3. Methodology A two-tier mixed integer non-linear programming model has been developed to solve the intra-cell and inter-cell layout sequentially at two different hierarchical levels, namely at the cellular and shop floor levels. Firstly, the intracell layout is being solved; secondly, and after the layout for all manufacturing cells have been finalized, the overall approach for the whole is being solved. Hence, the leader problem is the layout at the cell level (i.e. intracellular), while follower is the layout for the whole shop (i.e. intercellular). It is important to note that initially when we run the FLP for each cell (leader problem), we define an upper limit for the length and width of each cell using constraints named as within-cell constraint. At leader level, the traffic at intra-level is the material flow among the machines (operations already assigned to machines) located in cell, which is still to be laid out and constructed. Moreover, the material flows in the follower level are inter-travel between cells. Since the Group Formation is done in advance, it is already known which machine is assigned to which cell; i.e., operations of part j processed in cell k are known ahead of time. Therefore, material flows between cells are actually the flows among the operations of parts on machines done in each cell. From computational and optimization points of view, it is important to note that dividing and conquering the FLP for CMS does not produce an optimal solution; i.e., the solution obtained would not really be the same exact global optimum solving the problem combined in one math model for the two different levels (that is if we are to assume the nonlinear model to be presented in this section is linearized). However, it is important to pinpoint that such models in the literature were complex enough that they were not really being attempted and solved for optimality using OR (Operations Research) exact methods and commercial OR software. Moreover, few of these models carried constraints that were formulated in a way that hindered the ability to solve them using these tools. To elaborate, one of the models had conditions on the decision variables associated with the overlap elimination constraints. Finally, some models even went further and overcomplicated the problem by introducing other elements such as the grouping and clustering that is needed ahead of time for cell formation, as well as the production scheduling of each cell. In our case, we find it far more efficient to go about the clustering problem beforehand using a robust developed simulated annealing metaheuristic. 4. Problem Statement The problem is to arrange facilities that are cells in follower problem and machine tools in leader problem in the continual planar site. The site has rectangle shape with specified length (L) and width (W). Moreover, there is an aisle in the site by the same length as of site, however with two different vertical dimensions ܻ௦ௗ௧ and ܻ௦௪ௗ௧ . Aisle divides the site to two sections, upper and lower. No facilities could be arranged in this area. The objective is minimizing total travel-flow cost by considering shape, size and geometric characteristics
89
Maral Zafar Allahyari and Ahmed Azab / Procedia CIRP 33 (2015) 87 – 92
constraints. Each facility has rectangle shape denotes its position by the coordinates of its centre and its predetermined length and width. Hence, the facilities consider as rigid blocks. Facilities are not allowed to overlap each other and have to be assigned in their related boundaries area, which is the site’s boundaries for lower level problem and of machines tools is the boundaries of their related cell. The traditional Cartesian Coordinate System used shown in Figure 1, represents the scheme used in this paper. The distance measurement used in this paper is Manhattan distance. The problem is formulated under the following assumptions: 1. 2. 3. 4. 5. 6.
GF is known in advanced. Machines are not in the same size. Machines must be located within a given area. Machines must not overlap with each other. Cell’s dimensions and orientation are predetermined. Each part type has a number of operations that must be processed based on its operation sequence readily available from the route sheet of parts. 7. The demand for each part type in known and is constant 8. Material handling devices moving the part between machines are assumed to carry only one part at a time 9. Inter and intra-cell movements related to the part types have different costs is related to the distance travelled. We assume that the rectangular distance between each pair of machines’ centroid. 10. In determining machine size and dimensions, the workspace required for operator usage and that needed to enforce between the different machines have been taken into account. Sets: ܲ ൌ ሼͳǡʹǡ͵ǡ ǥ ǡ ܲሽ index set of part types ܯൌ ሼͳǡʹǡ͵ǡ ǥ ǡ ܯሽ index set of machine types ܥൌ ሼͳǡʹǡ͵ǡ ǥ ǡ ܥሽ index set of cell types ܱ ൌ ሼͳǡʹǡ͵ǡ ǥ ǡ ሽ index set of operations indices for part p Parameters L: Horizontal dimension of shop floor W: vertical dimension of shop floor ܻ௦ௗ௧ : vertical dimension of upper side of aisle ܻ௦௪ௗ௧ : vertical dimension of lower side of aisle ܺ௦௧ௗ௧ : horizontal dimension of left side of aisle ܺ௦ோ௧ௗ௧ : horizontal dimension of right side of aisle ݈ : length of machine i ݓ : width of machine i ݈ : length of cell c ݓ : width of cell c ܣܥ : intracellular transfer unit cost for part j ܧܥ : intercellular transfer unit cost for part j ܦ : demand quantity for part j ܷ : 1, if operation o of part p can be done by machine i, otherwise 0 ܷ : 1, if operation o of part j can be done by machine which is located in cell c, otherwise 0 ܳ : 1, if machine i is assigned in cell c
Decision variable ݔ : X coordinate of machine i (Machine is denoted by its centroid) ݕ : Y coordinate of machine i (Machine is denoted by its centroid) ݔ : X coordinate of cell c (Cell is denoted by its centroid) ݕ : Y coordinate of cell c (Cell is denoted by its centroid) ܼ௨ : 1, if machine u is arranged in the same horizontal level as machine i, and 0 otherwise ܹƴ : 1, if cell ܿ is arranged in the same horizontal level as cell ܿƴ and 0 otherwise ܼ : 1, if cell ܿ is arranged in out of aisle horizontal boundaries and 0 otherwise ܹ : 1, if cell ܿ is arranged in out of aisle vertical boundaries and 0 otherwise The continuous bi-level programming problem is defined as: The intra-cell layout mathematical formulation to layout the different machines (machines here are the facilities) of every cell c at a time is as follows: ିଵ ெ σǡ௨ୀଵ ܷ ܷାଵ௨ ሺȁݔ െ ݔ௨ ȁ ȁݕ െ ݊݅ܯσୀଵ σୀଵ ஷ௨
ݕ௨ ȁሻܣܥ ܦ s.t. ȁݔ െ ݔ௨ ȁ ܼ௨ ሺ݈ ݈௨ ሻȀʹ ݅ǡ ݑൌ ͳǡ Ǥ Ǥ ǡ ܯ ȁݕ െ ݕ௨ ȁ ሺͳ െ ܼ௨ ሻሺݓ ݓ௨ ሻȀʹ ݅ǡ ݑൌ ͳǡ Ǥ Ǥ ǡ ܯ ݅ ൌ ͳǡ Ǥ Ǥ ǡ ܯ ݔ ܿܮ ଶ
(1) (2) (3) (4)
ݔ െ Ͳ
݅ ൌ ͳǡ Ǥ Ǥ ǡ ܯ
(5)
ݕ
݅ ൌ ͳǡ Ǥ Ǥ ǡ ܯ
(6)
ଶ ௪ ଶ ௪
ܹܿ
ݕ െ Ͳ ݅ ൌ ͳǡ Ǥ Ǥ ǡ ܯ (7) ଶ ݔ ,ݕ are integer, ܼ௨ are binary ݅ ൌ ͳǡ Ǥ Ǥ ǡ ܯ (8) Finally, the intercellular layout problem tries to layout the different cells (cells here are the facilities) of the entire shop floor is as follows: ିଵ σǡƴୀଵ ܷ ܷାଵƴ ሺȁݔ െ ݔƴ ȁ ݊݅ܯσୀଵ σୀଵ ஷƴ
(9)
ȁݕ െ ݕƴ ȁሻ ܧܥ ܦ s.t ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ ݔ ܮ ݔ െ ݕ
ଶ
ଶ ௪ ଶ ௪
(10)
Ͳ
ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ
(11)
ܹ
ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ
(12)
ݕ െ Ͳ ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ ଶ ȁݔ െ ݔƴ ȁ ܹƴ ሺ݈ ݈ƴ ሻȀʹ ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ ȁݕ െ ݕƴ ȁ ሺͳ െ ܹƴ ሻሺݓ ݓƴ ሻȀʹ ܿ ൌ ͳǡ Ǥ Ǥ ǡ ܥ Aisle Constraints: Horizontal Aisles: ሺݕ ݓ Ȁʹሻ െ ܻ௦௪ௗ௧ Mܼ ܻ௦ௗ௧ െ ሺݕ െ ݓ Ȁʹሻ Mሺͳ െ ܼ ሻ Vertical Aisles: ሺݔ ݈ Ȁʹሻ െ ܺ௦௧ௗ௧ ܹܯ
(13) (14) (15) (16) (17) (18)
ܺ௦ோ௧ௗ௧ െ ቀݔ െ ቁ ܯሺͳ െ ܹ ሻ (19) ଶ ݔ ǡ ݕ are integer, ܹƴ , ܼ , ܹ are binary ܿ ൌ ͳǡ Ǥ Ǥ ǡ (ܥ20) Equations 1 and 9 are the objective function of upper-level and lower level program respectively. The first objective function minimizes the total intracelll transportation cost of parts and the second one minimize the intercelll transportation cost of parts. Equations 4 to 7 and set of equations 10 to 13 are within-site constraints that ensure each machine tool and
90
Maral Zafar Allahyari and Ahmed Azab / Procedia CIRP 33 (2015) 87 – 92
cell are assigned within the boundaries of its corresponding cell and shop floor respectively. Equations 2 and 3; and 14 and 15 forces overlap elimination for machine tools and cell respectively. Equations 16 to 19 in leader problem ensure that no cells would be assigned in aisle boundaries. Equations 8 and 20 specify that the decision variables are binary and integer variables.
distributes metalworking tools. The company is dedicated to develop specialized Carbide, PCD (Polycrystalline diamond) and CBN (Cubic Boron Nitride) inserts, as well as multitask tooling for the aerospace, automotive and mold-die industries. The current layout of the company has is far from efficient and optimal. The facility at the company was classified as a CMS. Group formation was performed, and machine tools were assigned to their respective cells, which followed a product layout. There are five different kinds of family of cutting insert tools produced. Each part has specific monthly demand. Since some of machine tools have more than one unit, the demand is being shared among the units of those machines. The main operation done on inserts is grinding done by different grinding machine tools. Some of the machine tools have identical copies on the shop floor to increase productivity. Moreover, there are three workstations such as inspection, wash, and packaging. The operations sequence for each cutting insert tool is different from others. In other words, all the operations are not being done for each part. The list of operations of each inserts and the machine tools using for those operation are shown in table 1.
݈ƴ ݓƴ
Cell ݔƴ
L
ݕƴ
Aisle
ݓ୧ୱ୪ୣ୮୮ୣ୰େ୭୭୰ୢ୧୬ୟ୲ୣ ݈
Cell
ݓ
ݓ୧ୱ୪ୣ୪୭୵ୣ୰େ୭୭୰ୢ୧୬ୟ୲ୣ
W
Figure 1: The scheme of shop floor
5. Case Study: The case study is a Carbide Tool Inc manufactures and Table 1. Machine tools characterisations and parts’ operations sequence ID Machine Dimension X
Y
M1
Double disk (1)
12.67
5
M2
Blanchard (2)
6
9.07
M3
Wendt (3)
8.5
6.1
6.8
9.45
M4
Polish (1)
6
5
M5
EWAG (1)
4.3
7.3
M6
Surface grinding (2)
7
6
M7
Surface grinding (1)
6
7.54
M8
Swing fixture (1)
8
Dog Bone
O1
Triangular
O1
O2
O2
6
O2
O4
O3
O5
V-bottom (1)
7
6
Wire-cutting (2)
7.8
6.7
7.8
5.7
Diamond-type 1
Diamond-type 2
O1
O1
O1
O1
O4
O2
O2
O2
O7 O3
M10
Top Notch
Diamondtype 3
O3 O4
M9
Cutting Insert tools
S Shape
O5
O3
O3 O3
O3
O4
O4
M11
Laser M/C (1)
7.6
9.74
M12
Brazing (1)
4
1.8
M13
ETCH (1)
3
4
ST1
Inspection (1)
4
3
O6
O6
O7
O5
O9
O8
O4
ST2
Wash (1)
5
3
O7
O7
O8
O6
O10
O9
O5
ST3
Packing (1)
18
9
Part Demand
O6 O5
O5
O6
O4
O6
O5
O1
O8
O7
O3
O8
O8
O9
O7
O11
O10
O6
200
400
500
250
200
200
200
*The number of units for each machine tools shown in bracket.
The company’s shop floor is a rectangular shape. The current layout they have is process layout. There is no special material handling device for transforming unfinished products among machine tools. A tray with capacity of 250 pieces arranged for inserts transformation. By considering table 1 would be obvious that the number of operations done on each
part inserts tool is large enough, hence the amount of travelling taking place every day on their floor. Additionally, all the raw materials are carried from the back side of the floor and transformed all the way to the front for starting the operations. This makes extra movement which leads to extra material handling cost. The inspection site and shipping department were not properly positioned with the current
91
Maral Zafar Allahyari and Ahmed Azab / Procedia CIRP 33 (2015) 87 – 92
layout. It is located at the end of hall and washing at the front. These two operations are the last operations have to be done for all parts. Transforming all finished part to washing site also is time and cost consuming. All the parts need the combination of these three operations: surface grinding, top and bottom grinding and periphery grinding that can be done by Double Disk, Blanchard and Wendt. Since they are using process layout, Wendt machines group are located in upper Table 2. GF results No. Name
side of hall and Blanchard and Double Disk machines arranged in the lower side. Therefore, it could be concluded there is much unnecessary traffic that is primarily due to following a non-convenient layout scheme, aside from all other factors listed that were not really optimized. GF done in advance, there are 4 cells with specific types of machine tools. The characterization of GF is as follows:
Machine tools / Work Station
1
Primary Cell
Double Disc (1)
Blanchard (2)
Polish (1)
2
Grinding Cell
Surface Grinding (2)
Swing Fixture (2)
V-Bottom (1)
Wendt (3)
3
Diamond Cell
Wire-cutting (2)
Surface Grinding (1)
EWAG (1)
Brazing (1)
4
Final Cell
ETCH (1)
Inspection (1)
Wash (1)
Packing (1)
Laser M/c (1)
*The number of units for each machine tools shown in bracket.
FICO Xpress Optimization Suit Software has been used to solve the continuous formulation of this paper. Since the mathematical formulation is nonlinear both Successive Linear Programing (SLP) and Non-linear Programming (NLP) solver have been used. 5.1. Intracellular Layout: For the leader problem the layout of the different machine tools and work stations in their respective cells are being solved. The intracellular travel cost per unit distance per one unit of each respective part are ¢10, ¢20, ¢20, ¢25, ¢30 respectively for Dog Bone, S Shape, Triangular, Top Notch and all types of Diamond. The results of intracellular layout represented in table 3 and the sketch of layout plan in each cell shown in figure 2. Table 3. Intercellular costs Cell’s Name
Dimension
Material Handling Cost
Length
Width
Primary Cell
35
25
¢737,175
Grinding Cell
26
30
¢151,625
Diamond Cell
30
16
¢659,100
Final Cell
34
18
¢862,875
Figure 2. Sketch of intercellular and intracellular layout plan
6. Conclusion and Future Work:
5.2. Intercellular Layout: In the follower problem, the four cells with the exact dimensions have been assigned in the whole shop floor with 97×60 feet length and width. The intercellular travel cost per unit distance per one unit of each respective part Dog Bone, S Shape, Triangular, Top Notch, all types of Diamond among cells are ¢10, ¢20, ¢40, ¢10, ¢15 respectively. The traveling flows costs among cells are ¢852617. The final sketch which is showing intercellular and intracellular layout has been shown as follows:
CMS is one of the most important facility layout problems, which does include group formation, intracellular layout, and intercellular layout. The focus of this paper is on the latter two aspects of the problem. The bi-level approach taken in this paper is to optimize the intracellular and intercellular facility layout problem in a sequential manner by solving a leader and follower problem at the cell- and shop floor-levels respectively. A continuous mixed-integer nonlinear programming mathematical formulation has been developed for each of the two levels. A novel overlap elimination constraint has been formulated. For the follower problem, vertical and horizontal aisle constraints have been developed. The develop model has been demonstrated and verified using a realistic industrial case study. The group formation has been done before and the results presented in this paper. Based on this results there are 4 main cells by special kind of machine tool. The results out of the model in this case study is intuitive. Generally speaking, solving both leader and follower problem, the layout generated was optimized in a sense to minimize travel-flow cost. It was found out that that for each of the facilities layout problems being solved at both the leader- and follower-level, the facilities have been lamped
92
Maral Zafar Allahyari and Ahmed Azab / Procedia CIRP 33 (2015) 87 – 92
together in a way to minimize the travel in between, and hence the travel- flow cost. For future work, the developed model will be linearized and solved for optimality.
7. References [1] Shayan, E. & Chittilappilly, A. Genetic algorithm for facilities layout problems based on slicing tree structure. International Journal of Production Research 2004; 42: 405567. [2] Kia, R., Baboli, A., Javadian, N., Tavakkoli-Moghaddam, R., Kazemi, M. & Khorrami, J. Solving a group layout design model of a dynamic cellular manufacturing system with alternative process routings, lot splitting and flexible reconfiguration by simulated annealing. Computers & Operations Research 2012; 39: 2642-58. [3] Tavakkoli-Moghaddam, R., Javadian, N., Javadi, B. & Safaei, N. Design of a facility layout problem in cellular manufacturing systems with stochastic demands. Applied Mathematics and Computation 2007; 184: 721-8. [4] Ariafar, S., Ismail, N., Ariffin, M. K. A. M., Tang, S. H. & Firoozi, Z. Illustration and model verification for an inter-cell and intra-cell facility layout problem in a cellular manufacturing system. Applied Mechanics and Materials, 2013; v263-266, 1310-14. [5] Drira, A., Pierreval, H. & Hajri-Gabouj, S. Facility layout problems: a survey. Annual Reviews in Control 2007; 31: 255-67. [6] Das, S. K. A facility layout method for flexible manufacturing systems. International Journal of Production Research 1993; 31: 279-97. [7] Wemmerlov, U. & Johnson, D. J. Cellular manufacturing at 46 user plants: implementation experiences and performance improvements. International Journal of Production Research, 35: 29-49. [8] Ariafar, S., Ismail, N., Tang, S. H., Ariffin, M. K. A. M. & Firoozi, Z. Design of a facility layout model in hybrid cellular manufacturing systems under variable demand. International Journal of Industrial and Systems Engineering 2011; 9: 373387. [9] Defersha, F. M. & Chen, M. A comprehensive mathematical model for the design of cellular manufacturing systems. International Journal of Production Economics 2006; 103, 767-783. [10] Armour ,G.C. & Buffa E.S. A heuristic algorithm and simulation approach to relative location of facilities, Management Science 1963 ; 9 (2): 294–309. [11] Meller, R. D. & Kai-Yin, G. The facility layout problem: recent and emerging trends and perspectives. Journal of Manufacturing Systems 1996; 15: 351-66. [12] Castillo, I., Westerlund, J., Emet, S. & Westerlund, T. Optimization of block layout design problems with unequal areas: A comparison of MILP and MINLP optimization methods. Computers and Chemical Engineering 2005; 30: 5469. [13] Kusiak, A. & Heragu, S. S. The facility layout problem. European Journal of Operational Research 1987; 29, 229-51.
[14] Chiang, W.-C. Visual facility layout design system. International Journal of Production Research 2001; 39, 18111836. [15]Yaman, R., Gethin, D. T. & Clarke, M. J. An effective sorting method for facility layout construction. International Journal of Production Research 1993; 31: 413-27. [16] Ketcham, R. L., Malstrom, R. M. & Mcroberts, K. L. A comparison of three computer assisted facilities design algorithms. Computers & Industrial Engineering 1989; 16: 375-86. [17] Heragu, S. S. & Kusiak, A. Efficient models for the facility layout problem. European Journal of Operational Research 1991; 53: 1-13. [18] Dunker, T., Radons, G. & Westkamper, E. A coevolutionary algorithm for a facility layout problem. International Journal of Production Research 2003; 41, 3479500. [19] Wu, X., Chu, C.-H., Wang, Y. & Yan, W. 2007. A genetic algorithm for cellular manufacturing design and layout. European Journal of Operational Research 1997; 181: 156-167. [20] Wu, X., Chu, C. H., Wang, Y. & Yan, W. Concurrent design of cellular manufacturing systems: A genetic algorithm approach. International Journal of Production Research 2006; 44: 1217-1241. [21] Arkat, J., Farahani, M. H. & Ahmadizar, F. Multiobjective genetic algorithm for cell formation problem considering cellular layout and operations scheduling. International Journal of Computer Integrated Manufacturing 2012a; 25: 625-635. [22] Arkat, J., Hosseinabadi Farahani, M. & Hosseini, L. Integrating cell formation with cellular layout and operations scheduling. International Journal of Advanced Manufacturing Technology 2012; 61: 637-47.