Open access peer-reviewed chapter

A Methodology to Design and Balance Multiple Cell Manufacturing Systems

Written By

Luis Valdivia and Pedro Palominos

Submitted: 11 May 2019 Reviewed: 02 September 2019 Published: 27 November 2019

DOI: 10.5772/intechopen.89463

From the Edited Volume

Mass Production Processes

Edited by Anil Akdogan and Ali Serdar Vanli

Chapter metrics overview

801 Chapter Downloads

View Full Metrics


Manufacturing cell formation and its balance in just-in-time (JIT) type production environments have usually been studied separately in the literature. This practice is unrealistic since both problems interact and affect each other when the cells are operating. This chapter proposes a methodology to design multiple manufacturing cells and simultaneously balance their workload. The cells considered are U-shaped and process mixed models of product families. A nonlinear integer programming mathematical model is proposed, which integrates cell formation and their balancing, considering various production factors. For illustration, the method is applied to the redesign of a rack manufacturing process.


  • manufacturing cells
  • assembly line balancing
  • N U-lines
  • mathematical programming
  • mixed model production

1. Introduction

Group technology (GT) can be defined as a manufacturing philosophy identifying similar parts and grouping them together to take advantage of their similarities in manufacturing and design [1, 2]. Cellular manufacturing (CM) is an application of GT and has emerged as a promising alternative manufacturing system [3]. When a productive system is changed to make it cellular, it implies solving the manufacturing cell formation problem (MCFP), which means identifying groups of machines and associating them with product families so that the intercellular traffic that the products can have within the productive system is minimized. This problem has been approached historically by analyzing the machine-product incidence matrix (A), where each row represents a machine and each column represents a product, with each element aij equal to one if machine i processes product j, and equal to zero otherwise. When this matrix is partitioned arbitrarily, it is usual to have products that remain outside the diagonal blocks (cells), which are called exceptional elements, since they carry out intercellular movements. Papaioannou and Wilson [3] reviewed the approaches between 1997 and 2008 to solve the above problem, proposing taxonomy based on the solution methodologies. It must be kept in mind that the latter approaches have started taking into account production factors other than the incidence matrix, like processing time, demanded production volumes, and operation sequences. Most recent works [4, 5, 6, 7] are oriented mainly to the heuristic and metaheuristic approach to solve the problem, without considering other factors that appear when there are exceptional products that carry out intercellular movements, implying that these different production cells are linked with one another due to the precedence restrictions of these exceptional products, without considering their cost and work under way. This means that in most real cases, it is not possible to analyze the manufacturing cells independently at the time of attempting to balance the load of their work stations, and furthermore, since those cells would be related with families of products, it is therefore possible to introduce the balancing concept of multiple manufacturing cells for mixed models. This happens because cellular manufacturing is commonly used in JIT-type productive systems, in which setup times can be reduced in such a way that each cell works by operating over a family of mixed product models. These cells can also be configured in U shape to use the advantages generated by this configuration. An approach to the above situation is the work of Kumar et al. [8], who propose implementing heuristic cell formation, having the capability to handle production data, operation sequence, production volume, and inter-cell cost simultaneously, taking up some of the previously described elements.

On the other hand, the problem of balancing N U-shaped lines has been studied mainly by Sparling [9] and Miltenburg [10], both of whom considered that all the cells operate with a common cycle time (C), but they assume that each cell is independent of the rest and furthermore process a single product. The problem of balancing U-shaped lines for mixed product models (denoted by MiMULBP) was proposed for the first time by Sparling and Miltenburg [11], who used the classical combined precedence graph proposed by Thomopoulos [12] and considered as cycle time (C) the quotient of the time period (T) and the total product demand (D). These authors focused only on the problem of balancing a single cell, making in the appendix the observation that it is possible to consider systems with multiple manufacturing cells, although once again they considered that those cells are independent of one another. More recently, in [13, 14], new heuristics are introduced to solve the problem, and Turkay [15] proposes models of integer linear programming (MILP) considering restrictions that express the precedence of the tasks.

A very recent work [16] proposes a novel configuration of assembly lines, namely parallel adjacent U-shaped assembly lines (PAUL), but none of the revised works integrate the balancing of U-lines with the design of the manufacturing cell. In the present chapter, what is being sought is to integrate these problems, proposing a methodology that delivers cells more applicable to reality, thereby introducing the problem of the formation and balancing of N in U-shaped cells for mixed models (denoting it by N-MiMUCFBP). The rest of this chapter is organized as follows: Section 2 introduces a methodology based on a mathematical model for the N-MiMUCFBP; Section 3 illustrates the proposed methodology using a real case, showing its results; and Section 4 gives the conclusions of the study.


2. Methodology

