Application of Layered Coding Genetic Algorithm in Optimization of Unequal Area Production Facilities Layout

Hou, Shiwang; Wen, Haijun; Feng, Shunxiao; Wang, Hui; Li, Zhibin

doi:https://doi.org/10.1155/2019/3650923

Computational Intelligence and Neuroscience

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 3650923 | https://doi.org/10.1155/2019/3650923

Application of Layered Coding Genetic Algorithm in Optimization of Unequal Area Production Facilities Layout

Shiwang Hou,^1,2Haijun Wen,³Shunxiao Feng,³Hui Wang,³and Zhibin Li³

Academic Editor: Maciej Lawrynczuk

Received08 Apr 2019

Revised25 May 2019

Accepted28 May 2019

Published19 Jun 2019

Abstract

Unequal area facilities layout problem (UA-FLP) is an inevitable problem in the process of new construction, reconstruction, and expansion of enterprises. The rationality of the facilities layout has a great influence on the operation performance of the production system. Finding the optimal solution of UA-FLP according to the requirement of production process is the main content of the plant design. The facilities were constrained by given areas and aspect ratio, respectively. By adopting the method of slicing tree, the layout space was divided into multiple regions for each facility. The genetic algorithm was developed by using layered coding to show the slicing process. Considering the production logistics cost as well as the adjacency relations between the facilities, the goal function was established and the optimal solution was obtained by running the proposed algorithm. Finally, the feasibility of the proposed approach was validated by a set of known problems. The comparison results show that it can provide decision support for rapid optimal layout of multifacilities.

1. Introduction

UA-FLP was proposed originally by Armour and Buffa in [1], and its objective is to determine the good locations for a given set of departments with different areas on some workshop floor to optimize the material handling cost and/or other objectives. The facility can be small or big according to the level of the facility layout, but it should be a physical entity with some function. There is about 20% to 50% of the processing cost used for material handling, and scientific and reasonable facilities layout can save at least 10% to 30% of the material handling fee [2, 3].

Classic FLPs tend to study equal-area facilities arrangement; i.e., all facilities have the same area and shape. In this case, the facility centroid is fixed and the overall closeness or distance between facilities will not change when switching the location of any two facilities. But, the equal-area hypothesis is impractical, and the facility area is often unequal. So, the centroid of each facility depends on its area and shape and has no regular distribution as equal-area case.

Heuristic algorithm based on the discrete model is proposed in reference [4] to deal with the UA-FLP. This approach divided the district into some small squares with a fixed area, and each facility was allocated some numbers of squares most close to its area by use of heuristic algorithm. As shown in Figure 1, the grid size determines the precision of the facility representation; the smaller the grid size, the more precise the facility representation. The facility shape was also represented well by this approach. But, the time consumption is great to calculate the interference between facilities during the detailed layout process. Also, this approach is easy to produce facility layout with irregular shape. By far, the most commonly used UA-FLP model is mainly based on block diagram with unequal area as depicted in [5].

Many researchers have studied the optimization methods for UA-FLP. The literature about the methods can be mainly divided into three categories. The first one is deterministic algorithm to calculate accurate solution, such as mathematical programming [6, 7], mixed-integer linear programming [8], and QAP as aforementioned. The second category is heuristic algorithm mentioned in reference [9–11]. Another category is intelligent algorithm, such as genetic algorithm (GA) [12, 13], particle swarm optimization (PSO) [14, 15], ant colony algorithm (ACA) [16, 17], simulated annealing algorithm (SAA) [18, 19], and tabu search algorithm [20].

Deterministic algorithm has high requirement in memory and CPU time, so it is often used to solve small-scale UA-FLP and reference [21] proved its effectiveness. Heuristic algorithm and intelligent algorithms have faster calculation speed and are suitable for large-scale UA-FLP. Furthermore, genetic algorithm is most widely used, and different genetic algorithms were developed for various UA-FLPs.

A genetic algorithm hybridized with local search to obtain the Pareto solutions set was proposed in [22], and the author adapted a random weight to combine values of two objectives. For unequal area facility layout problems, a genetic algorithm based on slicing structure was developed in [23], and four objective functions such as material handling costs, aspect ratio, closeness, and distance requests were considered simultaneously by use of a Pareto-based evolutionary approach. In order to improve the performance of premature convergence, lack of diversity, and high computational cost, an island model genetic algorithm was proposed in [24]. The compared results showed the proposed approach has great improvement on the above aspects. In [25], a multiobjective interactive genetic algorithm was proposed by considering both quantitative aspects and subjective features, which allowed the interaction between the expert designer and the algorithm. A biased random-key genetic algorithm to determine the placement order and dimensions of each facility was proposed in [26], and the results showed its better performance for 19 of the 28 benchmark facility layout problems.

