# An Optimization Method for the Remanufacturing Dynamic Facility Layout Problem with Uncertainties

Authors: Lingling Li; Congbo Li; Huijie Ma; Ying Tang
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

Remanufacturing is a practice of growing importance due to increasing environmental awareness and regulations. Facility layout design, as the cornerstone of effective facility planning, is concerned about resource localization for a well-coordinated workflow that leads to lower material handling costs and reduced lead times. However, due to stochastic returns of used products/components and their uncontrollable quality conditions, the remanufacturing process exhibits a high level of uncertainty challenging the facility layout design for remanufacturing. This paper undertakes this problem and presents an optimization method for remanufacturing dynamic facility layout with variable process capacities, unequal processing cells, and intercell material handling. A dynamic multirow layout model is presented for layout optimization and a modified simulated annealing heuristic is proposed toward the determination of optimal layout schemes. The approach is demonstrated through a machine tool remanufacturing system.


---

## Body

## 1. Introduction

Remanufacturing, as an industrial process of restoring discarded products/components back to their useful lives, is of growing importance due to the emerging pressure of legislation and increasing awareness of environmental conservation. As the ultimate form of recycling, remanufacturing maintains much of the value added from original manufacturing to material and reprocessing conservation, leading to lower production costs and improved firm profits [1]. A recent Wall Street Journal revealed that large scale remanufacturing in the United States employs more than 500,000 people and contributes to approximately $100 billion of goods sold each year [2]. More successful industry examples can be seen in Kodak, BMW, IBM, DEC, and Xerox [3]. However compared to manufacturing, remanufacturing is more complicated in the way that (1) the supply of returned products is unpredictable in timing and quantities; (2) the quality and composition of returned products vary; and (3) the process routings are not necessarily fixed but rather adapt to the actual conditions of products/components [4]. Such characteristics have led the facility layout problem (FLP) in remanufacturing to be much more complicated than that in traditional manufacturing.Facility layout problem (FLP), as the cornerstone of effective facility planning, is concerned about resource localization for a well-coordinated workflow that directly or indirectly leads to lower material handling costs, less work-in-process inventory, reduced lead times, and so forth [5]. Research in this area is commonly categorized as a static or dynamic FLP, depending on the nature of the input requirements and the time periods under consideration. Chaieb and Korbaa [6] dealt with the intracell machine layout problem in a circular bidirectional configuration to minimize both the movement distance and total transportation time. Chan et al. [7] incorporated the machine-part grouping decision into the static cellular layout problem to minimize the intercellular movement costs. Such efforts on static layout are usually performed for a single time period and assume flow between machines, product demand, and level of product mix constant.Given the nature of manufacturing environments where uncertainty is inevitable, periodic layout reconfiguration with respect to changes is necessary as the best facility layout configuration for one time period may not be efficient for successive others. A relayout or dynamic facility layout problem (DFLP) is usually modeled by discretizing the time into several planning periods where the localization of an existing facility in each period is a decision variable. Dharmalingam et al. [8] proposed a cellular layout method for a manufacturing system with dynamic demands. With the assumptions of equal size cells and a fixed number of machines per cell, their objective is to minimize the number of intercell movements and to perform greatest proportion of part operations within a single cell. Sharing the same assumption that all departments have the same size, Baykasoglu et al. [9] developed a dynamic facility layout model to minimize both material handling and reconfiguration costs. In the authors’ own claim, it was the first time that budget constraints were considered in DFLP. McKendall Jr. and Hakobyan [10] tackled DFLP with unequal-sized, free oriental departments. However, their design focused on open-field layout design with much loose constraints. By integrating the decisions of cell formation and group layout together, Kia et al. [11] developed a dynamic layout model that was incorporated with several design features, such as alternative process routings, operation sequences, machine capacity, lot splitting, multirow layout of equal area facilities, and flexible reconfigurations. The same group extended the work to propose a multiobjective model with more optimization constraints, including machine capability and capacity, part demand satisfaction, cell size, and location assignment [12].Summarizing the findings of the above literatures, few efforts have been made towards the optimization of DFLP with consideration of remanufacturing natures. In fact, the high level of uncertainty in remanufacturing challenges its layout design to account for variations. For instance, the capacity of each remanufacturing process, the dimension of processing cells, and the scales of material handing between two cells are nonstationary but change with time. Thus, there is a strong case for focusing on DFLP in remanufacturing. In this paper, an attempt has been made to investigate such a problem. In particular, a dynamic multirow layout problem is studied with consideration of variable process capacities, unequal processing cells, and intercell material handling. The rest of the paper is organized as follows. Section2 analyzes the uncertainties of DFLP in remanufacturing. Section 3 presents the dynamic facility layout optimization model. Section 4 presents a modified simulated annealing heuristic for the determination of optimal layout schemes. An example is provided in Section 5, followed by the conclusion in Section 6.


## 2. Uncertainties of FLP in Remanufacturing

In a remanufacturing system, there are a number of uncertain characteristics that significantly complicate its remanufacturing activities [4, 13]. Firstly, the product return process is highly uncertain with respect to timing, when the products are available for remanufacturing, and quantity, how many products are available. The uncertain characteristic in timing and quantity of returns is a reflection of the uncertain nature of the life of a product, which is commonly influenced by a number of factors including the life-cycle stage of a product and the rate of technological changes [14].Product/component returned in remanufacturing exhibits highly uncontrolled variability with respect to product condition, ranging from slightly used with minor blemishes to significantly damaged and requiring extensive repair [14]. In a remanufacturing facility, used products are first manually disassembled into what is called “cores” that are disassembled parts or components. Then, inspection and sorting of returns are usually implemented in the system to identify core’s quality and organize necessary activities for their recovery. According to inspection results, the cores are then classified into three groups (i.e., “reusable,” “remanufacturable,” and “material recoverable”). For the remanufacturable cores, their process routings are not fixed but rather adapted to their specific condition [15–17]. Figure 1 gives an example of the remanufacturing routings for cores (e.g., spindles, guide ways, worm gears, and lead screws) disassembled from a used machine tool. For instance, if a spindle is detected with severe abrasion, two typical process routings would be “grinding → chromium electroplating” or “grinding → cold welding → fine grind” to regain the surface accuracy. Among these routings, some remanufacturing operations may be probabilistic since the quality attributes of returned cores are characterized in the form of a set of classification probabilities [15].Figure 1
The remanufacturing process routings of machine tool cores.Another uncertain characteristic associated with remanufacturing is the high degree of variability in remanufacturing processing times. For a failed core, the amount of time required to perform the necessary remanufacturing operation is strongly contingent upon the core’s quality condition. In particular, the processing time is a function of the core, causing the core’s failure, and the extent of the failure, which is unpredictable in prior.From a facility layout perspective, the high level of uncertainty in remanufacturing challenges its layout design to account for variations. For instance, the remanufacturing processes are subject to extreme fluctuations in workloads because of uncertainty in stochastic returns, probabilistic process routings, and highly variable processing times. Accordingly, the required processing capacity of each remanufacturing process (i.e., the number of machines required to perform the specific remanufacturing operation) is nonstationary but changes with time. As a result, the dimensions of individual processing cells to accommodate a certain number of machines vary as well. Besides, due to the stochastic returns of used components and their stochastic process routings, the frequency of intercell material handling between two remanufacturing processing cells varies constantly over time. Given so much stochastic nature in remanufacturing, the best facility layout configuration obtained in a given time period may not be efficient for the successive others. Therefore, there is a strong case for focusing on DFLP in remanufacturing by discretizing the time into several planning periods. This paper undertakes this challenge and incorporates such variability in the optimization of a remanufacturing DFLP (R-DFLP) with the aim of achieving overall high efficiency and minimum cost.