The proposed methodology is based on formulating a new model for the problem of balanced formation of production cells. From the viewpoint of formation of the cells, the model must consider the aspects associated with their design, such as processing and preparation of machines, inefficiencies in the handling of materials or inventories of products being processed, and cell imbalance [17], which includes processing times, sequences, and production volumes, directly related with the mixed model assembly line balancing problem (MiMALBP). Because of this, and in agreement with the heuristic proposed by [12] for the MiMALBP, first it is necessary to group the precedence graphs of each final product in a single combined precedence diagram. This is possible, because in this type of mixed model assembly lines, the products that are processed in each cell have only small differences in processing times or in the elimination or addition of activities but always keeping the consistency among the precedence of these activities for the different product models. Therefore, the idea is for each manufacturing cell to process a family of products. In turn, using this combined precedence diagram will allow reducing the number of variables at the time of tackling the model for the formation of manufacturing cells. The methodology proposed for the N-MiMUCFBP problem considers five consecutive stages that are presented in Figure 1.

Figure 1.

Stages in the methodology proposed for the N-MiMUCFBP problem.

2.1 Calculation and assignment of the required machines

First it is necessary to determine the number of machines, qm, required per type of machine m. This number is obtained from (Eq. (1)).

q m j = 1 n d j t mj CAP m E1

where qm, number of machines of type m; j product number; dj production volume demanded of each product j within the planning horizon (in units); tmj unit processing time for product j that is processed in a type m machine (in hours per unit produced); CAPm capacity of each type of machine m within the planning horizon (in available machine hours).

It should be noted that when a machine of a certain type can work simultaneously on a product together with another machine of the same type and qm > 1, the machines will operate as a single “virtual machine,” i.e., that qm machines of type m (denoted by i = m) will work simultaneously in it and that the unit processing time of each of these machines will be tij = tmj/qm, provided that the products processed on that virtual machine are similar to each other. On the other hand, for that type of machine in which qm > 1 and which furthermore can operate simultaneously on a product together with another machine of its same type, it is necessary to assign first which products will be processed on each machine of type m. For that purpose, let us use the following notation:

  • CAP i Useful = capacity used by machine i.

  • CAP m Level = j = 1 n d j t mj q m = leveled reference capacity used by qm machines of type m.

  • wmj = workload to be assigned of model of product j on the type of machine m.

  • d j i = fraction assigned to machine i of the demand for product j.

  • Bm = arrangement that contains the models of products similar to each other processed on type m machine, arranged in decreasing order according to wmj.

  • tij = Unit processing time on machine i due to product j.

Then the following procedure, shown in Figure 2, is proposed as a formal assignment rule.

Figure 2.

Procedure for assigning products to machines.

In particular, when qm > 1 and the type m machines cannot work simultaneously, the above procedure aims to assign similar products to each of the qm machines, so that they have a workload as close as possible to the leveled load for that type, and it also attempts to have the demand for each product as little fractionated as possible, so that each product is assigned to a single machine when it has sufficient capacity.

This procedure is followed in order to not incorporate directly in the mathematical model the alternative processes and routes, because in this way its complexity and number of variables are reduced. Special care must be taken when assigning similar products (with respect to their precedence relations) to the different machines, to minimize probable intercellular motions.

2.2 Preparation of the extended combined precedence diagram

With the machines required to satisfy the capacity restrictions, and the assignment of each product to them, an extended combined precedence diagram called GUG’ must be created, and the weighted average processing times for each machine must be calculated. The process for preparing this diagram will be described now by combining the precedence diagrams of each model in a single precedence diagram where the nodes represent the operations and the arcs represent the precedence restrictions between the operations. A formal description of the combination of n product models in a combined precedence diagram was made by Macaskill [18]. This procedure, adapted to our problem, is summarized as follows: Represent the precedence diagram of product model j by means of the graph Gj = (Vj, Ej, tj), where the set of nodes Vj represents the set of tasks of product model j, the set of arcs Ej represents the precedence relations (a, b) between tasks a, bVj, and the weighting vector tj contains the processing times tij of task iVj.

As an example, in Figure 3 the precedence diagrams for six models are represented, remarking the virtual machines in which more than one machine operate simultaneously on the products.

Figure 3.

Precedence diagrams for six product models.

Furthermore, by specifying the demanded volumes of each product model within the planning horizon (dj), it is possible to determine the demand fractions df j of each model j with respect to the total demand D of the product mix, where 0 ≤ df j ≤ 1 is fulfilled, and they are calculated by the following equation (Eq. (2)):

df j = d j D = d j j = 1 n d j E2

Therefore, the combined precedence diagram can be represented by the graph G = V E t ¯ , which is derived from the following definitions (Eqs. (3)(5)):

V = j = 1 n V j E3
t ¯ i = j = 1 n df j d j i t ij i V E4
E = j = 1 n E j redundant arcs E5

As a prerequisite for the generation of the combined set of nodes V in Eq. (3), the tasks that are common to different models, even though they have different processing times, receive a consistent number of nodes for all the models. This prevents assigning these tasks to different stations, which otherwise would need multiple investments in the resources required at each station in which a duplicate task has been assigned. Tasks that are not required by a product model receive a processing time (weight of the node) equal to zero, so the average processing times t ¯ i can be calculated simply by weighting every specific task fractionated time according to model d j i t ij with its corresponding demand portion df j of the model in Eqs. (4) and (5), which determines the combined precedence restrictions by joining the arc sets of each model. This can lead to redundant arcs (a, b), which represent the transitive precedence relations. An arc is redundant and can therefore be deleted without loss of information, if there is another way from node a to node b by means of more than one arc.