For large-scale UA-FLP problems, a genetic algorithm combined with a decomposition strategy was proposed in [27]. Compared with basic genetic algorithm, the experiments in the paper showed that the proposed approach had an average solution improvement of 6% or 7% for large-scale instances with 90 or 100 facilities.

This paper put forward a method of LCGA (layered coding genetic algorithm) for slicing-based plane splitting to lay out facilities. This approach can generate feasible solution rapidly with the help of the layered coding method. It can provide larger-scale UA-FLP solving a new thinking. The rest of this paper is organized as follows: Section 2 provides a description of the UA-FLP dealt with in this paper. The design process of the proposed layered coding genetic algorithm is presented in Section 3. Section 4 compares the proposed approach with some known UA-FLPs. The conclusions are presented in Section 5.

2. UA-FLP Description

2.1. Location Relations of UA-FLP

Assume that all facilities are rectangular blocks with a given area. Considering the building module, facilities need to satisfy a given aspect ratio constraint to avoid producing an approach with long narrow shape facilities layout. All facilities must be located in a given area and cannot overlap between them. Take a UA-FLP involving facilities i and j; for example, the relative position relations are as shown in Figure 2.

Figure 3(a) shows the projection polygon of facility space. Figure 3(b) shows the projection polygon of facility and its necessary space, including transport corridor, operational space, and maintenance space for workers. To facilitate the facilities layout, Figure 3(c) considers the rectangular envelope of Figure 3(b) as the objects of UA-FLP. So, the necessary horizontal and vertical spacing between facilities shown in Figure 2 can be decomposed to the corresponding facility (as shown in Figure 4). By doing so, the amount of computation to solve UA-FLP can be decreased without losing the accuracy of result.

(a)

(b)

(c)

2.2. Objective Function and Constraints

Satisfying the abovementioned basic hypothesis, the UA-FLP becomes the following optimization problem: arranging n facilities to a specified area to minimize the material handling cost and maximize the closeness scores.

2.2.1. Minimizing the Material Handling Cost

where n is the number of facilities, c_ij is the per unit handling cost between facility i and facility j, f_ij is the logistics quantity between facility i and facility j, and d_ij is the distance between facility i and j.

2.2.2. Maximizing the Closeness Scores

where n is the number of facilities and is the ratio of distance between facility i and facility j to the maximum distance between facilities in a given layout approach, and is used to represent their closeness factor between two facilities. r_ij is the quantitative score of the nonlogistics relationship level determined by the SLP method (Table 1).

This problem belongs to multiobjective optimization problem (MOOP) since there is more than one objective function to be optimized simultaneously. Optimal decisions need to be taken in the presence of tradeoffs between two or more conflicting objectives. Taking UA-FLP as example, material handling cost is minimized while closeness scores are maximized while locating facilities. For a nontrivial multiobjective optimization problem, no single solution exists that simultaneously optimizes each objective. There exist a (possibly infinite) number of Pareto optimal solutions; i.e., under these solutions, none of the objective functions can be improved in value without degrading some of the other objective values. All these Pareto optimal solutions are considered equally good if there is no additional subjective preference information.

There are many kinds of forms of solutions for different goals, such as a representative set of Pareto optimal solutions and/or a single solution that satisfies the subjective preferences of decision maker. In this paper, we adopt the latter in order to obtain a well-determined layout approach.

Let M be the maximum of r_ij and and be the weight coefficients of material handling cost and nonlogistics closeness, respectively. The above two goals can be synthesized into a minimizing objective function:subject to where , , , and denote the x-coordinate and y-coordinate of upper-right corner and lower-left corner of facility i, respectively. Formulas (4) and (5) ensure the aspect ratio of facility i in a given range , and formula (6) constrain the area of facility i equal to given value .

Manhattan distance is adapted to calculate the distance between facilities, namely:where and denote the x-coordinate and y-coordinate of the centroid. The exact coordinate position of facility i, , can be obtained by the following formula:

3. Algorithm Design

3.1. Basic Layout Principle for UA-FLP Using Slicing Tree

In order to locate all facilities into a given area, the area should be divided into subareas with the same number of facilities. This paper adopted plane segmentation method to generate slicing tree in order to describe the relative position relationships of facilities during locating them.

