Abstract

One of the main concepts in group technology (GT) is the cellular manufacturing system (CMS) with three main problems of cell formation (CF), cell layout (CL), and cell scheduling (CS). This paper studies the cell layout problem (CLP), aiming to find the optimal layout of machines within each cell (intracellular layout) and the optimal layout of cells in each workshop (intercellular layout). To adapt to reality, the dimensions of the cells and machines (inside each cell) were considered unequal, and also the cells and machines could rotate. We believe that a cellular layout that assumes unequal dimensions of the cells and machines can be used for batch production. This kind of production has a wide variety of low to medium demand. Furthermore, a cellular layout can be applied in CMSs and also in noncontinuous industries that have a job shop layout. Our main contribution is considering the possibility of rotating the cells and machines inside the cells. For this purpose, a mixed nonlinear programming model was developed to solve the CLP with the minimum cost of intracellular and intercellular material flows. The proposed nonlinear model was first converted into a linear model, and then a problem was generated and solved with GAMS software to validate the resulting linear model. This model finds the best layout of cells within the workshop and the best layout of machines inside each cell. Then, because of the NP-hardness of the CLP and the fact that even exact methods cannot solve large-scale examples in an acceptable computational time, an imperialist competitive algorithm (ICA) was designed and used to solve the problem. To evaluate the efficiency of the proposed algorithm, its numerical results in small dimensions were compared with the results of GAMS software. In large dimensions, 30 random problems were created, and the results of ICA were compared with the results of the particle swarm optimization (PSO) algorithm and genetic algorithm (GA). Finally, the parameters of the three meta-heuristic algorithms were set by the Taguchi method. Numerical results indicated that ICA was superior to both the PSO algorithm and GA. It could also achieve efficient solutions in a shorter computational time.

1. Introduction

One of the main purposes of manufacturing systems is to convert raw materials, capital, information, labor, and other resources into more value-added goods or services. Manufacturing systems consist of several specifically organized machines to construct different kinds of products. They have different types based on facility layout, flow shop system, job shop system, fixed position layout, and cellular manufacturing system (CMS). In a flow shop system, production facilities are set sequentially. Examples of this flow are metal production industries, refineries, and petrochemical industries. In a job shop production system, similar machines are located next to each other as in automotive and home appliance industries. In fixed position layout, the product is fixed, and each department performs the essential processes. Examples of such systems are aircraft, ship building, and construction industries. The last type of production system, i.e., CMS, is an application of group technology (GT) and comparing to other types of manufacturing systems like the flow shop system, has higher flexibility in the production of different types of products. Also comparing to job shop production system, which is used only for low production volume with high production and storage cost, CMS enjoys good efficiency with a higher production rate [1].

A cell may contain a group of similar machines in which the components must have the maximum movement within each cell and the minimum movement relating to other cells. This advantage causes in the reduction of material handling costs, process and setup times, and inventory, which ultimately leads to better quality, customer satisfaction, and on-time delivery of the product. For the effective implementation of CMS, there are three important tasks:(1)Grouping of parts with similar design features or processes in the family of parts, and grouping of machines into machine cells, which is called cell formation (CF).(2)Finding the optimal layout of machines within each cell (intracellular layout) and the optimal layout of cells in the workshop (intercellular layout), which is collectively called cell layout problem (CLP).(3)Determining the order of part families in each cell and the order of tasks in each part family, which is called cell scheduling (CS) [2].

The main objective of this paper is to develop a new mixed nonlinear mathematical model for CLP while the cells and machines have unequal dimensions and can rotate. The main questions of the paper are: “How to extend a new mixed nonlinear mathematical model for CLP?” and “How to solve it due to the Np-hardness of the model?” Accordingly, the main contributions of this research can be outlined as follows:(i)Unlike previous research findings implying that do not rotate, using the model developed in this research, cells and machines can rotate in different directions.(ii)In order to conform to the reality, unequal dimensions for cells and machines have been considered in this study.(iii)The main application of cellular layout is in the automotive industries like the Toyota manufacturing industry.(iv)In the present study, a new mixed nonlinear mathematical model has been proposed for CLP.

The following phases have been done in the present study. The related literature is reviewed in Section 2. Problem statement is explained in Section 3. A mathematical model of the problem is presented in Section 4 to minimize the intercellular and intracellular material flow. Since CLP is very complex and its computational time increases extremely as the size of the problem increases, exact algorithms for solving large-size problems lose their efficiencies; therefore, Section 5 compares the results of ICA with those of (PSO) and (GA). Finally, conclusions, study limitations, and future works presented in Section 6.