The combined precedence diagram for the example is shown in Figure 4. Note that the redundant arcs are denoted by dotted lines.

Figure 4.

Combined precedence diagram for the example.

A particular action should be considered if there is no consistency among the precedence of the activities, i.e., if there are conflictive precedence relations between different models that lead to a cyclic (that repeats itself over and over) combined precedence graph. To allow a single sequence of task operations, those loops must be deleted by means of one of the following actions proposed by Ahmadi and Wurgaft [19]:

  • The models must be separated into subsets so that two or more acyclic precedence graphs can be formed. In practice, this leads to machine preparation operations that must be performed every time the production changes from one subset of models to another.

  • With the purpose of assigning the tasks to a single station, the loops in the precedence graphs can be deleted duplicating these nodes. To minimize the number of duplicate nodes and, in this way, reduce the danger of assigning equal tasks to different stations, an optimization problem must be solved [19].

To model what is related to the problem of balancing the N U-shaped lines, we will consider the concepts developed by Urban [20] to formulate the problem mathematically. For this we set up an auxiliary graph, connecting it with the original combined precedence graph. This is illustrated in Figure 5, denoting with dotted lines the auxiliary combined precedence graphs. If we start in the middle of this extended graph, it is possible to perform assignments to stations forward through the original graph, backward through the auxiliary graph, or simultaneously in both directions, and in this way, it is possible to create stations that have machines at the beginning and at the end of the “U” line. Special care must be taken when joining the auxiliary precedence diagram with the original. For example, final task 15 is joined only with initial task 5, because task 15 is finished only from task 5 for model 6 (see Figure 3).

Figure 5.

Diagram of extended combined precedence.

2.3 Calculation of the parameters needed for the mathematical model

From the input information and the extended combined precedence diagram produced in the previous point, we must calculate the total costs for intercellular transport between machines i and i’ (c ii), which are determined by means of Eq. (6):

c ii ' = j = 1 n co j e ii ' j d j d j i + d j i ' 2 a ij a i ' j E6


  1. cii’ = cost of intercellular transport between machines i and i’ (within the planning horizon).

  2. coj = intercellular transport cost of product j.

  3. e ii ' = 1 if machines i and i ' are directly or indirectly related to the GOUGA graph 0 otherwise a ij = 1 if machine i processes product j 0 otherwise

  4. a i ´ j = 1 if machine i ´ processes product j 0 otherwise

  5. dj = production volume demanded by product model j.

  6. d j i = fraction assigned to machine i of the demand for product j.

  7. d j i ´ = fraction assigned to machine i´ of the demand for product j.

Other parameters to be used in the model, some of which must be calculated, are the following:

  1. cs = unit cost per work station (within the planning horizon).

  2. V = set of machines of the combined precedence diagram, V = j = 1 n V j .

  3. E = set of precedence relations between the machines that belong to V, E = 1 e E . For example, e = (a, b) is the ordered pair that indicates that machine a precedes machine b immediately,

E = j = 1 n E j redundant arcs .

  1. T = period of time available for planning.

  2. tij = processing time on machine i of product model j.

  3. dj = production volume demanded by product model j.

  4. D = total demand for the product models, D = j = 1 n d j .

  5. C ¯ = common average cycle time, C ¯ = T / D .

  6. dfj = fraction of total demand for product model j, df j = d j / D .

  7. θ = maximum number of production cells to which a multicellular station can belong 2 θ K max .

  8. Kmax = maximum number of production cells (Kmax ≤ m). It can be specified based on diagram G, or the maximum bound (K´) proposed by Al Kattan [21] can be used as reference:

K ' = W max i = 1 , , m w i K max = K '


  1. wi = marginal workload of machine i, w i = j = 1 n df j d j i t ij

  2. W = total workload, W = i = 1 M w i .

  3. Kmin = minimum number of production cells. It can be specified based on diagram G, or a modification of the bound proposed by Al Kattan [21] can be considered:

K min = max i = 1 , , m w i W / m .

  1. t ¯ i = weighted average processing time of machine i,

t ¯ i = j = 1 n df j t ij .

  1. GO = original combined precedence diagram. G = (V, E, t ¯ i ).

  2. GA’ = auxiliary combined precedence diagram. G´ = (V, E′, t ¯ i ).

  3. Mmin = minimum number of machines that a cell must contain, M min = m / K max .

  4. Mmax = maximum number of machines that a cell must contain, M max = m / K min .

  5. Smin = minimum number of work stations required, S min = i = 1 m t i / C ¯ .

  6. Smax = maximum number of work stations required (Smaxm).

2.4 Mathematical model

2.4.1 Model assumptions