## 3. R-DFLP Optimization Model

### 3.1. Problem Statement

A typical remanufacturing facility considered in the paper is a multirow shop floor, each row of which has a set of designated departments with an equal size. The facility accommodates a finite set of remanufacturing processes to restore components of a specific type of used products. A processing cell (i.e., grouping a number of machines to perform a specific remanufacturing operation) is dedicated to each remanufacturing process, which occupies one or more than one departments according to its dimension to be determined over time. Figure2 gives an example of a facility layout configuration with 3 rows, 27 departments, and 6 processing cells. The problem is formulated under the following assumptions.(a)
A dynamic facility layout in remanufacturing is considered by discretizing the time into several planning periods(t=1,2,…,T). That is, the remanufacturing facility is allowed to be reconfigured at each given time period according to the dynamics of the system.(b)
The required capacity of a processing cell is not fixed but adapted to its workload in a given unit of time. Thus, the number of machines assigned to the cell might be different from one period to another. WhenNit<Ni(t-1), the redundant machine would be removed from the cell. When Nit>Ni(t-1), additional machines will be added to the ith processing cell. In either case, relocation cost would be charged for machine equipment handling regarding machine removal, machine transiting, and machine installation.(c)
There are an infinite number of machines in the facility. The purchase of new machines is not necessary and considered.(d)
The processing cells in the same row have the sameY-coordinate.(e)
Each processing cell only occupies departments in the same row.(f)
The dimension of the shop floor is given.(g)
The dimensions of machines are given but vary with respect to the processing cells they serve (i.e., the machines in a processing cell have the same dimensions).(h)
The intracell material handling cost is negligible.(i)
The width of a department is always larger than the length of machines in any processing cell.(j)
The orientations of machines in a processing cell are kept the same.Figure 2
An example of a facility layout configuration with 3 rows, 27 departments, and 6 processing cells.


### 3.2. Notations and Variables Involved in R-DFLP

The notations and variables used in the model are given below.


#### 3.2.1. Index

t: time period index (t=1,2,…,T), where T is the total number of periods;
i
,
j: processing cell index (i,j=1,2,…,Z), where Z is the total number of processing cells in the shop floor; in particular, i=1,2,3, Z denote disassembly, cleaning, inspection, and reassembly, respectively, while i=4,…,Z-1 represent the remanufacturing operations;
k: core type index;
r: remanufacturing processing routing index;
f: facility row index (f=1,2,…,F), where F is the number of rows in the shop floor;
h: department index in a given row (h=1,2,…,H), where H is the number of departments per row.


#### 3.2.2. Model Parameters

D
t: time duration (in days) of the period t;
Q
k
t: number of the cores of type k in period t;
p
k
r
t: probability of the cores of type k being assigned to the rth processing routing in period t;
c
i: machine availability (i.e., the available machine hours per workday of an individual machine) of the ith processing cell (hours/day);
ϖ
i: machine efficiency (in percentage) of the ith processing cell corresponding to machine failures;
ψ
i
k: average processing time (in minutes) for the cores of type k in the ith processing cell;
M
l
i: length of the machines in the ith processing cell (Mli≥Mwi);
M
w
i: width of the machines in the ith processing cell;
S
l: length of the remanufacturing facility (Sl≥Sw);
S
w: width of the remanufacturing facility;
ε: the width of main road designed in the facility;
δ
i: relocation cost for the machines in the ith processing cell;
θ
k: material handling cost of cores of type k per unit of distances;
Δ
t: total material handling cost in period t;
B
t: relocation budget in period t;
Γ
i
j
k
r: binary variable that has a value of 1 if the processing cell i precedes the processing cell j in the rth processing routing for the cores of type k and 0 otherwise;
α
i
k
r: binary variable that has a value of 1 if the ith processing cell is in the rth remanufacturing processing routing for the cores of type k and 0 otherwise.Among the above-mentioned parameters, there are three uncertain parameters that significantly complicate the R-DFLP and need to be elaborated as follows.(1) The Arrivals of Used Products. In a remanufacturing system, arrivals (Qkt) related to the failures of products are uncertain in timing and quantities and are commonly modelled as a Poison process [9].(2) The Process Routings. The process routings are not necessarily fixed but rather adapt to the actual conditions of cores. The likelihood (Pkrt) of a remanufacturing process routing being assigned to a certain type of cores is randomly generated to represent the stochastic features [4].(3) The Processing Times. The processing times (ψik) required to perform the necessary remanufacturing operations for a certain part are highly variable. An exponential distribution is used for the processing times in order to simulate the large range of possible values [9].


#### 3.2.3. Decision Variables