For a UA-FLP with n facilities, the slicing tree contains n leaf nodes and (n − 1) internal nodes. The information of the plane segmentation mode is contained in internal nodes, i.e., horizontal split (labeled H) or vertical split (labeled V). Taking a UA-FLP with 5 facilities as an example, a feasible plane segmentation approach is shown in Figure 5.

The process of plane segmenting is the process of facilities layout. The change of coordinates for each partition after each plane segmenting is described as follows. The area has only two points: and , when there is no facility located. The coordinate of upper-right corner is marked with the maximum limit value; i.e., x-coordinate is the sum of width of all facilities supposing they are side-by-side arranged horizontally, and y-coordinate is the sum of length of all facilities supposing they are side-by-side arranged vertically. Taking the first partition as an example, the splitting process is as follows. Suppose AB is the first cutting line. The plane area is divided into two partitions. The upper part is for facilities 1 and 2, and the lower part is for facilities 3 to 5. After the first segmentation, the coordinates of point A and B are , .

Similarly, the other four divisions can be done. Finally, all facilities are located in different parts, and the precise coordinates of the facilities are as follows: facility 1{A, ()}, facility 2{F, ()}, facility 3{O, ()}, facility 4{G, ()}, and facility 5{D, ()}. So, a feasible layout approach is provided. Based on this layout idea, this paper studies the hierarchical coding genetic algorithm in order to realize the layout scheme iterative optimization.

3.2. Genetic Algorithm Design

For a UA-FLP with n facilities, there are n! possible combinations of their position and the number of combination will be larger if the shape or orientation of each facility is considered. Above all, there is a large number of local optimum in this huge solution space. So, it is a NP-hard problem, and it is advisable to find a suboptimal solution in acceptable cost (money, time, or computing resource) by use of some heuristics search algorithms.

GA is a bionic algorithm for searching the optimal solution based on the principle of biological evolution. It simulates the natural process of gene recombination and evolution and compiles the parameters into binary-code or decimal-code (or other codes) genes. Several genes constitute a chromosome (individual), and many chromosomes carry out operations similar to natural selection, pairing crossover, and mutation. The final optimization result is obtained after repeated iterations (that is, generation inheritance).

3.2.1. Layered Coding Approach

Coding is to map the phenotype data in solution space into genotype data in genetic structure. During iterations of GA, a coding string represents a solution and genetic operations are done by operating the bits of this string. So, the coding method also affects the genetic operators.

There are mainly two coding methods: real number coding and binary coding. The former uses a real number as a gene, is easy to understand, and does not need decoding process, but it is also easy for premature convergence, thus falling into local optimum. The latter uses a binary string with specific length as a gene and has higher stability, larger population diversity, and better performance for global search. In this paper, we adopted the binary coding method. The number of bit is determined by the accuracy of the solution to be achieved.

For example, suppose an x-coordinate ranging in [0, 4] and the solution is exactly 4 decimal places behind the decimal point. The solution space is divided into (0-1) (1e + 4) = 10,000 equal fractions. It takes 14 bits of binary to represent a solution; i.e., the coding of a solution is a 14 bit binary string since 2¹³ < 10000 < 2¹⁴. The decoding process is as follows:

Generally, for , the value of x after decoding is

The coded chromosome string should represent the following information simultaneously: facility sequence for layout, splitting point sequence, and splitting mode. The coding approach will be detailed below.

Facility sequence code is in the first layer. N facilities are coded by n different integers in the interval [1, n] by use of integer coding. The coding string can be any sequences of n integers in the interval [1, n] to allow different facility locating orders. Figure 6 shows a code string [1–13] of 13 facilities, f1, …, f₁₃, with a random order of integers from 1 to 13, and the value of each bit of code string denotes the number of the facility it represents.

Splitting point code lies in the second layer. A feasible splitting point lies between every two bits of facility coding string. For every two adjacent bits of facilities coding string, there is a splitting point. For a n-bit facility coding string, there are n − 1 splitting points. We denote every splitting position as an integer between 1 to n − 1. For the facility code string in Figure 6, its corresponding splitting position is shown in Figure 7.

The facility set is divided into two different parts at the splitting point selected firstly. The remainder splitting points are contained in these two subsets of facilities. After each splitting operation, the number of subsets of facilities increases by 1. After n − 1 splitting, the facility set with n facilities will be divided into n single facilities locating in n different area blocks. For the code string of splitting point, we code them as a sequence from 1 to (n − 1) with random order corresponding to different plane segmentation approaches. The value of each bit of code string denotes the number of the splitting position in the facility code string. Taking the string shown in Figure 6 as an example, there are 12 positions that can be set as splitting points when carrying out plane segmentation. Suppose we produce a splitting point code string as [1–12], the splitting operation will begin with the 7^th splitting position, then 4^th, and finally, 6th. The splitting process is shown in Figure 8.