CLPs are generally divided into three types: intercellular layout, intracellular layout, and intra- and inter-cellular layout. In the first case, only the layout of cells in the workshop or intercellular layout is considered. The objective function, in this case, is to minimize the cost of material flow between the cells in the workshop. In the second case, the cells’ location is already known, and only the layout of machines within each cell or intracellular layout is considered. The objective function of this type of problem is to minimize the cost of material flow between the machines inside the cells. In the third case, the layout of cells in the workshop and the layout of machines inside the cells are considered, which is called intra- and inter-cellular layout. The objective function, in this case, is to minimize the cost of material flow between the cells and the cost of material flow between the machines within each cell. In the following, almost all the work done in the field of CLPs are reviewed. Some studies considered intra- and inter-cellular layout problems under dynamic situations and assumed equal dimensions for cells and machines. Others considered intracellular layout problems. The third group of studies mentioned intercellular layout problems with similar dimensions for cells and machines. Then a mathematical model and a solution method are proposed for the mentioned problems, and numerical results are presented.

Rosenblatt [3] studied factory layout problems under dynamic conditions and equal dimensions for machines. Vakharia and Wemmerlov [4] considered the problem of intracellular layout with equal dimensions and used a coefficient to assess the similarity of parts and the formation of a family of parts and cells. Two years later, Venugopal and Narendran [5] proposed the problem of cell formation (CF) to minimize the intercellular movements of parts, balance the capacity of cells, and solve the problem using GA. Meanwhile, Jajodia et al. [6] developed a new method for intracellular and intercellular layout problems with the equal dimensions for cells and machines. Arvindh and Irani [7] proposed the idea of designing a CMS in intracellular and intercellular layout with equal dimensions for cells and machines and used the production flow technique as a solution. Kim and Kim [8] assumed a model in which the shape of the machine was regular or irregular, its dimensions were fixed, the material handling system was based on an open field, and the objective function had a single criterion.

Vakharia and Cheng [9] developed the CF problem in an intercellular layout. The solution method was based on simulated annealing (SA) and tabu search (TS) algorithms. Consequently, Liang and Zolfaghari [10] proposed a problem for grouping machines with equal dimensions under uncertain and static conditions. Wang et al. [11] proposed a model for solving equipment layout problems in CMSs in which the demand rates changed over the product’s life cycle. Lozano et al. [12] developed machine grouping and intercellular layout problems and considered a mathematical programming model to minimize the cost of intercellular translocation. Solimanpur et al. [13] developed intercellular layout problems as a quadratic assignment. Anders and Lozano [14] proposed CF and intercellular layout problems by assuming a criterion for grouping machines in cells and used PSO algorithm to solve them. Tavakkoli-Moghaddam et al. [15] addressed intracellular and intercellular layout problems by assuming the possible demands and considered the fuzzy linear integer programming model. Hu et al. [16] inspected the intercellular layout problem and developed a model to minimize the cost of intercellular material flow. Chan et al. [17] proposed intercellular and intracellular layout problems by assigning machines within cells and considering equal dimensions for cells and machines. Tavakkoli-Moghaddam et al. [18] considered the problem of dynamic CF and intercellular layout under dynamic conditions. Mahdavi et al. [19] concentrated on the problem of intercellular and intracellular movements and CF in CMS. The following year, Rafiee et al. [20] proposed a CMS considering intracellular and intercellular layout problems in dynamic conditions. Ariafar et al. [21] applied intracellular and intercellular layout problems to CMS. Mahdavi et al. [22] developed a mathematical model for the simultaneous combination of CF and intracellular and intercellular layout. Shiyas and Pillai [1] considered the problem of intercellular layout and used a mathematical model to design manufacturing cells. Asl and Wang [23] developed an intracellular layout problem under uncertain and probabilistic conditions to minimize material flow costs. Mehdizadeh and Rahimi [24] proposed an integrated mathematical model for solving the dynamic CF problem, considering operator allocation and intracellular and intercellular layout problems. Ghosh et al. [25] considered the problem of intercellular layout and proposed a quadratic assignment planning model for such problems. Derakhshan Asl and Wang [26] developed the problem of intracellular layout in static and dynamic conditions, supposing unequal machine dimensions and the possibility of machine rotation. Rabbani et al. [27] hypothesized the problems of CF and intracellular and intercellular layout simultaneously in dynamic conditions and grouping of machines and parts. Feng et al. [28] investigated the CL and CF problems simultaneously and used GA and SA to solve them. Golmohammadi et al. [29] developed intracellular and intercellular layout problems and assumed machines, cells, parts, and batch sizes for transporting parts.