N
i
t: number of machines in the ith processing cell in period t;
β
i
t: binary variable that has a value of 1 if the machines in the ith processing cell are oriented vertically and 0 otherwise;
L
i
t: length of the ith processing cell in period t;
η
i
t: number of departments occupied by the ith processing cell in period t;
π
i
h
f
t: binary variable that has a value of 1 if the ith processing cell occupies the department in the hth row and ith column as its first location in period t;
(
x
i
t
,
y
i
t
): the centroid of the ith processing cell in period t;
d
i
j
t: distance (in meters) between the two processing cells i and j in period t;
τ
i
t: binary variable that has a value of 1 if the centroid of the processing cell i in period t is different from that in period t-1 and 0 otherwise.As stated earlier, the variability in remanufacturing regarding stochastic returns, probabilistic process routings, and highly variable processing times complicates its facility layout design. Thus the decision variables in the R-DFLP, including the number of machines, the dimensions of individual processing cells, and the detailed location (i.e., the centroid) of each processing cell in a given time period, will be mathematically elaborated in the following.(1) Number of Machines in Processing Cells. The number of machines dedicated to the ith processing cell in period t, Nit, is a function of the number of cores needed to be processed at the time, the average processing times, and machine availability and efficiency:
(1)Nit=∑k∑rQkt·pkrt·ψik·αikrϖi·ci·Dt,
where ci denotes the machine availability (i.e., the available machine hours per workday of a certain machine type) of the ith processing cell (hours/day) and ϖi is the machine efficiency (in percentage) of the ith processing cell corresponding to machine failures.(2) Dimensions of Processing Cells. It is assumed that an individual processing cell should be equipped with a series of machines of the same type to perform a specific remanufacturing process. The dimension of a processing cell as formulated in (2) is pertinent to the numbers of machines dwelt in and their orientations. The number of departments occupied by the ith processing cell is given in (3):
(2)Lit=Mli·1-βit+Mwi·βit·Nit,(3)ηit=ceilLitSl/H,
where Sl/H (i.e., dividing Sl by H) gets the length of each department designed in each row of the floor shop and ceiling (x) returns a value that is the smallest integer not less than x.(3) The Centroids of Processing Cells. For each processing cell, with the required number of departments (ηit) and the first department it occupies (πihft) being determined, its centroids in period t then can be determined by (4). The distance between the centroids of the processing cell i and j in period t is given in (5):
(4)xit=h-0.5·Sl2H+ηit-1·Sl2H,yit=Sw2F+f-1·SwF,xxxxxxxxxxxx∀πihft=1,(5)dijt=xit-xjt+yit-yjt.


### 3.3. Optimization Model

The R-DFLP considered in this paper aims to balance the trade-off between the intercell material handling cost and the increased relocation cost of inefficient layouts.(1) Cost of Intercell Material Handing. While cores go through various processes to be restored back to their useful lives, material handling cost is associated with their logistic activities (i.e., moving from one processing cell to another). The total amount of material handling cost in a given unit of time is an incremental linear function of the total distances of cores moving in the facility, as defined below:
(6)Δt=∑k∑r∑i∑jQkt·pkrt·Γijkr·dijt·θk,
where θk denotes the material handling cost of cores of type k per unit distance.(2) Relocation Cost. Since the best facility layout scheme for one time period may not be efficient for the successive others, the floor shop needs to be reconfigured accordingly by changing the locations of a series of processing cells in a given time period. In this case, machine equipment in these processing cells needs to be removed, transited, and reinstalled which necessarily incur a certain amount of relocation costs, as formulated in the following:(7)Ωt=∑iτit·Nit·δi,∀t>1,τit=0,ifxit=xi(t-1),yit=yi(t-1),1,otherwise,∀t>1,
where τit is a binary variable that has a value of 1 if the centroid of the processing cell i in period t is different from that in period t-1 and 0 otherwise; δi is the relocation cost per machine of the ith processing cell.The R-DFLP problem is then formulated as a nonlinear mixed-integer programming model to minimize the total cost of intercell material handling and machine relocation:(8)Minimizetotal cost=∑t=1TΔt+∑t=2TΩt(9)Subject  to∑t=2T∑iτit·δi≤Bt(10)∑iηit≤F·H,∀t,(11)∑i=3∑j=4Z-1∑k∑rQkt·pkrt·Γijkr=∑i=4Z-1∑j=Z∑k∑rQkt·pkrt·Γijkr,∀t.
Constraint (9) assures that the total relocation cost in any period does not exceed the given budget. Constraint (10) guarantees that the total number of departments occupied by processing cells is not larger than the total available in the shop floor. Constraint (11) represents material flow conservation between disassembly and reassembly.


## 3.1. Problem Statement

A typical remanufacturing facility considered in the paper is a multirow shop floor, each row of which has a set of designated departments with an equal size. The facility accommodates a finite set of remanufacturing processes to restore components of a specific type of used products. A processing cell (i.e., grouping a number of machines to perform a specific remanufacturing operation) is dedicated to each remanufacturing process, which occupies one or more than one departments according to its dimension to be determined over time. Figure2 gives an example of a facility layout configuration with 3 rows, 27 departments, and 6 processing cells. The problem is formulated under the following assumptions.(a)
A dynamic facility layout in remanufacturing is considered by discretizing the time into several planning periods(t=1,2,…,T). That is, the remanufacturing facility is allowed to be reconfigured at each given time period according to the dynamics of the system.(b)
The required capacity of a processing cell is not fixed but adapted to its workload in a given unit of time. Thus, the number of machines assigned to the cell might be different from one period to another. WhenNit<Ni(t-1), the redundant machine would be removed from the cell. When Nit>Ni(t-1), additional machines will be added to the ith processing cell. In either case, relocation cost would be charged for machine equipment handling regarding machine removal, machine transiting, and machine installation.(c)
There are an infinite number of machines in the facility. The purchase of new machines is not necessary and considered.(d)
The processing cells in the same row have the sameY-coordinate.(e)
Each processing cell only occupies departments in the same row.(f)
The dimension of the shop floor is given.(g)
The dimensions of machines are given but vary with respect to the processing cells they serve (i.e., the machines in a processing cell have the same dimensions).(h)
The intracell material handling cost is negligible.(i)
The width of a department is always larger than the length of machines in any processing cell.(j)
The orientations of machines in a processing cell are kept the same.Figure 2
An example of a facility layout configuration with 3 rows, 27 departments, and 6 processing cells.


## 3.2. Notations and Variables Involved in R-DFLP

The notations and variables used in the model are given below.


### 3.2.1. Index

t: time period index (t=1,2,…,T), where T is the total number of periods;
i
,
j: processing cell index (i,j=1,2,…,Z), where Z is the total number of processing cells in the shop floor; in particular, i=1,2,3, Z denote disassembly, cleaning, inspection, and reassembly, respectively, while i=4,…,Z-1 represent the remanufacturing operations;
k: core type index;
r: remanufacturing processing routing index;
f: facility row index (f=1,2,…,F), where F is the number of rows in the shop floor;
h: department index in a given row (h=1,2,…,H), where H is the number of departments per row.


### 3.2.2. Model Parameters