The last layer of coding provides the information of splitting mode, horizontally or vertically. The splitting result of these two ways is different, and the facility layout is also different. We use 0 for horizontal and 1 for vertical. So, the splitting mode code is a binary string that has the same number of bits as the splitting point code. For the example as shown in Figure 8, assume a splitting mode code is [0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0]. Figure 9 shows a complete three-layer coding string or a chromosome of a feasible solution for a UA-FLP with 13 facilities. For the first step, the splitting mode is horizontal and the process of plane segmentation is similar to the description in Section 3.1.

So, for a UA-FLP with n facilities, the three-layer coding string can be expressed as shown in Figure 10.

3.2.2. Crossover Operation

In genetic algorithms and evolutionary computation, crossover, also called recombination, is a genetic operator used to combine the genetic information of two parents to generate new offspring. It is one way to stochastically generate new solutions from an existing population. Newly generated solutions are typically mutated before being added to the population.

The coding method determines the data structures to store genetic information and also affects the crossover operators.

Due to the aforementioned layered coding structure, genetic operation must be carried out by the segment to ensure the feasibility of the new code string. Taking the chromosome gene string in Figure 7 as an example, crossover points are selected in three layers, respectively. The two parent individuals swap the gene segments before and after the crossover points in three layers, respectively.

A crossover operation of layered coding string of 5 facilities is shown in Figure 11.

It can be found that some bit value of new code strings will lose or repeat in the first two layers after crossover operation. So, repair operation is necessary. Sort the missing value in ascending order and replace the repeated value of corresponding layer code strings. Taking the above offspring as example, the layer 1 code of offspring is [5 3 1 5 1]. The reappeared values are 5 and 1, and the missing values are 2 and 4. We repair the offspring code of layer 1 as [5 3 1 2 4]. After similar repairment, the offspring code of layer 2 is [4 1 2 3].

3.2.3. Mutation

Mutation changes one or more gene values in a chromosome from its initial state in order to maintain genetic diversity from one generation of a population to the next. This can also prevent the population of chromosomes from becoming too similar to each other, thus slowing or even stopping evolution. For the layered code, we also adopted different mutation operators.

A mutation operator involves a probability P_m that an arbitrary bit in a genetic sequence will mutate from its original state. A common method of implementing the mutation operator involves generating a standard uniform-distributed random number R for each bit in a sequence.

The corresponding bit will be modified if . This single-point mutation is suitable for third-layer binary code of an individual.

For the first- and second-layer code, mutation operation is implemented by interchanging two genes of code string in order to ensure the feasibility of new individual after mutation. The gene bits to mutate are selected randomly and also use a standard uniform-distributed random number R to determine whether or not the selected bits will be swapped.

3.3. Processing for UA-FLP with Empty Space

The method proposed herein split the area into subareas that are consistent with the number of facilities, each of which accommodates a facility. For a UA-FLP with empty space, we firstly need to turn it into a UA-FLP without empty space. The detailed process is explained as follows.

Suppose the ratio of the original region for all the facilities is and the area of i-th facility is A_i, we can arrange all the facilities in a new region, marked by the red border in Figure 12, with the area equal to the sum of all facilities since all necessary spaces for each facility has been included in its area as illustrated in Figure 3. The width W′ and height H′ of this new region should satisfy the following equations:

In general, solutions found in this way will be better than solutions found in other methods because the horizontal and vertical coordinates of each facility will be smaller.

4. Case Study

In order to validate the proposed approach, a set of problems described in the literature was used in this section. All the tested problems are shown in Table 2.

The algorithm is coded in Matlab 2015b. The computer’s configuration running the algorithms was as follows: Intel Core i5-4460 (3.20 GHz), 8 GB RAM, and a Windows 10 operating system.

The algorithm parameters settings for the tested problems are listed in Table 3, and the algorithm proposed in this paper is described in detail in Section 4.1.

4.1. Parameters Setting and Pseudocode of Algorithm