Mahmoodian et al. [30] developed a “smart” PSO algorithm for CLPs. Danilovic and Illic [31] worked on a hybrid algorithm to solve CLPs. Paramasamy et al. [32] developed a GA to form families of machine parts and cells and solved CLPs.

Rahimi et al. [33] considered CF, intercellular layout, and cell scheduling problems for CMSs. Zhao et al. [34] proposed a new layout methodology for the multifloor linear cellular manufacturing layout and a mathematical model. Neufeld et al. [35] pointed out several characteristics that specify the distinctiveness of cell scheduling. Golmohammadi et al. [36] developed a biobjective optimization model to integrate CF and intercellular or intracellular layouts. Rostami et al. [37] presented a multiobjective mathematical model for the simultaneous integration of virtual cellular manufacturing with the supply chain and new product development. Ayough et al. [38] integrated job assignment and job rotation scheduling problems and presented a novel multiperiod nonlinear mixed integer model.

Goli et al. [39] designed a fuzzy mixed-integer linear programming model for CF problems in CMS. Salimpour et al. [40] developed CF, intracellular and intercellular layout problems, and machines of unequal dimensions. Golmohammadi et al. [41] considered the problem of facility layout by integrating CF in continuous space and intercellular and intracellular layouts. Mondal et al. [42] hypothesized the problems of intercellular layout and CF in CMS and offered a model for minimizing cell load variation and reducing intercellular movements. Razmjoei et al. [43] hypothesized an intercellular layout in CMS. Hazarika [44] proposed an intercellular layout in CMS and surveyed the problem of machine CF. Mansour et al. [45] studied cellular manufacturing and intracellular and intercellular layout problems at the same time. Fakhrzad et al. [46] simultaneously addressed CF and intercellular layout problems in terms of scheduling, labor assignment, and limited financial resources. Finally, Forghani and Fatemi Ghomi [47] addressed CF and group layout problems and proposed a mixed integer programming model. A summary of the results of the above works are presented in Table 1.

As shown in Table 1, the current research was conducted on CLPs, including intracellular and intercellular layouts where the dimensions of cells and machines are considered unequal and rotation of cells and machines are allowed. However, in the previous research studies, these six factors were not considered simultaneously.

2.1. Research Gap

As explained above, a comprehensive review of the literature within the last three decades clarifies that no research has examined CLPs under unequal dimensions of cells and machines within each cell and the possibility of rotation of cells and machines. Hence, this study was designed to examine these types of problems.

3. Problem Statement

In CLP problem, to find the layout of cells in the workshop and that of machines within each cell, first the number of cells and machines and then the dimensions of the workshop, cells, and machines are determined. Next, the number of machines inside each cell is specified. Lastly, the material flow between the cells and between the machines is determined. For example, assume we have three cells and ten machines. The cells are positioned in a workshop (1.5 m in length and 1.5 m in width). The length of the first cell is 50 cm and its width is 45 cm, and the length of the second cell is 45 cm and its width is 40 cm. The length and width of the third cell are both 40 cm, and the dimensions of the machines are also known. The first, second, third, and fourth machines are inside the first cell, the fifth, sixth, and seventh machines are inside the second cell, and the eighth, ninth, and tenth machines are inside the third cell. There is a material flow between the first, second, and third cells called intercellular material flow. There is also a material flow between the machines inside each cell, called intracellular material flow. The layout of the first, second, and third cells in the workshop area is called intercellular layout, and the machine’s layout inside each cell is called intracellular layout. It is also possible to rotate all machines and cells. Figure 1 schematically displays an integrated intercellular and intracellular layout.

This study examines CLPs when the cells are located in the workshop and do not go beyond the workshop area, and the machines are located inside the cells and do not go outside. In addition, the cells and machines should not overlap. As you know, in the real world, cells and machines have unequal dimensions; so, in this study, in order to get closer to reality, the dimensions of cells and machines inside each cell have been considered to be unequal. Another considered aspect in this research is that, unlike previous studies where cells and machines did not rotate, here, they can rotate in different directions. The present study suggests a new mixed nonlinear mathematical model for CLPs. Because the proposed model is nonlinear, to solve small CLPs using GAMS software, the model was linearized, and the results were compared with the ICA results.

Since CLPs are NP-hard, we used ICA and other meta-heuristic algorithms like PSO and GA to solve large-scale problems. Finally, the results of ICA were compared with the results of the two other algorithms (PSO and GA).