D
t: time duration (in days) of the period t;
Q
k
t: number of the cores of type k in period t;
p
k
r
t: probability of the cores of type k being assigned to the rth processing routing in period t;
c
i: machine availability (i.e., the available machine hours per workday of an individual machine) of the ith processing cell (hours/day);
ϖ
i: machine efficiency (in percentage) of the ith processing cell corresponding to machine failures;
ψ
i
k: average processing time (in minutes) for the cores of type k in the ith processing cell;
M
l
i: length of the machines in the ith processing cell (Mli≥Mwi);
M
w
i: width of the machines in the ith processing cell;
S
l: length of the remanufacturing facility (Sl≥Sw);
S
w: width of the remanufacturing facility;
ε: the width of main road designed in the facility;
δ
i: relocation cost for the machines in the ith processing cell;
θ
k: material handling cost of cores of type k per unit of distances;
Δ
t: total material handling cost in period t;
B
t: relocation budget in period t;
Γ
i
j
k
r: binary variable that has a value of 1 if the processing cell i precedes the processing cell j in the rth processing routing for the cores of type k and 0 otherwise;
α
i
k
r: binary variable that has a value of 1 if the ith processing cell is in the rth remanufacturing processing routing for the cores of type k and 0 otherwise.Among the above-mentioned parameters, there are three uncertain parameters that significantly complicate the R-DFLP and need to be elaborated as follows.(1) The Arrivals of Used Products. In a remanufacturing system, arrivals (Qkt) related to the failures of products are uncertain in timing and quantities and are commonly modelled as a Poison process [9].(2) The Process Routings. The process routings are not necessarily fixed but rather adapt to the actual conditions of cores. The likelihood (Pkrt) of a remanufacturing process routing being assigned to a certain type of cores is randomly generated to represent the stochastic features [4].(3) The Processing Times. The processing times (ψik) required to perform the necessary remanufacturing operations for a certain part are highly variable. An exponential distribution is used for the processing times in order to simulate the large range of possible values [9].


### 3.2.3. Decision Variables

N
i
t: number of machines in the ith processing cell in period t;
β
i
t: binary variable that has a value of 1 if the machines in the ith processing cell are oriented vertically and 0 otherwise;
L
i
t: length of the ith processing cell in period t;
η
i
t: number of departments occupied by the ith processing cell in period t;
π
i
h
f
t: binary variable that has a value of 1 if the ith processing cell occupies the department in the hth row and ith column as its first location in period t;
(
x
i
t
,
y
i
t
): the centroid of the ith processing cell in period t;
d
i
j
t: distance (in meters) between the two processing cells i and j in period t;
τ
i
t: binary variable that has a value of 1 if the centroid of the processing cell i in period t is different from that in period t-1 and 0 otherwise.As stated earlier, the variability in remanufacturing regarding stochastic returns, probabilistic process routings, and highly variable processing times complicates its facility layout design. Thus the decision variables in the R-DFLP, including the number of machines, the dimensions of individual processing cells, and the detailed location (i.e., the centroid) of each processing cell in a given time period, will be mathematically elaborated in the following.(1) Number of Machines in Processing Cells. The number of machines dedicated to the ith processing cell in period t, Nit, is a function of the number of cores needed to be processed at the time, the average processing times, and machine availability and efficiency:
(1)Nit=∑k∑rQkt·pkrt·ψik·αikrϖi·ci·Dt,
where ci denotes the machine availability (i.e., the available machine hours per workday of a certain machine type) of the ith processing cell (hours/day) and ϖi is the machine efficiency (in percentage) of the ith processing cell corresponding to machine failures.(2) Dimensions of Processing Cells. It is assumed that an individual processing cell should be equipped with a series of machines of the same type to perform a specific remanufacturing process. The dimension of a processing cell as formulated in (2) is pertinent to the numbers of machines dwelt in and their orientations. The number of departments occupied by the ith processing cell is given in (3):
(2)Lit=Mli·1-βit+Mwi·βit·Nit,(3)ηit=ceilLitSl/H,
where Sl/H (i.e., dividing Sl by H) gets the length of each department designed in each row of the floor shop and ceiling (x) returns a value that is the smallest integer not less than x.(3) The Centroids of Processing Cells. For each processing cell, with the required number of departments (ηit) and the first department it occupies (πihft) being determined, its centroids in period t then can be determined by (4). The distance between the centroids of the processing cell i and j in period t is given in (5):
(4)xit=h-0.5·Sl2H+ηit-1·Sl2H,yit=Sw2F+f-1·SwF,xxxxxxxxxxxx∀πihft=1,(5)dijt=xit-xjt+yit-yjt.


## 3.2.1. Index

t: time period index (t=1,2,…,T), where T is the total number of periods;
i
,
j: processing cell index (i,j=1,2,…,Z), where Z is the total number of processing cells in the shop floor; in particular, i=1,2,3, Z denote disassembly, cleaning, inspection, and reassembly, respectively, while i=4,…,Z-1 represent the remanufacturing operations;
k: core type index;
r: remanufacturing processing routing index;
f: facility row index (f=1,2,…,F), where F is the number of rows in the shop floor;
h: department index in a given row (h=1,2,…,H), where H is the number of departments per row.


## 3.2.2. Model Parameters

D
t: time duration (in days) of the period t;
Q
k
t: number of the cores of type k in period t;
p
k
r
t: probability of the cores of type k being assigned to the rth processing routing in period t;
c
i: machine availability (i.e., the available machine hours per workday of an individual machine) of the ith processing cell (hours/day);
ϖ
i: machine efficiency (in percentage) of the ith processing cell corresponding to machine failures;
ψ
i
k: average processing time (in minutes) for the cores of type k in the ith processing cell;
M
l
i: length of the machines in the ith processing cell (Mli≥Mwi);
M
w
i: width of the machines in the ith processing cell;
S
l: length of the remanufacturing facility (Sl≥Sw);
S
w: width of the remanufacturing facility;
ε: the width of main road designed in the facility;
δ
i: relocation cost for the machines in the ith processing cell;
θ
k: material handling cost of cores of type k per unit of distances;
Δ
t: total material handling cost in period t;
B
t: relocation budget in period t;
Γ
i
j
k
r: binary variable that has a value of 1 if the processing cell i precedes the processing cell j in the rth processing routing for the cores of type k and 0 otherwise;
α
i
k
r: binary variable that has a value of 1 if the ith processing cell is in the rth remanufacturing processing routing for the cores of type k and 0 otherwise.Among the above-mentioned parameters, there are three uncertain parameters that significantly complicate the R-DFLP and need to be elaborated as follows.(1) The Arrivals of Used Products. In a remanufacturing system, arrivals (Qkt) related to the failures of products are uncertain in timing and quantities and are commonly modelled as a Poison process [9].(2) The Process Routings. The process routings are not necessarily fixed but rather adapt to the actual conditions of cores. The likelihood (Pkrt) of a remanufacturing process routing being assigned to a certain type of cores is randomly generated to represent the stochastic features [4].(3) The Processing Times. The processing times (ψik) required to perform the necessary remanufacturing operations for a certain part are highly variable. An exponential distribution is used for the processing times in order to simulate the large range of possible values [9].