(1)Coding Method. Integer coding was adopted in this case. The number n of facilities was assigned to n_f. The layout approach and the splitting point sequence were denoted by two interpermutations which took values from [1, n] and [1, n − 1], respectively. The splitting mode of each splitting operation was represented by a thirteen bits encoded string, and each bit can take value 0, splitting horizontally, or 1, splitting vertically.(2)Parameters Setting. See Table 3.(3)The Generation of Initial Population. According to the above approach, initial facility population, splitting position string, and splitting mode string can be generated by use of Algorithm 1.

	Input: Popu_size, n_f
	Output: F(facility population), Sp[Popu_size × n_f] (Splitting position string), Sm[1 × n_f − 1] (splitting mode string)
(1)	for I ⟵ 1 to Popu_size
(2)	F(i,:) ⟵ permutation of n_f integers valued from 1∼n_f;
(3)	Sp(i,:) ⟵ permutation of (n_f − 1) integer valued from 1∼(n_f − 1);
(4)	Sm(i,:) ⟵ permutation of (n_f − 1) number valued 1 or 0;
(5)	end for
(4)	Establishment of fitness function.

Algorithm 2 carries out the establishment of fitness function and its value calculation.

	Input: (x, y) ⟵ coordinate of facility centroid
	D ⟵ Matrix of logistics quantity between facilities;
	C ⟵ Matrix of unit logistics cost between facilities;
	r ⟵ Matrix of closeness score between facilities;
	r_max ⟵ maximum of closeness score between facilities;
	a1 ⟵ weight of logistics cost; a2 ⟵ weight of nonlogistics factors;
	Output: goal_v (value of fitness function)
(1)	for i ⟵ 1 to n_f
(2)	for j ⟵ 1 to n_f
(3)
(4)	end for
(5)	end for
(6)	G1 ⟵ sum according to the row (sum according to the column(D.flow_f.flow_c));
(7)	G2 ⟵ sum according to the row (sum according to the column(D./D_max.relation_v));
(8)	goal_ ⟵ a1G1 + a2(r_max-G2);
(5)	Genetic manipulation

The crossover operation and the corresponding repair approach of facility level are presented in Algorithm 3.

	Input: Popu_size, Pc, n_f
	Output: newPop_f (newfacility population after crossover peration and repairment)
(1)	f_popu_cro ⟵ individuals in facility level for crossover;
(2)	so_popu_cro ⟵ individuals in splitting position level for crossover;
(3)	sp_popu_cro ⟵ individuals in splitting mode level for crossover;
(4)	newPop_f ⟵ zeros(popu_sizePc, n_f);
(5)	r1 ⟵ sorted array of r1 from lowest to highest;
(6)	while r1(1) = r1(2)
(7)	Do r1 ⟵ regenerate two random integers within [1, n_f];
(8)	end while
(9)	r1 ⟵ sorted array of r1 from lowest to highest;
(10)	for j ⟵ 1 TO popu_sizePc/2
(11)	newPop_f(2j − 1,:) ⟵ [f_popu_cro(j + popu_n0.8/2,1 : r1(1)),f_popu_cro(j, r1(1) + 1 : r1(2) − 1),
f_popu_cro(j + popu_n0.8/2, r1(2) : end)];
(12)	newPop_f(2 j,:) ⟵ [f_popu_cro(j, 1 : r1(1)), f_popu_cro(j + popu_n0.8/2, r1(1) + 1 : r1(2) − 1),
f_popu_cro(j, r1(2) : end)];
(13)	end for
(14)	ch_temp ⟵ integer order sequence taking value within [1, N_f];
(15)	for i ⟵ 1 TO popu_nPc
(16)	ind1 ⟵ locating the position of element with multiple occurrences in newPop_f;
(17)	temp_re1 ⟵ all elements that appear in newPop_f;
(18)	temp_re2 ⟵ temp_re1(ind1);
(19)	ind2 ⟵ locating the position of element in ch_temp that do not appear in newPop_f;
(20)	for j ⟵ 1 TO length(temp_re2)
(21)	ind3 ⟵ locating the position of element in ith row of newPop_f that equals to the jth
element of temp_re2;
(22)	newPop_f (i,ind3) ⟵ ch_temp(ind2(j));
(23)	end for
(24)	end for

As for the crossover and repair operation of the splitting point level, the method is similar to the above ideas; just replace the f_popu and n_f with so_popu and (n_f − 1). There is no need for repair operation in the splitting mode level, and its crossover operation is the same as the facility level.

The mutation operation of the facility level is shown in Algorithm 4.

	Input: f_popu, n_f
	Output: f_popu after mutation operation