The assumptions of the proposed model are presented below:

  • Multiple similar models of a product are produced if possible in a single manufacturing cell.

  • The processing times of the tasks are known and constant.

  • Each task of the combined precedence diagram is performed for at least one product model.

  • The average time of each task is not greater than the average cycle time C ¯ .

  • Each machine is assigned to a single work station for each product.

  • The operators can work in or out of the “U” cell.

  • In a manufacturing cell, the precedence restrictions are consistent among the different product models produced in it, i.e., if task a precedes task b in some model, then there is no other model in the cell in which task b precedes task a.

  • The precedence graphs of the product models are not fractionated, i.e., all the tasks to produce a product model are joined together directly or indirectly.

  • The mix of models, i.e., the demands for models within the planning horizon, is known with certainty (static problem).

  • There are no buffers between the work stations, so it is not possible for the work stations to operate at different production rates.

  • The workers are capable of performing any task in the manufacturing cell, i.e., they are capable of operating any machine.

  • The setup times of the machines are not significant.

  • The displacement times of the workers in a manufacturing cell are not significant, but not so out of them.

  • The machines related by some intercellular product movement must belong, if possible, to the same work station.

2.4.2 Notation

The following points present the notation and the proposed model:

  1. Indexes

    h = work station (h = 1,…, Smax).

    i, i’ = machine (i = 1,…, m).

    j = product model (j = 1,…, n).

    k = production cells (k = 1,…, Kmax).

  2. Decision variables

    xik = 1 if machine i is assigned to cell k 0 otherwise .

    yk = 1 if production cell k is used , i . e . , if it is assigned machines 0 otherwise .

    uihG = 1 if machine i of graph G is assigned to work station h 0 otherwise .

    rh = 1 if work station h is used , i . e . , it is assigned machines 0 otherwise .

    fkh = 1 if cell k is used by work station h 0 otherwise .

    gh = 1 if work station h is multicellular 0 otherwise .

2.4.3 Model Objective function

The objective of this formulation is to minimize the total cost of intercellular transport between machines, which will appear every time there are finished products between machines i and i’ and they belong to different production cells. Decision variables xik define the set of groups of machines, while the product families will be defined after the solution of this model, supported by the information of the extended combined precedence graph (GOUGA):

min Z = i = 1 m 1 i ' = i + 1 m k = 1 K max c ii ' 1 x ik x i ' k + cs h = S min + 1 S max r h E7

Restrictions for the formation of manufacturing cells:

k = 1 K max x ik = 1 ; i = 1 , , m E8
i = 1 m x ik M max y k ; k = 1 , , K max E9
i = 1 m x ik M min y k ; k = 1 , , K max E10

Restrictions for the balance of U-shaped cells:

h = 1 S max u ihGO + u ihGA = 1 i = 1 , , m E11
i = 1 m t ¯ i u ihGO + u ihGA C ¯ r h h = 1 , , S max E12
h = 1 S max S max h + 1 u ahGO u bhGO 0 a b GO E13
h = 1 S max S max h + 1 u bhGA u ahGA 0 a b GO E14

Linking restrictions between the formation and the balance of the cells:

i = 1 m u ihGO + u ihGA x ik m f kh k = 1 , , K max ; h = 1 , , S max E15
h = 1 S max g h f kh 2 k = 1 , , K max E16
k = 1 K max f kh θ 1 g h + 1 h = 1 , , S max E17

Restrictions for defining binary variables:

x ik , y k , u ihG , r h , f kh , g h 0 1 h , i , j , k E18

The set of restrictions (8) restricts each machine to a single cell. The set of restrictions (9) restricts each created cell to a maximum of M max machines, while the set of restrictions (10) restricts them to a minimum of M min machines; the value of yk is equal to one for the first K min restrictions, since it is known that these cells are required. As to the balance of the lines, the objective is to minimize the cost per required work station in addition to the theoretical minimum, avoiding the need to have rh variables for stations 1 through S min . The set of restrictions (11) ensures that each machine is assigned to only one station, either in the original precedence graph or in the auxiliary one [18]. The set of restrictions (12) ensures that for every station, the sum of the weighted average processing times of their assigned machines does not exceed the average cycle time; the values of rh are equal to one for the first S min restrictions. The set of restrictions (13) and (14) force the precedence restrictions between the machines; these relations are reversed for the auxiliary graph.

The set of restrictions (15) makes each variable fkh be equal to one when cell k is used by station h. The set of restrictions (16) allows a maximum of two multicellular work stations for each manufacturing cell, so that there are not many interferences between stations [9]. The set of restrictions (17) limits to θ the number of cells to which a multicellular work station can belong, and it also makes every variable gh equal to one when h is a multicellular station. The set of restrictions (18) defines the decision variables xik, yjk, uih, rh, fkh, and gh as binary.

2.5 Assigning the product models to the obtained production cells

With the groups of machines obtained, we must now assign the product models to each resultant cell. For this, let:

Pkj = Number of machines of cell k that process product j,

P kj = i = 1 m a ij x ik

Ω j = Set of cells that have the maximum number of machines that process j,

Ω j = k P kj = max k = 1 , , K P kj