## 3.2.3. Decision Variables

N
i
t: number of machines in the ith processing cell in period t;
β
i
t: binary variable that has a value of 1 if the machines in the ith processing cell are oriented vertically and 0 otherwise;
L
i
t: length of the ith processing cell in period t;
η
i
t: number of departments occupied by the ith processing cell in period t;
π
i
h
f
t: binary variable that has a value of 1 if the ith processing cell occupies the department in the hth row and ith column as its first location in period t;
(
x
i
t
,
y
i
t
): the centroid of the ith processing cell in period t;
d
i
j
t: distance (in meters) between the two processing cells i and j in period t;
τ
i
t: binary variable that has a value of 1 if the centroid of the processing cell i in period t is different from that in period t-1 and 0 otherwise.As stated earlier, the variability in remanufacturing regarding stochastic returns, probabilistic process routings, and highly variable processing times complicates its facility layout design. Thus the decision variables in the R-DFLP, including the number of machines, the dimensions of individual processing cells, and the detailed location (i.e., the centroid) of each processing cell in a given time period, will be mathematically elaborated in the following.(1) Number of Machines in Processing Cells. The number of machines dedicated to the ith processing cell in period t, Nit, is a function of the number of cores needed to be processed at the time, the average processing times, and machine availability and efficiency:
(1)Nit=∑k∑rQkt·pkrt·ψik·αikrϖi·ci·Dt,
where ci denotes the machine availability (i.e., the available machine hours per workday of a certain machine type) of the ith processing cell (hours/day) and ϖi is the machine efficiency (in percentage) of the ith processing cell corresponding to machine failures.(2) Dimensions of Processing Cells. It is assumed that an individual processing cell should be equipped with a series of machines of the same type to perform a specific remanufacturing process. The dimension of a processing cell as formulated in (2) is pertinent to the numbers of machines dwelt in and their orientations. The number of departments occupied by the ith processing cell is given in (3):
(2)Lit=Mli·1-βit+Mwi·βit·Nit,(3)ηit=ceilLitSl/H,
where Sl/H (i.e., dividing Sl by H) gets the length of each department designed in each row of the floor shop and ceiling (x) returns a value that is the smallest integer not less than x.(3) The Centroids of Processing Cells. For each processing cell, with the required number of departments (ηit) and the first department it occupies (πihft) being determined, its centroids in period t then can be determined by (4). The distance between the centroids of the processing cell i and j in period t is given in (5):
(4)xit=h-0.5·Sl2H+ηit-1·Sl2H,yit=Sw2F+f-1·SwF,xxxxxxxxxxxx∀πihft=1,(5)dijt=xit-xjt+yit-yjt.


## 3.3. Optimization Model

The R-DFLP considered in this paper aims to balance the trade-off between the intercell material handling cost and the increased relocation cost of inefficient layouts.(1) Cost of Intercell Material Handing. While cores go through various processes to be restored back to their useful lives, material handling cost is associated with their logistic activities (i.e., moving from one processing cell to another). The total amount of material handling cost in a given unit of time is an incremental linear function of the total distances of cores moving in the facility, as defined below:
(6)Δt=∑k∑r∑i∑jQkt·pkrt·Γijkr·dijt·θk,
where θk denotes the material handling cost of cores of type k per unit distance.(2) Relocation Cost. Since the best facility layout scheme for one time period may not be efficient for the successive others, the floor shop needs to be reconfigured accordingly by changing the locations of a series of processing cells in a given time period. In this case, machine equipment in these processing cells needs to be removed, transited, and reinstalled which necessarily incur a certain amount of relocation costs, as formulated in the following:(7)Ωt=∑iτit·Nit·δi,∀t>1,τit=0,ifxit=xi(t-1),yit=yi(t-1),1,otherwise,∀t>1,
where τit is a binary variable that has a value of 1 if the centroid of the processing cell i in period t is different from that in period t-1 and 0 otherwise; δi is the relocation cost per machine of the ith processing cell.The R-DFLP problem is then formulated as a nonlinear mixed-integer programming model to minimize the total cost of intercell material handling and machine relocation:(8)Minimizetotal cost=∑t=1TΔt+∑t=2TΩt(9)Subject  to∑t=2T∑iτit·δi≤Bt(10)∑iηit≤F·H,∀t,(11)∑i=3∑j=4Z-1∑k∑rQkt·pkrt·Γijkr=∑i=4Z-1∑j=Z∑k∑rQkt·pkrt·Γijkr,∀t.
Constraint (9) assures that the total relocation cost in any period does not exceed the given budget. Constraint (10) guarantees that the total number of departments occupied by processing cells is not larger than the total available in the shop floor. Constraint (11) represents material flow conservation between disassembly and reassembly.


## 4. The Simulated Annealing Algorithm for R-DFLP

Simulated annealing (SA) is a stochastic gradient method for global optimization that has been successfully applied to a wide variety of sophisticated combinatorial problems [18–20]. While it conducts a local search, it is capable of exploring the solution space stochastically and effectively to prevent being trapped in a local optimum. The SA procedure starts from an initial basic solution s0 and generates a neighbourhood solution s′ by a suitable mechanism. To accept s′ over s0, s′ either is a better neighbourhood solution or has a certain probability to escape from a local minimum though s′ is worse than s0 [21]. The acceptance probability AP is often defined as follows:
(12)AP=1,ifΛ<0,exp⁡-ΛΩ,otherwise,
where Λ=f(s′)-f(s0), Ω is called temperature, a global time-varying parameter, and f is the optimization objective function. At each iteration, the parameter Ω is reduced by a factor ϑ and the chance of choosing an inferior solution decreases as well. The search process, which is generating a neighbour solution, rejecting or accepting it with the above probability, continues until the stopping criterion is met. Given SA is a well-established metaheuristic, the major steps are modified to suit our specific application, followed by the proposed algorithms.


### 4.1. Solution Representation