3.1. Assumptions

In this research, a new mixed nonlinear mathematical programming model is proposed to solve intracellular and intercellular layout problems with unequal dimensions for cells and machines. The assumptions of the proposed mathematical model are as follows:(1)Manufacturing cells should be enclosed inside the workshop with exact dimensions.(2)Location of the cells inside the workshop is not predetermined.(3)Dimensions of the cells and machines are known and not equal (i.e., it is not an assignment problem).(4)The cells are not located over each other and do not overlap.(5)The machines used to produce various parts are enclosed inside the cells.(6)The machines are not located over each other and do not overlap.(7)Structural cells can rotate.(8)Machines can also rotate.(9)Material entry and exit points are always in the center of the cells and machines.

3.2. Notations
3.2.1. Indices
i, k: Cell indexj, p: Machine index
3.2.2. Parameters
: Length of the workshop: Width of the workshop: Number of cells: Length of the cell : Width of the cell : Length of the machine in the cell : Width of the machine in the cell : Number of the machines in the cell : Material flow between the and cells , : Material flow between the and machines in cell
3.2.3. Decision Variables
: Center coordinates of the cell along the x-axis : Center coordinates of the cell along the y-axis : Center coordinates of the machine in the cell along the -axis : Center coordinates of the machine in the cell along the -axis : Rectilinear distance between and cells, a function of and : and
3.3. Preliminaries

This model aims to find the best layout of cells in the workshop and the best layout of machines inside each cell.

3.4. Mixed Nonlinear Programming Model

Based on the above information, the mixed nonlinear programming model for intracellular and intercellular layout problems is as follows:

In the proposed model, the objective function consists of two parts. The first part shows minimizing the cost of material flow between cells, and the second part states minimizing the cost of material flow between machines located inside the cell. Equation (1) indicates minimizing the sum of material flow costs between cells, and equation (2) represents minimizing the sum of material flow costs between machines. In equations (1) and (2), the cost is considered one per unit distance. Equation (3) indicates intracellular and inter-cellular cost minimization. Inequalities (4) and (5) show that each cell is enclosed inside the workshop along the x- and y-axes. Inequalities (6) and (7) express that the machines of each cell are enclosed in the corresponding cell along the x- and y-axes. Inequality (8) indicates the nonoverlap of cells, while Inequality (9) expresses the nonoverlap of machines inside the cells. Inequalities (10) and (11) express that the center coordinates of the cell along the x- and y-axes are equal to or greater than 0. Inequalities (12) and (13) indicate that the center coordinates of the machine located in the cell, along the x- and y-axes, are equal to or greater than 0. Constraints (14) and (15) ensure that the decision variables, and are binary variables.

3.5. Linearization of the Mixed Nonlinear Programming Model

Since the proposed model is nonlinear, we should linearize the model. For this purpose, the following variables were changed in the objective function and constraints as shown in Table 2.

After applying the changes in the above variables, the objective function and model constraints were transformed into equations (16)–(21):

3.6. Solution Method

Ballakur and Steudel [50] described that CLPs are complex and increase in complexity when the cells and machines are unequally sized. For this reason, exact methods cannot solve the problems explained above during an acceptable computational time. Therefore, to test our proposed model, meta-heuristic algorithms should be used. To our knowledge, ICA has not been used until now in CLPs. Hence, this algorithm was implemented on CLPs, and its numerical results were compared with the numerical results of PSO and GA.

ICA creates an initial set of random solutions, known as chromosomes in GA, as particles in PSO, and as countries in ICA. Convergence rates and achieving near-optimal solutions are other advantages of ICA [51, 52].

3.6.1. Implementation of ICA to Solve the Problem

ICA was first introduced by Atashpaz-Gargari and Lucas in 2007 [51]. It is designed based on population. The original population is generated randomly, each member of which is known as a country. It steadily improves the initial answers (countries), and finally, creates the suitable answer to the optimization problem (optimal country). In the concept of optimization, some of the best members of the population that have the lowest cost are selected as imperialist countries, and the rest of the population is allocated as a colony to the imperialist country. Colonies in the original population are divided between the imperialist countries based on their power. The power of any imperialist is inversely related to the fitness (cost) of that country. Colonies and imperialist countries form empires [53, 54]. The steps to implementation of ICA are as follows.