W j k = Total workload (in hours) of cell k due to product j,

W j k = i = 1 m d j t ij x ik

Ψ j = Set of cells having maximum total workload due to product model j and belonging to Ω j ,

Ψ j = k W j k = max k Ω j W j k

Then, to assign to which production cell each product j belongs, the following formal procedure is defined, where three cases can occur:

  1. Case 1: If Ω j = 1 j J k P kj = max k = 1 , , K P kj ; assign each product to the cell where it will be processed by more machines.

  2. Case 2: If Ω j > 1 Ψ j = 1 j J k W j k = max k = 1 , , K W j k ; if there is a tie it must be assigned to the cell in which the product spends most processing time.

  3. Case 3: If Ω j > 1 Ψ j > 1 , assign j arbitrarily to a cell that belongs to Ψ j ; if a tie occurs, assign the product to the cell with a smaller number of machines or randomly.


3. Application of the methodology: illustrative case

The company in which the proposed methodology will be applied is of the metalworking type, making storage products (racks) for the retail industry. The final products are assembled at the customer’s facilities, so the work orders are divided considering the final product’s components, which are generally pillars, beams, struts, slotted angles (ANRA), and accessories. Each of them can have various modifications in size, processing times, and complexity in its operations (precedence restrictions), so it is possible to identify them previously as families of products, therefore complying with the observations proposed by Burbidge [22], who says that a system can be naturally susceptible to be transformed into one of the cellular manufacturing type.

The company has a factory in the commune of Quilicura, in the Metropolitan Region of Chile, and has 31 machines that can be classified into 17 types, processing 67 different product models. The place where the methodology will be applied has a job shop-type configuration.

This study will consider a time planning horizon of 3 months, in which the plant operates 16 h/d from Monday to Saturday. According to the company’s policies, 1 day per month is devoted to preventive maintenance operations of each machine, so it will be considered that each machine has a capacity of 1248 h, within a planning horizon of T = 1296 h. Tables 1 and 2 present general information with respect to the types of machines and the different product models, respectively.

3.1 Step 1: calculating and assigning the required machines

Applying Eq. (1) to determine the number of machines of type (qm) needed to satisfy the capacity restrictions, the values presented in Table 3 are obtained.

m (type of machine) Identification characteristic Present number Simultaneous machine
1 Sheet metal cutter 3 No
2 Strippit 1 No
3 Punching machine 1 No
4 25 tons press 2 No
5 45 tons press 2 No
6 55 tons press 2 No
7 90 tons press 1 No
8 Pneumatic press 160 tons 2 No
9 Sheet metal bender 2000 m 2 No
10 Sheet metal bender 3000 m 2 No
11 Sheet metal bender 4000 m 2 No
12 Sheet metal bender 160 ton 2 Yes
13 Welder accessories 2 No
14 Beam welder 2 No
15 Connector welder 2 No
16 Pillar welder 2 No
17 Forming machine ANRA 1 No

Table 1.

Information on the types of machines.

J (products) Identification characteristic
1–4 Strut
5–29 Beam
30–51 Accessory
52–59 Pillar
60–67 Slotted angle (ANRA)

Table 2.

Information on the types of products.

m 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
qm 3 1 2 2 2 2 1 2 2 2 2 2 2 2 1 2 1

Table 3.

Values of qm obtained for the industrial problem.

Comparing these results with the actual values of the number of machines of a given type present in the plant, we get that they are equal for almost all the values of m, except in the following two cases:

  1. m = 3, which means that one more machine must be bought to fulfill the required capacity of this type of machine.

  2. m = 15, where there is currently one extra machine.

Case (a) explains in some way why machine 3 punching machine has become a bottleneck in the plant when “accessory”-type products are made. Maintain, in this study use will be made of the current excess capacity produced by the two type 15 machines (connector welder) reflected in case (b), assigning products to both machines.

Applying the proposed assignment procedure of Figure 2, we get the assignment of machines (i) to the (m) machine types and of products (j), which are presented in Table 4, where the fact that machine 12 will function as a “virtual machine” in which two machines will operate simultaneously is pointed out.

i m j i m j
1 1 1–4, 10–19, and 60–67 17 17 60–67
2 2 60–67 18 1 5–9 and 20–34
3 3 30–38 19 1 35–51
4 4 15–23 20 3 39–47
5 5 30–34 21 4 22–29
6 6 5–9 and 20–23 22 5 43–51
7 7 60–67 23 6 22–29
8 8 52–55 24 8 56–59
9 9 30–34 25 9 43–51
10 10 30–38 26 10 39–51
11 11 1–4 and 10–23 27 11 5–9 and 22–29
12 12 52–59 28 13 39–51
13 13 30–38 29 14 5–9 and 22–29
14 14 10–23 30 15 43 and 44
15 15 45–47 31 16 56–59
16 16 52–55

Table 4.

Results of the assignment of machines and products.

3.2 Step 2: preparation of the extended combined precedence diagram