To address the challenge of R-DFLP with unequal processing cells, this paper proposes a special solution representation that uses two-stage matrices for remanufacturing layout schemes.The first matrixA=[atw]T×W is called sequence matrix that defines the sequence of processing cells to be allocated departments in any given period, where W=Z+F·H-∑ηit. atw is a nonnegative integer, indicating the index of a processing cell that is the wth in sequence to be assigned departments in period t. If atw=0, a department is left unoccupied in the wth action. For each unoccupied department, an intangible machine is assumed to be allocated in it. That is, ∀atw=0, and there exists N[atw]t=1, Ml[atw]=SL/H, and Mw[atw]=SW-(F-1)ε/F. Figure 3(a) gives an example of A that represents allocation sequences in four time periods. For instance, in the first time period, the allocation order is the processing cells 1, 3, 2, and 4, followed by an unoccupied department.(a) Sample matrixA; (b) sample matrix E; (c) sample matrix G derived from A and E.
(a)
(b)
(c)The second matrixE=[etw]T×W has the same dimensions as matrix A that defines the orientation of processing cells in period t. etw takes a binary value of zero or one, indicating the processing cell in the wth allocation sequence is oriented vertically (1) or horizontally (0). In this case, matrix E is determined according to matrix A as defined in (13). It should be noted that if atw>0, then the machine can be located either vertically or horizontally (i.e., etw∈{0,1}). If atw=0, then etw=1 which denotes that the intangible machine can only be oriented vertically in the unoccupied department. Figure 3(b) gives an example of matrix E that aligns to matrix A in Figure 3(b):
(13)etw=β[atw]t∈0,1,atw>0,1,atw=0.To validate the generated solutions, another matrix,G=[gfht]F×H×T, is defined to provide the detailed layout scheme for all processing cells in remanufacturing facility. gfht is the index of the processing cell that occupies the department in the fth row and hth column in period t. If gfht=0, the department in the fth row and the hth column is unoccupied. Matrix G is not part of the solution representation but can be directly derived from matrices A and E. Given the matrices A and E, Lit and ηit can be first calculated, and then matrix G can be determined according to Procedure 1. A sample layout matrix is provided in Figure 3(b) that is directly derived from the sample matrix A in Figure 3(a).Procedure 1:Determination of the layout matrix G.
Input:
A: Sequence matrix
E: Orientation matrix
Output:
Set  h = 1; f = 1;
for  t=1 to T
forw=1 to n
if (atw=0) {
gfht=atw;
if (h≤H)  h=h+1;
else { h=1; f=f+1;}
}
else  {
ξ=η[atw]t;
while (ξ≠0) do {
ξ=ξ-1;
gfht=atw;
if(h≤H)   h=h+1;
else{h=1; f=f+1;}
}
}For each processing cell, with its dimension (i.e.,Lit, ηit) and allocation sequence (i.e., in matrix A) determined, Procedure 1 first selects a corresponding department as its starting location. The selection continues with the right neighbours of the departments from the left to the right until the dimension of the processing cell is met. With the way of allocating departments to processing cells and the assumption that each processing cell only occupies departments in the same row, the layout matrix should meet the following condition:
(14)gfHt≠gf+11t,∀f=1,2,…,F-1.With a feasible matricG obtained, the first department that the ith processing cell occupies can be determined by (15). Then using (4)–(8) the total cost can be calculated accordingly:
(15)πihft=1,∀f=min⁡∀gfht=if,0,otherwise,∀i,t.


### 4.2. Initial Solution Generation

The initial solution is important to the performance of SA as the quality of the final solution relies heavily on that of the initial solution [22, 23]. The procedure of generating an initial solution is given below.


### 4.3. Procedure: SA Initialization

(1)
Generate two sets, the acceptable solution setAS={Av,Ev,v=1,2,…,O} and the invalid solution set IS={(Au,Eu)}. Initialize AS and IS to be empty.
(2)
Randomly generate a solution, 2-tuple{(A*,E*)}. If (A*,E*)∉IS, go to (3); otherwise, restart (2).
(3)
Derive the layout matrixG* based on (A*,E*). If G satisfies constraints (6) and (10), AS=AS∪(A*,E*). If AS=O, go to (4); otherwise, IS=IS∪(A*,E*), go back to (2).
(4)
Calculate the total cost of each feasible solution inAS, and choose the one with the minimal total cost as the initial solution s0 for the SA algorithm (Algorithm 1).Algorithm 1:Proposed SA.
Step  1. Set the annealing parameters: Ω and ϑ and initialize the outer loop counter (ol) and the inner loop counter (il) zero.
olmax⁡, ilmax⁡, and σmax⁡ are maximum outer loop, inner loop, and identical solution, respectively.
Step  2. Apply Procedure SA initialization to generate an initial solution s0. Assign sol=s=s0.
Step  3. While (ol≤olmax⁡ and σ≤σmax⁡) do:
for (il=0; il≤ilmax⁡; il=il+1)
Apply Procedure Neighbour_S(s)
Calculate Λ=f(s′)-f(s);
If (Λ<0) or (exp⁡{-Λ/Ω}≥Random(0,1))
s=s′;
else reject s′;
ol=ol+1;
sol = s;
if (sol = s(ol-1)) σ=σ+1;
elseσ = 0;
Step  4. Output the matrix G based on sol


### 4.4. Neighbourhood Solution Generation

In order to obtain the neighbourhood solution from the current solution, the standard SA neighbourhood structure that features various types of move in work [23] is adopted. The detailed procedure is given below.


### 4.5. Procedure: Neighbour_S(s)

(1)
Randomly select a rowt and a column w in matrix A, and use one type of moves (i.e., m-opt, m∈{1,2,…,n-w}) to exchange atw with at(w+m) and formulate a new matrix A′. Then randomly select a etw>0 in matrix E, and update etw=1-etw to formulate a new matrix E′. A potential neighbour solution s′=(A′,E′) is then formulated.
(2)
Outputs′ as the neighbourhood solution if it satisfies the constraints (6) and (10); otherwise, go back to (1).


## 4.1. Solution Representation