(1)	Pm ⟵ Mutation probability;
(2)	for i ⟵ 1 TO N_f
(3)	if Pm > random in [0,1];
(4)	f_popu_mu ⟵ f_popu(i,:);
(5)	r4 ⟵ unidrnd(N_f, 1, 2);
(6)	while r4(1) = = r4(2)
(7)	Do r4 ⟵ unidrnd(N_f, 1, 2);
(8)	end while
(9)	f_popu_mu(1, r4(1)) ⟵ f_popu(i, r4(2));
(10)	f_popu_mu(1, r4(2)) ⟵ f_popu(i, r4(1)); f_popu(i,:) ⟵ f_popu_mu;
(11)	end if
(12)	end for

The mutation operation of the splitting point level or splitting mode level can be completed by replacing the f_popu and N_f of Algorithm 4 with so_popu and (n_f − 1) or sp_popu and (n_f − 1).

4.2. Computational Results and Analysis

According to the above parameters setting, all the tested problems were solved by the proposed layout algorithm; the average optimal value of objective function of each tested problems was obtained as shown in Table 4 by running the proposed algorithm 30 times for each problems.

The optimal layout approach of all tested cases obtained by the proposed approach is presented in Figures 13(a)–13(l).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 13

Optimal layout approach for the test cases. (a) Optimal layout approach for O7. (b) Optimal layout approach for O8. (c) Optimal layout approach for O9. (d) Optimal layout approach for F10. (e) Optimal layout approach for VC10. (f) Optimal layout approach for MB11. (g) Optimal layout approach for Ba12. (h) Optimal layout approach for Ba14. (i) Optimal layout approach for AB20. (j) Optimal layout approach for SC30. (k) Optimal layout approach for SC35. (l) Optimal layout approach for P62.

From Table 4, we found that the proposed approach can search better solutions than did the methods from literature for 8 of 12 tested problems. As shown, the proposed approach has more improvement when the facility number becomes bigger; for example, the improvement present can be up to 4–7% when the facility number reaches 30 and 35. Moreover, the average running time for finding the best solutions and the average total CPU running times have reduced drastically than did the approach of other literature. But, for the test case of P62 with no spare space for all the facilities, the proposed approach has a large gap from the optimal result of [35].

5. Conclusions and Prospects

Considering the constraints of area and aspect ratio of facilities, this paper proposed a slicing-tree-based binary plane segmentation method. The given area is divided into some blocks whose number equal to the number of facilities waiting for arrangement. The optimal solution was found by use of the layered coding genetic algorithm with the goal of maximizing the closeness relationship score and minimizing the material handling cost between facilities at the same time. The results of comparison of above 12 known problems between other literature methods and proposed approach show the effectiveness of the plane segmentation layout strategy and the reliability of layered coding genetic algorithm for solving the problems.

We draw the following conclusions based on the above study:(i)It is reasonable and effective to partition the facilities layout area by use of the binary plane segmentation method, which can arrange a reasonable block area for each facility. So, the feasibility of solution was guaranteed during the iterative process of genetic algorithm.(ii)By use of the layered coding genetic algorithm, the optimal splitting approach, i.e., optimal facility layout scheme, will be found during multiple iterative process. The result can provide decision support for actual production facility layout.(iii)The plane splitting process was expressed well by the layered coding approach. When the UA-FLP problem changed, i.e., the facility number, facility area, and also the aspect ratio, these changed values can be assigned in the form of parameters; as a result, the corresponding optimal layout approach can be output quickly. Also, the layout approaches of different parameters can be compared easily in order to find the key influence factors.

The optimal layout approach of UA-FLP can be obtained fast by use of the plane segmentation method and layered coding genetic algorithm; the corresponding objective function value and the minimal area needed for locating facilities can also be provided. For multiproduct facility layout problem, different layout approaches for different products can be output by running the proposed algorithm in terms of the different requirement of product for logistics cost and closeness relationship between facilities, i.e., realizing a dynamic optimization for multiproduct facility layout. The facility number has a significant impact on performance of proposed algorithm, and there is the possibility of fall into local optimum, so the large-scale facility layout optimization algorithm and the combination scheme with the local search algorithm can be chosen as the direction of further research.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by China Scholarship Fund (CSC No. 201708430274); Scientiﬁc Research Fund of Hunan Provincial Education Department (Grant No. 16B208); Humanities and Social Sciences Project of Hunan Province (Grant No. 16YBX018); and MOE (Ministry of Education in China) Project of Humanities and Social Sciences (Grant No. 13YJC630049).

References