With the determination of the number of types of machines, and the assignment of products to machines obtained in the previous stage, the precedence restrictions deliver the extended combined precedence diagram shown in Figure 6, where the virtual machine 12 differs from the rest because in it the two benders with the largest capacity operate simultaneously in the production of pillar-type products. These products have a larger size, so it is advisable to use both benders one next to the other to mechanize the product and, in this way, reduce the processing time of this operation.

Figure 6.

Extended combined precedence diagram for the industrial problem.

3.3 Step 3: calculation of the parameters needed for the mathematical model

The intercellular transport costs are shown in Table 5. The rest of the parameters to be used in the mathematical model were the following: average cycle time C ¯ = 0.15 gives a value of S min = 6, and a value of S max = 10 was also considered. The number of cells obtained from the combined precedence diagram gave values of K min = 4 and K max = 6, which are adequate for the groups of machines that can be visualized previously, giving values of M max = 8 and M min = 5, respectively. All this information is presented in Tables 68.

i′ i
cii’ (MM$)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
2 0.123
3 0 0
4 0.093 0 0
5 0 0 0.165 0
6 0 0 0 0.075 0
7 0.123 0.123 0 0 0 0
8 0 0 0 0 0 0 0
9 0 0 0.165 0 0.165 0 0 0
10 0 0 0.308 0 0.165 0 0 0 0.165
11 0.353 0 0 0.168 0 0.075 0 0 0 0
12 0 0 0 0 0 0 0 0.165 0 0 0
13 0 0 0.308 0 0.165 0 0 0 0.165 0.308 0 0
14 0.243 0 0 0.168 0 0.075 0 0 0 0 0.318 0 0
15 0 0 0 0 0 0 0 0 0 0 0 0 0 0
16 0 0 0 0 0 0 0 0.165 0 0 0 0.165 0 0 0
17 0.123 0.123 0 0 0 0 0.123 0 0 0 0 0 0 0 0
18 0 0 0.165 0.075 0.165 0.225 0 0 0.165 0.165 0.075 0 0.165 0.075 0
19 0 0 0.143 0 0 0 0 0 0 0.143 0 0 0.143 0 0.058
20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.058
21 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.058
23 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
24 0 0 0 0 0 0 0 0 0 0 0 0.11 0 0 0
25 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.058
26 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.058
27 0 0 0 0 0 0.15 0 0 0 0 0 0 0 0 0
28 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.058
29 0 0 0 0 0 0.15 0 0 0 0 0 0 0 0 0
30 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
31 0 0 0 0 0 0 0 0 0 0 0 0.11 0 0 0
16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
17 0
18 0 0
19 0 0 0
20 0 0 0 0.253
21 0 0 0.203 0 0
22 0 0 0 0.185 0.1 0
23 0 0 0.203 0 0 0.203 0
24 0 0 0 0 0 0 0 0
25 0 0 0 0.185 0.1 0 0.185 0 0
26 0 0 0 0.338 0.253 0 0.185 0 0 0.185
27 0 0 0.353 0 0 0.203 0 0.203 0 0 0
28 0 0 0 0.338 0.253 0 0.185 0 0 0.185 0.338 0
29 0 0 0.353 0 0 0.203 0 0.203 0 0 0 0.353 0
30 0 0 0 0.043 0.043 0 0.043 0 0 0.043 0.043 0 0.043 0
31 0 0 0 0 0 0 0 0 0.11 0 0 0 0 0 0

Table 5.

Intercellular transport costs between machines i and i´ (Cii´).

i 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
t ¯ i 0.033 0.03 0.03 0.031 0.033 0.034 0.031 0.021 0.02 0.031 0.034 0.028 0.029 0.03 0.014
i 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
t ¯ i 0.022 0.031 0.034 0.028 0.023 0.034 0.034 0.034 0.014 0.02 0.034 0.034 0.026 0.033 0.011

Table 6.

Weighted average processing times ( t ¯ i ) for each machine i.

K min K max θ M max M min M C ¯ A S min S max cs (MM$) cb (MM$)
4 6 4 8 5 31 0.15 0.005 6 10 0.45 0.05

Table 7.

Parameters to be used in the model.

E (8,12), (24,12), (12,16), (12,31), (1,2), (1,4), (2,7), (4,11), (7,17), (11,14), (18,4), (18,6), (18,21), (18,23), (18,5), (6,11), (6,27), (21,27), (23,27), (27,29), (5,10), (3,10), (3,9), (10,13), (9,13), (19,3), (19,22), (19,20), (22,26), (20,26), (20,25), (25,30), (25,15), (26,28), (30,28), (15,28)

Table 8.

Set of precedence relations between machines.

3.4 Step 4: statement and solution of the mathematical model

With the parameters calculated above, we state the mathematical model. The model was then solved using the Extended LINGO© version 8.0 software, getting an optimal global solution in 2 h and 32 min, after 121 iterations, giving as a result 5 manufacturing cells and 9 work stations, 3 of them multicellular, because they process exceptional products. This solution is represented in Figure 7.

Figure 7.

Representation of the resultant cells and stations.