To address the challenge of R-DFLP with unequal processing cells, this paper proposes a special solution representation that uses two-stage matrices for remanufacturing layout schemes.The first matrixA=[atw]T×W is called sequence matrix that defines the sequence of processing cells to be allocated departments in any given period, where W=Z+F·H-∑ηit. atw is a nonnegative integer, indicating the index of a processing cell that is the wth in sequence to be assigned departments in period t. If atw=0, a department is left unoccupied in the wth action. For each unoccupied department, an intangible machine is assumed to be allocated in it. That is, ∀atw=0, and there exists N[atw]t=1, Ml[atw]=SL/H, and Mw[atw]=SW-(F-1)ε/F. Figure 3(a) gives an example of A that represents allocation sequences in four time periods. For instance, in the first time period, the allocation order is the processing cells 1, 3, 2, and 4, followed by an unoccupied department.(a) Sample matrixA; (b) sample matrix E; (c) sample matrix G derived from A and E.
(a)
(b)
(c)The second matrixE=[etw]T×W has the same dimensions as matrix A that defines the orientation of processing cells in period t. etw takes a binary value of zero or one, indicating the processing cell in the wth allocation sequence is oriented vertically (1) or horizontally (0). In this case, matrix E is determined according to matrix A as defined in (13). It should be noted that if atw>0, then the machine can be located either vertically or horizontally (i.e., etw∈{0,1}). If atw=0, then etw=1 which denotes that the intangible machine can only be oriented vertically in the unoccupied department. Figure 3(b) gives an example of matrix E that aligns to matrix A in Figure 3(b):
(13)etw=β[atw]t∈0,1,atw>0,1,atw=0.To validate the generated solutions, another matrix,G=[gfht]F×H×T, is defined to provide the detailed layout scheme for all processing cells in remanufacturing facility. gfht is the index of the processing cell that occupies the department in the fth row and hth column in period t. If gfht=0, the department in the fth row and the hth column is unoccupied. Matrix G is not part of the solution representation but can be directly derived from matrices A and E. Given the matrices A and E, Lit and ηit can be first calculated, and then matrix G can be determined according to Procedure 1. A sample layout matrix is provided in Figure 3(b) that is directly derived from the sample matrix A in Figure 3(a).Procedure 1:Determination of the layout matrix G.
Input:
A: Sequence matrix
E: Orientation matrix
Output:
Set  h = 1; f = 1;
for  t=1 to T
forw=1 to n
if (atw=0) {
gfht=atw;
if (h≤H)  h=h+1;
else { h=1; f=f+1;}
}
else  {
ξ=η[atw]t;
while (ξ≠0) do {
ξ=ξ-1;
gfht=atw;
if(h≤H)   h=h+1;
else{h=1; f=f+1;}
}
}For each processing cell, with its dimension (i.e.,Lit, ηit) and allocation sequence (i.e., in matrix A) determined, Procedure 1 first selects a corresponding department as its starting location. The selection continues with the right neighbours of the departments from the left to the right until the dimension of the processing cell is met. With the way of allocating departments to processing cells and the assumption that each processing cell only occupies departments in the same row, the layout matrix should meet the following condition:
(14)gfHt≠gf+11t,∀f=1,2,…,F-1.With a feasible matricG obtained, the first department that the ith processing cell occupies can be determined by (15). Then using (4)–(8) the total cost can be calculated accordingly:
(15)πihft=1,∀f=min⁡∀gfht=if,0,otherwise,∀i,t.


## 4.2. Initial Solution Generation

The initial solution is important to the performance of SA as the quality of the final solution relies heavily on that of the initial solution [22, 23]. The procedure of generating an initial solution is given below.


## 4.3. Procedure: SA Initialization

(1)
Generate two sets, the acceptable solution setAS={Av,Ev,v=1,2,…,O} and the invalid solution set IS={(Au,Eu)}. Initialize AS and IS to be empty.
(2)
Randomly generate a solution, 2-tuple{(A*,E*)}. If (A*,E*)∉IS, go to (3); otherwise, restart (2).
(3)
Derive the layout matrixG* based on (A*,E*). If G satisfies constraints (6) and (10), AS=AS∪(A*,E*). If AS=O, go to (4); otherwise, IS=IS∪(A*,E*), go back to (2).
(4)
Calculate the total cost of each feasible solution inAS, and choose the one with the minimal total cost as the initial solution s0 for the SA algorithm (Algorithm 1).Algorithm 1:Proposed SA.
Step  1. Set the annealing parameters: Ω and ϑ and initialize the outer loop counter (ol) and the inner loop counter (il) zero.
olmax⁡, ilmax⁡, and σmax⁡ are maximum outer loop, inner loop, and identical solution, respectively.
Step  2. Apply Procedure SA initialization to generate an initial solution s0. Assign sol=s=s0.
Step  3. While (ol≤olmax⁡ and σ≤σmax⁡) do:
for (il=0; il≤ilmax⁡; il=il+1)
Apply Procedure Neighbour_S(s)
Calculate Λ=f(s′)-f(s);
If (Λ<0) or (exp⁡{-Λ/Ω}≥Random(0,1))
s=s′;
else reject s′;
ol=ol+1;
sol = s;
if (sol = s(ol-1)) σ=σ+1;
elseσ = 0;
Step  4. Output the matrix G based on sol


## 4.4. Neighbourhood Solution Generation

In order to obtain the neighbourhood solution from the current solution, the standard SA neighbourhood structure that features various types of move in work [23] is adopted. The detailed procedure is given below.


## 4.5. Procedure: Neighbour_S(s)

(1)
Randomly select a rowt and a column w in matrix A, and use one type of moves (i.e., m-opt, m∈{1,2,…,n-w}) to exchange atw with at(w+m) and formulate a new matrix A′. Then randomly select a etw>0 in matrix E, and update etw=1-etw to formulate a new matrix E′. A potential neighbour solution s′=(A′,E′) is then formulated.
(2)
Outputs′ as the neighbourhood solution if it satisfies the constraints (6) and (10); otherwise, go back to (1).


## 5. An Example