(1) Creating Initial Solutions and Forming Empires. ICA is implemented in four phases, and each phase is executed to the maximum number of repetitions. In each iteration, the amount of material flow cost is calculated for all population members and selected as the number of empires from the countries with the lowest material flow cost, as imperialist countries. The selection of the best countries in our case was made using Roulette’s wheel method [55]. The remaining countries were then assigned as colonies to the imperialist country, and primary empires were created.

(2) Solution Representation. Previous studies have shown that displaying solutions as continuous is efficient and provides better solutions in terms of facility layout [23]. Therefore, in the present study, solutions were displayed in a continuous and linear manner. To implement the algorithm, in displaying shadow solutions for each decision variable, a number between 0 and 1 was mentioned. Four parts were used to represent the solution. The first part consists of 9 numbers; the first three of which are the shadow of the center coordinates of the cells along the x-axis (), the second three are the shadow of the center coordinates of the cells along the y-axis (), and the last three numbers show the shadow of rotation or nonrotation of the cells in relation to their original direction (). The second part contains the next 15 numbers; the first five numbers in this section are the shadow of the center coordinates of the machines in the first cell along the x-axis (), the second five numbers are the shadow of the center coordinates of the machines in the first cell along the y-axis (), and the last five numbers are the shadow of rotation or nonrotation of the machines in the first cell (). The third part includes the next 12 numbers; the first four numbers in this section are the shadow of the center coordinates of the machines in the second cell along the x-axis (), the second four numbers are the shadow of the center coordinates of the machines in the second cell along the y-axis (), and the last four numbers are the shadow of the orientation of the machines in the second cell (). The fourth part comprises the numbers 36 to 45; the first three numbers in this section are the shadow of the center coordinates of the machines in the third cell along the x-axis (), the second three numbers are the shadow of the center coordinates of the machines in the third cell along the y-axis (), and the last three numbers are the shadow of the orientation of the machines in the third cell (). Then, by inserting the obtained shadow coordinates in equations (22)–(24), we can determine the center coordinates of the machines along the x-axis (, the center coordinates of the machines along the y-axis () in the corresponding cell, and their orientation relative to their original direction ():

In this type of solution representation, a matrix with 1 row and 3 × M columns is applied. The above process is also used to display the solution for the cells. For example, in displaying the solution for three cells and 12 machines, a matrix with one row and 45 columns is utilized. Figure 2 displays the shadow of the center coordinates of machines and cells, while Figure 3 displays their actual center coordinates.

(3) Improving the Solution Space. At this phase, the position of the imperialist countries is changed using the absorption operator. By using the revolution operator, random changes are made in the material flow cost of the imperialist and colony countries to attain better and more diverse solutions. The movement of the colonies towards the imperialist country (following a policy of absorption) is done using the following equation:where is the secondary position of the colony, x is the primary position of the colony, β is the simulation coefficient, t is the value of the imperialist or the objective value, and 0 ≤ β ≤ 2. The absorption policy about colonies is shown in Figure 4.

(4) Replacement of the Colony and the Imperialist with Each Other in case of Meeting Relevant Conditions. Next, the material flow cost of the colonies is compared with the relevant imperialist country, and if the material flow cost of the colony is better than the material flow cost of the imperialist, the position of the colony and the imperialist country will be changed (the intragroup competition).

(5) Transfer of Colonies and Imperialists to Another Empire in case of Meeting Relevant Conditions. At this step, empires are first evaluated using equation (26) in order to obtain the index for each empire:where refers to the value of the imperialist objective function, denotes the average value of the objective function for the colonies, and ξ is the zeta coefficient, the value of which is set by parameter setting.

The weakest colony is removed from the weakest empire, and accidently, moved to a more powerful empire (intergroup or inter-empire competition). If the weakest empire has no colony, the imperialist, as a colony, is transferred to another empire. Then the value of the objective function of the empires is updated, and the same process is repeated.

(6) Reporting the Best Solution Found. In this section, the best solution found in the output of MATLAB software is shown.

(7) Checking the Stop Condition. Although there are several conditions for stopping the algorithm, in this paper, reaching a certain number of iterations or reaching the maximum iteration (MaxIt) in implementing the ICA is considered a stop condition. The pseudocode of the proposed ICA is presented in Figure 5.

4. Results

4.1. Computational Results

Here, numerical examples as well as the results of running them with GAMS (win32 24.1.3) and MATLAB R2015 are presented.

4.1.1. Solving a Small-Scale Problem to Assess the Proposed Model’s Validity

To evaluate the model’s validity, a small-scale problem was solved with GAMS software. The data of this example were taken from the Asl and Wong article [23]. This example comprised three cells and 12 machines. The cells were set in a 100 cm wide × 100 cm long workshop. Other parameters of the problem are presented in Tables 39.

The small-scale problem was solved by the above model in GAMS software. Since CLPs are very complex and exact solution methods cannot solve them in a reasonable computational time, we solved them using ICA, and the results were compared with those of GAMS software as shown in Table 10. The layout of cells in the workshop and the layout of machines inside each cell are shown in Figure 6.

4.2. Parameter Setting for Meta-Heuristic Algorithms Using the Taguchi Method

In meta-heuristic algorithms, several parameters must be set. Parameter setting is important because better solutions can be created using the right parameter setting. In this paper, parameter setting for ICA, PSO, and GA was performed using the Taguchi method [56]. In the Taguchi method, the factors and parameters of each algorithm are categorized into two groups of controllable and uncontrollable, and the effect of uncontrollable factors is minimized. Then, the parameters of the assumed algorithms are adjusted by considering the output of the Taguchi method, which is the main effects plot for the signal-to-noise ratios.

4.2.1. Parameter Setting for the ICA

ICA has two operators and eight parameters. The following symbols are used to set the parameters of this algorithm: MaxIt, npop, nEmp, alpha, beta, PRevolution, mu, and zeta, where MaxIt is the highest number of iterations, npop is the population size, nEmp is the number of empires, alpha is selection pressure, beta is similarity coefficient, PRevolution is the probability of revolution operator, mu is the revolution operator rate, and zeta is the average cost factor of colonies. The parameters of ICA are set according to Figures 7 and 8, and the results are presented in Table 11.

4.2.2. Parameter Setting for PSO

There are five parameters for PSO, and MaxIt, nPop, , , and W symbols are applied to regulate its parameters: MaxIt is the highest number of iterations, nPop is the number of particles, is the acceleration coefficient of the best personal answer, is the acceleration coefficient of the best overall answer, and W is the inertia weight coefficient. The results are presented in Table 12.

4.2.3. Parameter Setting for GA

GA has five parameters, and MaxIt, nPop, Pc, Pm, and Mu symbols are used to adjust these parameters: MaxIt is the highest number of iterations, nPop is the population size, Pc is the percentage of intersection operator for population, Pm is the percentage of mutation operator for population, and Mu is the percentage of mutation rate for selected chromosomes. The results are presented in Table 13.

4.3. Generating Random Problems

To evaluate the mixed nonlinear programming model and meta-heuristic algorithms, 30 random problems were designed in MATLAB R2015 software and solved by (ICA), (PSO), and (GA). Then, the results were statistically analyzed using analysis of variances. The mean values of the objective function related to three meta-heuristic algorithms are shown in Figure 9 and Table 14.

4.3.1. Results of Running Random Problems

After designing the random problems, we executed them and the center coordinates of cells in the workshop and the center coordinates of machines inside each cell for the 30th random problem after implementation of ICA, PSO, and GA were determined, and the results are shown in Tables 1517, respectively. Also, the layout of cells in the workshop and the layout of machines inside each cell for the three meta-heuristic algorithms are shown in Figures 1012, respectively. It is worth noting that the dimensions of problems 26 to 30 are presented in Tables 18 and 19.

To demonstrate a more accurate estimate of the criteria for meta-heuristic algorithms, the mean relative percentage difference () of layout cost and computational time for the 30 random problems created in Section 4.3 are shown in Table 20. Because each random problem was accomplished 5 times using meta-heuristic algorithms, the relative percentage difference (RPD) is calculated for the cost and computational time criteria. Then, the average values or () of the performances were computed and applied for comparison among the three proposed algorithms. The results are presented in Table 20. The RPD index was acquired using the following equation:where i is the random problem number, and j represents the algorithm used (j = 1 is ICA, j = 2 is PSO, and j = 3 is GA). Besides, is the best value achieved from the implementation of algorithm number j for random problem i, and is the best result obtained from implementing the algorithms on random problem i.

Figures 13 and 14 show the of implementation of three meta-heuristic algorithms in terms of layout cost and computational time criteria, respectively. Figure 13 shows the relative difference of the layout cost of ICA, PSO, and GA with the symbols COICA, COPSO, and COGA, respectively. Figure 14 denotes the RPD values of ICA, PSO, and GA with the symbols TOICA, TOPSO, and TOGA, respectively. The mean symbol represents the RPD of the layout cost, and NP symbol represents the number of random problems in Figures 13 and 14.

4.3.2. Normality Test of Three Meta-Heuristic Algorithms

The results of performing the normality test for the three proposed algorithms are displayed in Table 21. The results related to the significant value show that the layout cost of PSO and the computational time of PSO and GA are normal. Regarding the central limit theorem, since the number of samples is 30, the setup cost of ICA and GA and the computational time of ICA are also normal.

4.4. Analysis of Variance Test for Assuming Equality of Means of Layout Cost and Computational Time of Three Meta-Heuristic Algorithms (ICA, PSO, and GA)

To compare the mean layout costs of three proposed meta-heuristic algorithms using the analysis of variance (ANOVA) test, four conditions must be met: samples from each group or community should have a normal distribution, samples from each group or community should be random, the three communities should be independent of each other, and the variances of the three communities should be homogeneous. Since the first three conditions were met, we used Levene’s statistic to test the fourth condition, and its results are shown in Table 22. Because of the high performance of 95% confidence level, the tests were performed at this confidence level.

According to Table 22, because the significance level is greater than 0.05, there is no significant difference between the variances of the three communities; in other words, the variances of these communities are homogeneous. Therefore, to test the equality of the means of the three communities in terms of layout cost, the ANOVA test can be used at a 95% confidence level. The analysis of variance compares the mean of multiple populations. The null hypothesis in the analysis of variance states that values of the three algorithms are equal according to the layout cost criterion, and there is no significant difference between the mean layout costs of the three algorithms. Hypothesis one states that the mean value of the three algorithms differs in terms of the layout cost criterion.

ANOVA results for assuming equality of the mean layout cost of three meta-heuristic algorithms are provided in Table 23. As shown, because the value of the significance level is greater than 0.05, there is no significant difference between the mean layout cost of the three communities, and the assumption of equality of means at 95% confidence level is accepted.

To suppose the equality of the mean computational time of the three communities, the above three conditions are met; to test the fourth condition, we used Levene’s statistic (Table 24).

Because the significance level is less than 0.05, the null hypothesis is rejected; so, the variances of three communities or three meta-heuristic algorithms are not equal. The results of the corresponding ANOVA are given in Table 25. As shown, since the significance level is less than 0.05, the null hypothesis is rejected (i.e., a significant difference exists between the mean computational times of the three communities). Therefore, the Tukey test should be executed to check the possible inter-group differences. The results of the Tukey test in Table 26 confirm that the mean computational time of all the three communities is different.

Since the variances are not homogeneous, the validity of the ANOVA test may be questioned. Therefore, Tamhane’s T2 test was used and the results are displayed in Table 27. Fortunately, Tamhane’s T2 test results manifested the previous results; i.e. the mean computational time of the three meta-heuristic algorithms is different significantly.

A comparison of the mean computational time of the three proposed meta-heuristic algorithms is shown in Figure 15. Numbers 1, 2, and 3 on the x-axis represent ICA, PSO, and GA, respectively. The mean computational time of ICA is much less than that of the other two algorithms. Therefore, based on the numerical results achieved from the ANOVA, it is obvious that ICA has much better results in terms of computational time comparing to the other two algorithms.

4.5. Discussion

In this paper, we first solved a small-scale problem with GAMS software to evaluate the model validity. Then, we compared the results of GAMS software with the results acquired from ICA. Numerical results presented the validity of the model as well as the efficiency of the proposed algorithm. Next, parameter settings for meta-heuristic algorithms were executed using the Taguchi method. In the next step, 30 random problems were created using MATLAB R2015 and performed by three meta-heuristic algorithms. After designing and implementing random problems with ICA, PSO, and GA, the center coordinates of the cells and machines were determined. To demonstrate a more accurate approximation of the criterion for meta-heuristic algorithms, we determined the values of layout cost and computational time. Also, the results of the three meta-heuristic algorithms in terms of layout cost and computational time criterion were obtained. Then the ANOVA was executed for assuming the equality of means of layout cost and computational time of three meta-heuristic algorithms aiming to compare their mean layout costs. The results denoted the equality of means of layout cost of the three proposed meta-heuristic algorithms. However, the computational time results showed the inequality of the means of computational time of all the three meta-heuristic algorithms. Therefore, the Tukey test was used to investigate the possible inter-group differences. According to the results, the mean computational time of all algorithms was different. Finally, the comparison of the mean computational time of the mentioned algorithms indicated that the mean computational time of ICA was much less than the other two algorithms; in other words, it demonstrated much better results in terms of computational time.

5. Managerial Insights and Practical Implications

Based on the research findings, a cellular layout that assumes unequal dimensions of cells and machines can be used for batch production of products. Such products have a wide variety and low to medium demand. In addition, cell layout can be used in CMSs for moving from traditional job shop layout to cellular layout. Likewise, cellular layout is used in noncontinuous industries that have a job shop layout. In continuous layout, for example, one may consider a factory that produces four parts by machines 1, 2, 3, and 4 using the same fixed procedure. In job shop layout, there is no need for the production procedure of all parts to be the same.

Other managerial insights are presented as follows:(i)The flow of materials as an influencing factor in the links between machines in different cells must be determined suitably because the incorrect input of this parameter has a tremendous negative effect on the formation of cells and increases the layout cost.(ii)The use of a meta-heuristic algorithm compared to exact solution methods is much better in terms of the speed of determining the layout, and managers are suggested to measure the efficiency of these algorithms in high-risk layouts while using multiple algorithms.(iii)Cell layout technology can reduce the layout cost for many production units with continuous material flow and increase the speed of material flow.(iv)Because CLPs that assume unequal dimensions for cells and machines are very difficult, entry and exit points of cells should be in their center to solve such problems.(v)Consequently and in order to solve CLPs, entry and exit points of machines should be considered in their center.(vi)The flow of materials between cells and machines can be nondeterministic; in other words, it can be probabilistic and dynamic.

These points can be determined in future works.

6. Conclusions

In this paper, a new mixed nonlinear programming model was presented for CLPs under the condition of unequal dimensions of cells and machines. Furthermore, it was assumed the cells’ location in the workshop and the machines’ location inside the cells are not predetermined, and the cells and machines can rotate. Based on the research findings, a cellular layout that assumes unequal dimensions of cells and machines can be used for products with batch production, and cell layout can be applied to CMSs aiming to move from traditional job shop layout to cellular layout. In the same way, cellular layout can be used in noncontinuous industries that have a job shop layout.

In brief, the following conclusions were made.

Due to nonlinearity of the developed model and failure to reach an optimal solution, the model was linearized. Then in order to assess its validity, a small-scale problem was solved by GAMS software and the results approved the model’s validity.(i)Since the problems considered in this research were continuous, (ICA), an efficient and continuous algorithm, was used to solve the small-scale problem. Comparison of the ICA and GAMS software showed that the presented model was valid.(ii)To evaluate the mixed nonlinear programming model and meta-heuristic algorithms, 30 random problems were designed and solved by three algorithms (ICA, PSO, and GA). Then statistical analysis was done on the results.(iii)Similarly, the mean values of the objective function related to three meta-heuristic algorithms were determined.(iv)The results related to the center coordinates of the cells in the workshop and the center coordinates of the machines within each cell were presented in Tables and Figures after implementation of ICA, PSO, and GA for the 30th random problem.(v)The values of the layout cost and computational time for all the three mentioned algorithms were determined.(vi)Normality and variance homogeneity tests were performed before doing the ANOVA test on the results obtained from the three meta-heuristic algorithms.(vii)Then ANOVA test results indicated no significant difference between the mean costs of the three meta-heuristic algorithms; however, there was a significant difference between the mean computational times among the three algorithms. Therefore, the Tukey test was used to test which two groups differed in terms of computational time. The Tukey test results also revealed that the mean computational times of all the three algorithms were different.(viii)Finally, it was found that ICA had reached the optimal solution in a very short computational time compared to the other two meta-heuristic algorithms. Thus, the superiority of ICA over PSO and GA is evident.

6.1. Study Limitations

There are always limitations in any research, and the current study is no exception. The limitations of the current research are as follows:(1)Considering rectangular shapes for the machines to solve cell layout problems.(2)Considering the entry and exit points of machines in their center.(3)Considering the entry and exit points of cells in their center.(4)Not considering stochastic condition for cells and machines.

6.2. Future Works

For future investigations, other researchers can further develop this research in the following directions:(1)The problems of cellular layout can be investigated in dynamic and probabilistic modes.(2)ICA provides encouraging results for CLPs when the material flow is considered at the center of the cells and machines. In the future, ICA can be applied for CLPs while considering entry and exit points for cells and machines.(3)ICA is an optimization method that can reach a globally optimal solution in less computational time. Some parts of this algorithm can be modified to reduce the total cost value.(4)ICA can cope with different types of optimization problems. In the future, new constraints and multiobjective cases can be considered.(5)The cells and machines in this research had a rectangular shape. In the future, irregular shapes can be considered for them.

Data Availability

The data used to support the findings of this study are available in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.