3.5 Step 5: assigning the product models to the production cells

With the groups of machines obtained thanks to the model’s solution, the product models were assigned to each of these production cells. This solution is represented in Table 9.

Table 9.

Assignment of products to cells.

The proposed cellular manufacturing system that results from the application of the methodology shows the resultant cells and the product families assigned to them, where it is seen that cell 1 processes products of the beam and slotted angle type, cells 2 and 3 process only accessories, cell 4 processes only beams, and cell 5 processes only pillars. In this way, it is easier to improve the obtained balance by sequencing mixed product models, because in a cell there are no different types of products that compete for the same resources (machines), but they rather process the same types of products, grouped in families. Only work stations 1, 2, and 8 turn out being multicellular, because they process exceptional products that undergo intercellular movements, most of which are associated with the operation of the cutter type machines, since they are those that are mostly shared by the different products.

To analyze the results, it is necessary to have performance measures of the solution, but since the present study has taken up a new problem that is part of the manufactured cell formation problem (MCFP) and the general assembly line balancing problem (GALBP), we must use measures of performance commonly employed for both problems separately, such as group capability index (GCI), for example, proposed by Seifoddini and Hsu [23], and the grouping efficacy (GE) proposed by Kumar and Chandrasekharan [24], as well as balancing measures like “total line imbalance” proposed by Thomopoulos [12]. Table 10 presents a comparison of results of the case study by the proposed methodology versus the model of Won and Currie [25], considering the latter as an MCFP-type problem. The model of [25] was chosen because it considers a number of production factors and a complexity level relatively similar to the proposed methodology, and it can be solved in Lingo© in a reasonable computation time.

Problem number (m × n) Computation time (min) Number of iterations GCI (%) GE (%)
Proposed methodology Won and Currie [23] Proposed methodology Won and Currie [23] Proposed methodology Won and Currie [23] Proposed methodology Won and Currie [23]
(31 × 67) 152 140 121 114 89.45 85.45 86.45 84.25

Table 10.

Comparison of the measures for the MCFP.

In Table 10, it is seen that the proposed methodology takes more time to get the global optimal solution as well as a greater number of iterations compared to the model of [25]. This can be explained because in the proposed model of the methodology the first component of the formulation is of a quadratic type, so the LINGO© software uses an algorithm specialized in these types of problems which makes each iteration take longer time; on the other hand, the model of [25] is of the p-media type, which although linear, in terms of the quality of the solution, it is seen that the proposed methodology gives a group capability index (GCI) greater than the model of [25]. We believe that this is due to the fact that the proposed formulation puts emphasis in minimizing the costs of intercellular product movements, in contrast with the model of Won and Currie [25], which aims to maximize the similarity between the machines that constitute the production cells. On the other hand, with respect to the grouping efficacy (GE), the proposed methodology is better than the model of [25]. Making a deeper analysis of the results and getting more general conclusions is risky because it is necessary to do further research with study cases in which the indicators and models used can be compared. In relation to the comparison of indicators such as the total imbalance of the line proposed by Thomopoulos [12] under the viewpoint of the general assembly line balancing problem (GALBP), it was not done, and it is expected to be dealt with in an extension of the present work.


4. Conclusions

This work has proposed a new way of approaching balancing and cell formation problems, which were previously studied independently. In this way, it is possible to consider aspects that were previously avoided, such as production volumes, processing times, and operation sequences for the MCFP, and the fact that the production cells are not established yet and that they also share information among themselves, for balancing the lines. In the proposed methodology, a model is presented that has advantages at the time of solving it with a commercial software, because it does not need as many variables as other proposals. This approach delivers cells that are more amenable to be used in practice, although it will remain for future research to deliver immediately the layout of the cells as well as to integrate the problem of sequencing mixed models and, in that way, to improve the balance of each station.

In relation to the preliminary results obtained by comparing with another mathematical model from the viewpoint of a type MCFP problem, it must be pointed out that although the grouping efficacy indicators are better, they are not conclusive if a comparative study with more study problems as well as a sensitivity analysis of the parameters of the proposed model is not made. The same comparative study must be made from the viewpoint of a type GALBP problem in relation to the balancing of the resultant line, in such a way that the indicators can be compared with the proposed methodology with the techniques or methods applied separately for the same sets of study problems.

Finally, we believe that the proposed methodology responds to a problem of integrating two problems like the MCFP and GALBP under a same approach which is perfectible insofar as the results can be validated in future comparative work, as well as extending the proposed problem, integrating in the design, and balancing of manufacturing cells configured in “U” the design of the family of products that one wishes to manufacture.