In this section, a sample machine tool remanufacturing facility with 2 rows and 16 processing cells is used to demonstrate the application of the proposed model and algorithm. The workflows of cores in this remanufacturing shop are shown in Figure2. In order to verify the efficiency of the proposed method, two cases are considered in the simulation for a 12-month period. In the static case, the system treats the 12 months as a single period and no reconfiguration is allowed. In the dynamic case, the facility layout is reconfigured at every 4 months based on the changes in the supplier. The proposed SA algorithm is applied to both cases using Matlab 2009. During the implementation of the simulation, the input data is populated as follows. The facility data is shown in Table 1.(1)
The arrival of used machine tools satisfies a Poisson distribution. Thus the arrival of each type of remanufacturable cores in periodt follows a Poisson distribution with the arrival rate, λ∈U[10,15]*π, where π∈{1.0,1.1,1.2,1.3,1.4,1.5} is the arrival rate expansion factor.(2)
The probabilities of each type of cores going through various process routings in a certain period are randomly generatedpkrt∈[0,1] and ∑rpkrt=1.(3)
The average processing time for the cores of typek in the ith processing cell follows an exponential distribution, ψik~Γ(1,λ), where λ∈[0.08,0.10].Table 1
Other data involved in the dynamic facility layout problem.
S
l
=
85 m
S
w
=
20 m
F
=
2
H
=
16
c
i = 8 hours/day
ϖ
I = 95%
B
1
=
B
2
=
B
3
=
B
4
=
1000
θ
1
=
θ
3
=
θ
4
=
θ
5
=
1, θ2=5
δ
i
=
100, ∀i∈{1,2,…,16}
ϖ = 80%
ψ
1
k
=
ψ
2
k
=
15 min
ψ
3
k
=
5 min; ψ16,k=20 min
M
l
1
=
M
l
2
=
M
l
16
=
M
l
14
=
M
l
12
=
M
l
6
=
M
l
8
=
M
l
10
=
6 mMli=3 m, ∀i∈{4,5,6,…}
M
w
1
=
M
w
2
=
M
w
16
=
M
w
14
=
M
w
12
=
3 mMw6=Mw8=Mw10=5 mMwi=2 m, ∀i∈{4,5,6,…}Figure4 compares the total cost of the dynamic layout scheme to that of the static one via twenty-four runs of the SA where the relocation cost δ=100 and π=1.0. It is clear that the dynamic layout scheme outperforms the static one with about 20% cost reduction. Figure 5 shows the optimal layout schemes for both dynamic layout and static one derived from the 8th simulation. Table 2 gives the detailed indication of the letters A–T in Figures 5 and 6.Table 2
The detailed indication of the letters A–T in Figures5 and 6.
A–E
Indication
F–J
Indication
K–O
Indication
P–T
Indication
A
Warehouse for used machine tool
F
Inspection center
K
Chromium electroplating
P
Polishing machine
B
Workshop office
G
Manually scraping center
L
Groove machine
Q
Gear hobbing machine
C
Warehouse for remanufactured machine tool
H
Rough grinding machine
M
Laser cladding
R
Lathe machine
D
Manually disassembly center
I
Cold welding machine
N
Press and straighten
S
Reassembly center
E
Automatic cleaning machine
J
Precision grinding machine
O
Planer machine
T
Unoccupied departmentFigure 4
Performance comparison of the dynamic layout versus the static one (δ=100, π=1.0).The dynamic layout scheme for a whole year (total cost = 83035).
(a)
The dynamic layout scheme for 1–4 months (material handling cost = 23224)
(b)
The dynamic layout scheme for 5–8 months (total cost = 32148)
(c)
The dynamic layout scheme for 9–12 months (total cost = 27663)Figure 6
The static layout scheme for a whole year (total cost = 113250).In order to determine any significance difference between performances of the dynamic layout and the static one, a set of simulations are implemented for different arrival rate level. The average total cost is obtained and compared among the dynamic layout and the static one, as shown in Table3. It is clear to see that the dynamic layout scheme outperforms the static one in terms of generating less cost. A 95% confidence interval for the cost reduction ratio of the dynamic layout is also conducted using the data in Table 3. The approximate 100(1-α)% confidence interval for ϕ is defined as
(16)ϕ^±t1-α/2n-1Sn,orϕ^-t1-α/2n-1Sn≤ϕ≤ϕ^-t1-α/2n-1Sn,
where ϕ^ is the sample mean of ϕ based on a sample of size n; S is the sample standard deviation; t1-α/2(n-1) is the 100(1-α)% percentage point of a t-distributed with n-1 degree of freedom.Table 3
Comparison of the average total cost of two cases.
Arrival rate level
The average total cost
Cost differenceϕ2
Cost reduction ratio(ϕ1-ϕ2)/ϕ1
The static layoutϕ1
The dynamic layoutϕ2
v
=
1.0
100915
78084
28831
0.226239
v
=
1.1
124947
90250
34697
0.277694
v
=
1.2
148068
109396
38672
0.261177
v
=
1.3
178295
135050
43245
0.242547
v
=
1.4
196572
147063
49509
0.251862
v
=
1.5
233200
178090
55110
0.236321
Sample mean
163666.2
122988.8
41677.33
0.249307
Sample standard deviation
48604.4
37507.91
9664.517
0.018437A 95% confidence interval for the cost reduction ratio of the dynamic layout is given by(17)0.249307±2.571×0.01843760.229955≤ϕ1-ϕ2ϕ1≤0.268659.The 95% confidence interval for the cost reduction ratio lies completely above zero, that is,(ϕ1-ϕ2)/ϕ1>0, which provides strong evidence that at 95% confidence there is significant difference between performances of the dynamic layout and the static one.The similar comparison is conducted with the relocation cost changing from 100 to 1900. Figure7 shows that when relocation cost is less than 1200, relocation of machines in dynamic layout would lead to much reduction in material handling and is significant to lower the overall cost. However, when relocation becomes extremely expensive (i.e., higher than 1200), the dynamic layout starts to lose its advantage. The trends in Figure 7 indicate that the dynamic layout is much superior to the static one especially for the facilities with higher flexibility in machine equipment handling. However if the facility has gained much resistance to machine equipment handling, the static layout would outperform the dynamic one with less cost.Figure 7
Performance of the dynamic layout schemes and the static ones with relocation cost ranging from 100 to 1900.


## 6. Conclusion

The remanufacturing process exhibits a high level of uncertainty due to stochastic returns of used products/components and their uncontrollable quality. Such uncertainty challenges the facility layout design for remanufacturing. This paper undertakes this problem and presents an optimization method for remanufacturing dynamic layout design with variable process capacities, unequal processing cells, and intercell material handling. The uncertainties inherent in the remanufacturing facility layout are explicitly analyzed. A dynamic multirow layout model is present for layout optimization with minimization of both material handling cost and machine relocation cost. A modified simulated annealing heuristic is proposed towards the determination of optimal layout schemes with unequal processing cells. The approach is validated through a layout design for a machine tool remanufacturing system and its effectiveness is demonstrated through simulations. Through simulations it is shown that while relocation cost incurs in dynamic layout reconfiguration, the reduction in material handling is significant to lower the overall cost. However, when relocation becomes extremely expensive, the dynamic layout starts to lose its advantage. In the future, factory date will be used to test our methodology through more complicated industrial cases.



---

*Source: 685408-2015-02-23.xml*