G. C. Armour and E. S. Buffa, “A heuristic algorithm and simulation approach to relative location of facilities,” Management Science, vol. 9, no. 2, pp. 294–309, 1963.
View at: Publisher Site | Google Scholar
J. A. Tompkins, Facilities Planning, John Wiley, New York, NY, USA, 4th edition, 2010.
F. Liu, M. Dong, W. H. Hou, and Y. Zhai, “Multi criteria evaluation for facility layout problems,” Journal of Shanghai Jiaotong University, vol. 4, pp. 664–668, 2007.
View at: Google Scholar
M.-J. Wang, M. H. Hu, and M.-Y. Ku, “A solution to the unequal area facilities layout problem by genetic algorithm,” Computers in Industry, vol. 56, no. 2, pp. 207–220, 2005.
View at: Publisher Site | Google Scholar
U. Haktanirlar and K. B. S. Kulturel, “An artificial immune system based algorithm to solve unequal area facility layout problem,” Expert Systems with Applications, vol. 39, no. 5, pp. 5384–5395, 2012.
View at: Publisher Site | Google Scholar
T.-Y. Wang, H.-C. Lin, and K.-B. Wu, “An improved simulated annealing for facility layout problems in cellular manufacturing systems,” Computers & Industrial Engineering, vol. 34, no. 2, pp. 309–319, 1998.
View at: Publisher Site | Google Scholar
K. L. Mak, Y. S. Wong, and F. T. S. Chan, “A genetic algorithm for facility layout problems,” Computer Integrated Manufacturing Systems, vol. 11, no. 1-2, pp. 113–127, 1998.
View at: Publisher Site | Google Scholar
M. C. Georgiadis, G. Schilling, G. E. Rotstein, and S. Macchietto, “A general mathematical programming approach for process plant layout,” Computers & Chemical Engineering, vol. 23, no. 7, pp. 823–840, 1999.
View at: Publisher Site | Google Scholar
K. R. Kumar, G. C. Hadjinicola, and T. I. Lin, “A heuristic procedure for the single-row facility layout problem,” European Journal of Operational Research, vol. 87, no. 1, pp. 65–73, 1995.
View at: Publisher Site | Google Scholar
A. Taghavi and A. Murat, “A heuristic procedure for the integrated facility layout design and flow assignment problem,” Computers & Industrial Engineering, vol. 61, no. 1, pp. 55–63, 2011.
View at: Publisher Site | Google Scholar
R. Kothari and D. Ghosh, “Insertion based Lin-Kernighan heuristic for single row facility layout,” Computers & Operations Research, vol. 40, no. 1, pp. 129–136, 2013.
View at: Publisher Site | Google Scholar
K. Krishnakumar and S. N. Melkote, “Machining fixture layout optimization using the genetic algorithm,” International Journal of Machine Tools and Manufacture, vol. 40, no. 4, pp. 579–598, 2000.
View at: Publisher Site | Google Scholar
S. Vallapuzha, E. C. De Meter, S. Choudhuri, and R. P. Khetan, “An investigation into the use of spatial coordinates for the genetic algorithm based solution of the fixture layout optimization problem,” International Journal of Machine Tools and Manufacture, vol. 42, no. 2, pp. 265–275, 2002.
View at: Publisher Site | Google Scholar
H. Samarghandi, P. Taabayan, and F. F. Jahantigh, “A particle swarm optimization for the single row facility layout problem,” Computers & Industrial Engineering, vol. 58, no. 4, pp. 529–534, 2010.
View at: Publisher Site | Google Scholar
S. Önüt, U. R. Tuzkaya, and B. Doğaç, “A particle swarm optimization algorithm for the multiple-level warehouse layout design problem,” Computers & Industrial Engineering, vol. 54, no. 4, pp. 783–799, 2008.
View at: Publisher Site | Google Scholar
M. Solimanpur, P. Vrat, and R. Shankar, “Ant colony optimization algorithm to the inter-cell layout problem in cellular manufacturing,” European Journal of Operational Research, vol. 157, no. 3, pp. 592–606, 2004.
View at: Publisher Site | Google Scholar
Y. H. Chen, “A new data structure of solution representation in hybrid ant colony optimization for large dynamic facility layout problems,” International Journal of Production Economics, vol. 142, no. 2, pp. 362–371, 2013.
View at: Publisher Site | Google Scholar
R. Matai, “Solving multi objective facility layout problem by modified simulated annealing,” Applied Mathematics and Computation, vol. 261, pp. 302–311, 2015.
View at: Publisher Site | Google Scholar
R. Şahin, “A simulated annealing algorithm for solving the bi-objective facility layout problem,” Expert Systems with Applications, vol. 38, no. 4, pp. 4460–4465, 2011.
View at: Publisher Site | Google Scholar
S. Kulturel-Konak, “A linear programming embedded probabilistic tabu search for the unequal-area facility layout problem with flexible bays,” European Journal of Operational Research, vol. 223, no. 3, pp. 614–625, 2012.
View at: Publisher Site | Google Scholar
X. H. Suo and Z. Q. Liu, “Modeling and solution algorithms for facility layout of manufacturing systems,” Computer Integrated Manufacturing, vol. 13, pp. 1941–1951, 2007.
View at: Google Scholar
M. J. Ye and G. G. Zhou, “The application of genetic algorithm in the Bi-criteria layout problem with aisles,” Journal of Systems Engineering—Theory & Practice, vol. 10, pp. 101–107, 2005.
View at: Google Scholar
G. Aiello, G. La Scalia, and M. Enea, “A multi objective genetic algorithm for the facility layout problem based upon slicing structure encoding,” Expert Systems with Applications, vol. 39, no. 12, pp. 10352–10358, 2012.
View at: Publisher Site | Google Scholar
J. M. Palomo-Romero, L. Salas-Morera, and L. García-Hernández, “An island model genetic algorithm for unequal area facility layout problems,” Expert Systems with Applications, vol. 68, pp. 151–162, 2017.
View at: Publisher Site | Google Scholar
L. García-Hernández, A. Arauzo-Azofra, L. Salas-Morera, H. Pierreval, and E. Corchado, “Facility layout design using a multi-objective interactive genetic algorithm to support the DM,” Expert Systems, vol. 32, no. 1, pp. 94–107, 2015.
View at: Publisher Site | Google Scholar
F. G. Paes, A. A. Pessoa, and T. Vidal, “A hybrid genetic algorithm with decomposition phases for the unequal area facility layout problem,” European Journal of Operational Research, vol. 256, no. 3, pp. 742–756, 2016.
View at: Publisher Site | Google Scholar
J. F. Gonçalves and M. G. C. Resende, “A biased random-key genetic algorithm for the unequal area facility layout problem,” European Journal of Operational Research, vol. 246, no. 1, pp. 86–107, 2015.
View at: Publisher Site | Google Scholar
R. D. Meller, V. Narayanan, and P. H. Vance, “Optimal facility layout design,” Operations Research Letters, vol. 23, no. 3–5, pp. 117–127, 1998.
View at: Publisher Site | Google Scholar
K. Y. Wong and Komarudin, “Solving facility layout problems using flexible bay structure representation and ant system algorithm,” Expert Systems with Applications, vol. 37, no. 7, pp. 5523–5527, 2010.
View at: Publisher Site | Google Scholar
B. Montreuil, N. Ouazzani, and E. Brotherton, “Coupling zone-based layout optimization, ant colony system and domain knowledge,” in proceedings of the 8th International Material Handling Research Colloquium, pp. 301–331, Graz, Cincinnati, OH, USA, June 2004.
View at: Google Scholar
D. J. van Camp, M. W. Carter, and A. Vannelli, “A nonlinear optimization approach for solving facility layout problems,” European Journal of Operational Research, vol. 57, no. 2, pp. 174–189, 1992.
View at: Publisher Site | Google Scholar
Q. Liu and R. D. Meller, “A sequence-pair representation and MIP-model-based heuristic for the facility layout problem with rectangular departments,” IIE Transactions, vol. 39, no. 4, pp. 377–394, 2007.
View at: Publisher Site | Google Scholar
Y. A. Bozer and R. D. Meller, “A reexamination of the distance-based facility layout problem,” IIE Transactions, vol. 29, no. 7, pp. 549–560, 1997.
View at: Publisher Site | Google Scholar
M. S. Bazaraa, “Computerized layout design: a branch and bound approach,” A I I E Transactions, vol. 7, no. 4, pp. 432–438, 1975.
View at: Publisher Site | Google Scholar
S. Kulturel-Konak and A. Konak, “Unequal area flexible bay facility layout using ant colony optimisation,” International Journal of Production Research, vol. 49, no. 7, pp. 1877–1902, 2011.
View at: Publisher Site | Google Scholar
Komarudin and K. Y. Wong, “Applying ant system for solving unequal area facility layout problems,” European Journal of Operational Research, vol. 202, no. 3, pp. 730–746, 2010.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Shiwang Hou 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.

PDF Download Citation

Download other formats

Order printed copies

Views

2649

Downloads

1408

Citations