The authors would like to thank for the support of the Universidad de Santiago de Chile (031817 PB DICYT-USACH) and also thank all the teams of the Smart City Lab Program Centre of the same university.


  1. 1. Selim HM, Askin RG, Vakharia AJ. Cell formation in group technology: Evaluation and directions for future research. Computers and Industrial Engineering. 1998;34(1):3-20. DOI: 10.1016/S0360-8352(97)00147-2
  2. 2. Deuse J, Konrad B, Bohnen F. Renaissance of group technology: Reducing variability to match lean production prerequisite. In: Proceedings of the 7th IFAC Conference on Manufacturing Modelling, Management, and Control International Federation of Automatic Control (IFAC); Saint Petersburg, Russia, 19–21 June. 2013. pp. 998-1003
  3. 3. Papaioannou G, Wilson J. The evolution of cell formation problem methodologies based on recent studies (1997–2008): Review and directions for future research. European Journal of Operational Research. 2010;206:509-521. DOI: 10.1016/j.ejor.2009.10.020
  4. 4. Kashan A, Karimi B, Noktehdan A. A novel discrete particle swarm optimization algorithm for the manufacturing cell formation problem. The International Journal of Advanced Manufacturing Technology. 2014;73(9–12):1543-1556. DOI: 10.1007/s00170-014-5906-4
  5. 5. Karoum B, Elbenani YB. Discrete cuckoo search algorithm for solving the cell formation problem. International Journal of Manufacturing Research. 2019;14(3). DOI: 10.1504/IJMR.2019.100991
  6. 6. Patel J, Patel S. Approaches to solve cell formation, machine layout and cell layout problem: A review. Transactions on Machine Learning and Artificial Intelligence. 2014;2(5):80-96
  7. 7. Anbumalar V, Sekar R. Methods for solving cell formation, static layout and dynamic layout cellular manufacturing system problems: A review. Asian Journal of Science and Technology. 2015;06(12):2107-2112
  8. 8. Kumar S, Sharma R. Cell formation heuristic procedure considering production data. International Journal of Production Management and Engineering. 2014;2(2):75-84. DOI: 10.4995/ijpme.2014.2078
  9. 9. Sparling D. Balancing JIT production units: The N U-line balancing problem. Information Systems and Operational Research. 1998;36(4):215-237. DOI: 10.1080/03155986.1998.11732360
  10. 10. Miltenburg J. Balancing U-lines in a multiple U-line facility. European Journal of Operational Research. 1998;109(1):1-23. DOI: 10.1016/S0377-2217(97)00169-0
  11. 11. Sparling D, Miltenburg GJ. The mixed-model U-line balancing problem. International Journal of Production Research. 1998;36(2):485-501. DOI: 10.1080/002075498193859
  12. 12. Thomopoulos N. Line balancing-sequencing for mixed-model assembly. Management Science. 1967;14(2):B59-B75
  13. 13. Mustafa F, Kursad A, Mustafa Y. A new algorithm for U-shaped two-sided assembly line. Transactions of the Canadian Society for Mechanical Engineering. 2010;34(2):225-241
  14. 14. Fathi M, Alvarez M, Rodriguez VA. New heuristic approach to solving U-shape assembly line balancing problems type. World Academy of Science Engineering and Technology. 2011;59:413-421
  15. 15. Turkay AF. On the MILP-model for the U-shaped assembly line balancing problems. European Journal of Operational Research. 2015;242:343-346. DOI: 10.1016/j.ejor.2014.10.036
  16. 16. Chutima P, Suchanun T. Productivity improvement with parallel adjacent U-shaped assembly lines. Advances in Production Engineering & Management. 2019;14(1):51-64. DOI: 10.14743/apem2019.1.311
  17. 17. Wemmerlöv U, Hyer N. Cellular manufacturing in the US industry: A survey of users. International Journal of Production Research. 1988;27(9):1511-1530. DOI: 10.1080/00207548908942637
  18. 18. Macaskill JLC. Production-line balances for mixed-model lines. Management Science. 1972;19(4):423-434
  19. 19. Ahmadi RH, Wurgaft H. Design for synchronized flow manufacturing. Management Science. 1994;40(11):1469-1483
  20. 20. Urban TL. Optimal balancing of U-shaped assembly lines. Management Science. 1998;44(5):738-741
  21. 21. Al Kattan I. Workload balance of cells in designing of multiple cellular manufacturing systems. Journal of Manufacturing Technology Management. 2005;16(2):178-196. DOI: 10.1108/17410380510576822
  22. 22. Burbidge JL. The Introduction to Group Technology. Nueva York: John Wiley & Sons; 1975 267p
  23. 23. Seifoddini H, Hsu CP. Comparative study of similarity coefficients and clustering algorithms in cellular manufacturing. Journal of Manufacturing Systems. 1994;13(2):119-127
  24. 24. Kumar KR, Chandrasekharan MP. Grouping efficacy: A quantitative criterion for goodness of block diagonal forms of binary matrices in group technology. International Journal of Production Research. 1990;28(2):233-243. DOI: 10.1080/00207549008942706
  25. 25. Won Y, Currie K. An effective p-median model considering production factors in machine cell-part family formation. Journal of Manufacturing Systems. 2006;25(1):58-64. DOI: 10.1016/S0278-6125(06)80033-6

Written By

Luis Valdivia and Pedro Palominos

Submitted: 11 May 2019 Reviewed: 02 September 2019 Published: 27 November 2019