Mathematical Problems in Engineering

Volume 2017, Article ID 7249876, 12 pages

https://doi.org/10.1155/2017/7249876

## A New Energy-Aware Flexible Job Shop Scheduling Method Using Modified Biogeography-Based Optimization

School of Information, Zhejiang University of Finance and Economics, Hangzhou 310018, China

Correspondence should be addressed to Wenyu Zhang; gs.ude.utn.e@gnahzyw

Received 22 May 2017; Accepted 19 July 2017; Published 22 August 2017

Academic Editor: Thomas Hanne

Copyright © 2017 Hua Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Industry consumes approximately half of the total worldwide energy usage. With the increasingly rising energy costs in recent years, it is critically important to consider one of the most widely used energies, electricity, during the production planning process. We propose a new mathematical model that can determine efficient scheduling to minimize the makespan and electricity consumption cost (ECC) for the flexible job shop scheduling problem (FJSSP) under a time-of-use (TOU) policy. In addition to the traditional two subtasks in FJSSP, a new subtask called speed selection, which represents the selection of variable operating speeds, is added. Then, a modified biogeography-based optimization (MBBO) algorithm combined with variable neighborhood search (VNS) is proposed to solve the biobjective problem. Experiments are performed to verify the effectiveness of the proposed MBBO algorithm for obtaining an improved scheduling solution compared to the basic biogeography-based optimization (BBO) algorithm, genetic algorithm (GA), and harmony search (HS).

#### 1. Introduction

Under the pressure of sustainable development, manufacturers today must consider not only production efficiency but also energy consumption. In recent years, the demand for energy and the investment in energy have continued to increase. Global energy consumption was 524 quadrillion Btu in 2010 and is expected to increase by 56.5%, to 820 quadrillion Btu, by 2040 [1]. Industry accounted for 52% of energy usage worldwide in 2010 [1]. As one of the most widely used industrial energies, electricity is an important sector that cannot be neglected. However, the rising cost of energy sources to generate electricity such as coal, natural gas, and nuclear energy leads to an increasing electricity consumption cost (ECC). Moreover, the characteristics of electricity make its storage inefficient. Although increasingly heavy investment has been made to support backup infrastructures, with the variable demand of consumers it remains difficult to achieve a trade-off between demand and supply. To address this issue, electricity suppliers have implemented demand response technology.

Time-of-use (TOU) electricity pricing, as the main method of demand response technology, means that the price of electricity is dependent on the electricity consumption at a specific time. That is, it is related to electricity demand. The aim of this method is to encourage consumers and manufacturers to reduce electricity consumption in high-peak periods. For example, the electricity price in high-peak periods can be double that in off-peak periods. Thus, manufacturers can exploit TOU pricing and shift electricity consumption from high-peak periods to off-peak periods. Consequently, ECC, a significant industry cost, can be somewhat reduced.

Three kinds of actions can be implemented to reduce ECC [2]. Manufacturing companies can purchase new energy-efficient machines; this is expensive and has minimal impact. New product forms can be designed. This method is difficult for medium and small-sized enterprises owing to large capital investment. Compared to the aforementioned two methods, production scheduling is more reasonable and applicable because of its minimal investment and strong practicality. The flexible job shop scheduling problem (FJSSP), which has two subtasks, is an extension of the classical job shop scheduling problem. Machine selection is an essential sector in FJSSP. Processing routes are arranged after selection of the appropriate machines. In this study, we assume that each candidate machine has adjustable speeds. Machines in high-speed mode can reduce processing time but increase energy consumption, whereas machines in low-speed mode are relatively energy saving but, more time-consuming. Reasonable arrangements of machines of different speeds in various periods can reduce ECC efficiently under the TOU pricing strategy that contributes to price fluctuations in different periods.

In this study, a new mathematical model is proposed to solve the biobjective problem, namely, makespan and ECC in FJSSP. Moreover, we present a modified biogeography-based optimization (MBBO) algorithm that adopts the extended migration model, combines elitism strategy, and employs variable neighborhood search (VNS). This kind of extension can improve the local search ability and accelerate the convergence of the original biogeography-based optimization (BBO) algorithm. A series of experiments are performed on extended well-known benchmark instances. We compare the MBBO algorithm with the basic BBO algorithm, genetic algorithm (GA), and harmony search (HS) to demonstrate the effectiveness of the proposed algorithm.

The remainder of this paper is organized as follows. In Section 2, we discuss some related works. In Section 3, the problem definition and mathematical model are presented. Section 4 describes the proposed MBBO algorithm for solving the biobjective problem that minimizes makespan and ECC in FJSSP. In Section 5, experiments are performed to validate the effectiveness of the proposed MBBO algorithm. In Section 6, we present conclusions and discuss future works.

#### 2. Related Works

Both industry and academia have addressed green-manufacturing in recent years, and research concentrating on energy-efficient scheduling is gradually increasing. In the single machine environment, Mouzon et al. [3] determined that there existed considerable energy waste, accounting for 80 percent of total energy consumption. Thus, they proposed a “Turn Off/Turn On” method in the idle phase of a single CNC machine to reduce energy consumption. In their follow-up work, they considered total tardiness and energy consumption and used a novel greedy heuristic search method to solve the biobjective optimization problem [4]. Considering variable electricity price, Shrouf et al. [5] utilized GA to reduce ECC in a production-scheduling problem of a single machine. They employed a “Turn On/Turn Off” strategy to determine when to shut down the machine to save electricity cost. Che et al. [6] developed a new mathematical model to address a single machine scheduling problem under TOU electricity tariffs to minimize the total electricity cost and proposed a greedy insertion heuristic to solve it. Ding et al. [7] investigated the unrelated parallel-machine scheduling problem to optimize ECC under the TOU tariff. They considered the TOU setting where electricity prices fluctuate frequently in short periods and proposed a column generation heuristic. Further, a time-interval-based mixed-integer programming formulation was developed allowing a more efficient solution. In the flow shop environment, Bruzzone et al. [8] established a mixed-integer programming model where energy consumption was considered without changing the original jobs’ assignment and sequencing provided by the reference schedule generated by an advanced planning and scheduling system. Lin et al. [9] proposed an integrated model for processing parameter optimization and flow shop scheduling considering makespan and carbon footprint and introduced three carbon-footprint reduction strategies. Lu et al. [10] proposed an energy-efficient permutation flow shop scheduling model considering controllable transportation times solved by backtracking search algorithm and developed a new energy saving strategy. In the job shop environment, Liu et al. [11] used a novel multiobjective GA to minimize the total nonprocessing electricity consumption and total weighted tardiness and employed a “Turn On/Turn Off” method to save electricity.

Diaz et al. [12] concluded that energy consumption is related to the cutting speed and cutting at a higher speed can save more energy than traditional-speed cutting. Fang et al. [13] proposed a new model for solving the flow shop scheduling problem to reduce peak power load, energy consumption, and carbon footprint. They assumed that jobs were operated at dynamic speeds. However, the commercial software they applied directly was only practical for a small-sized problem. Thus, there could exist a gap between the theory and application in industry. To address a parallel-machine scheduling problem considering the total weighted job tardiness and power cost, Fang and Lin [14] attempted to adjust CPU frequencies to force jobs to be processed at a variable machine speed. They assumed that higher machine speed saved time but increased power cost, whereas reduced machine speed sacrificed completion times but led to higher energy cost. They proposed two heuristic algorithms and designed a specific encoding scheme for a particle swarm optimization algorithm for solving the problem. Luo et al. [15] proposed an extended ant colony optimization metaheuristic to solve a hybrid flow shop scheduling problem under the TOU tariff. According to a parameter analysis, the combination of a high-power machine and a low-power machine could save more energy than using two middle-power machines. Sharma et al. [16] proposed a new econological scheduling model under TOU pricing where machine speeds were allowed to vary. Both economic and ecological benefits were simultaneously achieved using a multicriteria metaheuristic optimization method. For the job shop environment, Salido et al. [17] presented a biobjective problem to minimize makespan and energy consumption where machines could work at different speeds.

However, literature concentrating on energy in FJSSP is still lacking. Moon and Park [18] employed constraint programming and mixed-integer programming approaches to minimize production cost in FJSSP considering electricity costs with distributed energy resources. Zhang et al. [19] employed the extended nondominated sorting genetic algorithm II (NSGA-II), considering makespan, machine workload, and carbon footprint, for solving FJSSP. He et al. [20] utilized the nested partitions algorithm to solve a new mathematical model formulated by a mixed-integer programming to reduce energy consumption in FJSSP. Liu and Tiwari [21] developed an optimization approach based on NSGA-II to make the scheduling plans for a carbon fiber reinforced polymer recycling workshop, considering both makespan and energy reduction under the circumstance of flexible job shop. Yang et al. [22] used NSGA-II for solving FJSSP considering makespan and total energy consumption under stochastic processing times. Lei et al. [23] proposed a shuffled frog-leaping algorithm to investigate FJSSP with the objective of minimizing workload balance and total energy consumption. Yin et al. [24] considered productivity, energy efficiency, and noise reduction in FJSSP where machining spindle speed affected the production time, power, and noise.

Metaheuristics and intelligent algorithms have been widely used in computer-integrated manufacturing. In our previous work [25], to determine the manufacturing resource allocation for supply chain deployment, extended GA was utilized for solving the multiobjective decision-making model. In our follow-up work [26], the teaching-learning-based optimization algorithm was improved to plan the distributed manufacturing resource allocation optimally. Pan et al. [27] proposed a chaotic HS to minimize makespan to solve the permutation flow shop scheduling problem considering limited buffers. Ribas et al. [28] proposed an improved artificial bee colony algorithm for solving the blocking flow shop problem by employing various strategies for each phase. Karthikeyan et al. [29] integrated discrete firefly algorithm with local search method to solve FJSSP considering makespan, the workload of the critical machine, and the total workload of all machines. Li and Gao [30] used GA for global search and tabu search for local search to achieve a balance between the intensification and diversification for solving FJSSP. As one of the metaheuristics, the BBO algorithm based on population used in this study, which aims to solve optimization problems, was first proposed by Simon [31] in 2008. Inspired by the biogeography theory, the BBO algorithm assumes that a suitable place for living has a high habitat suitability index (HSI). It is clear that HSI is related to many factors such as rainfall, temperature, and soil. These factors are defined as suitability index variables (SIVs) in the BBO algorithm. To our knowledge, a habitat with a high HSI usually has a relatively significant diversity of species, whereas species of the habitat with low HSI are minimal. Therefore, if the number of species is overly large, then the habitat with a high HSI cannot persist and migration will subsequently occur. In addition to the migration operator, mutation is another main operator that can randomly change SIVs, accordingly affecting the HSI of the specific habitat. A habitat can evolve constantly based on the two operations. The BBO algorithm has been applied and improved in scheduling problems. Rahmati and Zandieh [32] utilized the BBO algorithm to solve FJSSP and performed a comparison between the BBO algorithm and GA for an improved understanding. Wang and Duan [33] improved the standard BBO algorithm by integrating chaos theory and a strategy called “searching around the optimum,” which can prevent the local optima effectively. Lin and Zhang [34] proposed a hybrid BBO algorithm, combined with new heuristics to minimize the makespan in the distributed assembly permutation flow shop scheduling problem. However, the original BBO algorithm has a poor local search ability and converges slowly. The proposed MBBO algorithm in this work adopts an extended migration model and is integrated with VNS, which can enhance the local search ability and accelerate the convergence. This kind of extension avoids the weaknesses of the BBO algorithm effectively.

#### 3. Problem Definition and Modeling

The definition of FJSSP can be described as independent jobs to be processed on a set of machines . A job consists of a sequence of operations . FJSSP aims to achieve the objectives by determining the appropriate machine assignment (MA) and operation sequence (OS).

FJSSP for reducing ECC in this study is an extension of classical FJSSP. It not only considers the traditional objective makespan but also ECC under the TOU tariffs. We assume that the machine frequency is adjustable and jobs can be operated at variable speeds. Thus, different from classical FJSSP with two subproblems containing MA and OS, another subproblem that considers the selection of speeds is included.

The objective of this study is to achieve a balance between makespan and ECC by assigning jobs to the appropriate machines with selected operating speeds. To make this scheduling problem more concise, there are some assumptions to be satisfied.(1)An operation, once started, cannot be interrupted until it is completed.(2)Each machine can process only one job at a time and each job can only be processed on one machine at a time.(3)The operations of each job have priority constraints.(4)All machines are available during the planning horizon and all jobs can be processed at time zero.(5)The setup and adjustment periods of the machines are negligible.

##### 3.1. Notations

To understand the mathematical model more clearly, the corresponding parameters and decision variables are explained as follows:

*Sets* : set of jobs : set of machines .

*Parameters* : index of jobs, : index of machines, : index of operation sequences, : index of time periods, : index of machine processing speeds, : the* j*th operation of : number of time periods of processing on machine at speed : start time of : finishing time of : makespan of the schedule : processing energy consumption on machine at speed at period : idle energy consumption on machine at period : the electricity cost at period ECC: total electricity cost of the schedule

*Decision Variables* : binary variable, if is processed on machine at speed , then ; otherwise, : binary variable, if is processed on machine at speed at period , then ; otherwise, : binary variable, if machine is at the processing status at speed during period , then ; otherwise : binary variable, if machine is at the standby status during period , then ; otherwise .

##### 3.2. Formulations

During the production process, the time and energy consumption of the setting up and transforming periods are relatively small compared to those of the processing and standby periods. In this study, we only consider three machine statuses: processing, standby, and off status.

From the following two equations, we can determine the makespan and total electricity cost:

Before we transform the biobjective problem, which aims to minimize the makespan and ECC into a monoobjective, the normalization method provided in [35] must be performed. Because the values of the makespan and ECC belong to intervals with different lengths, a domination could occur leading to an inefficiency of the evaluation. Thus, the normalized makespan value and normalized total electricity cost value can be obtained. By setting weights and manually in the interval between zero and one, a balance can be achieved between time and ECC. Weight coefficients summed to one are used to reflect the respective importance of the objectives. For example, if the decision maker considers time more important than the other objectives, the weight coefficient of time can be set higher manually. This kind of method is more flexible for different decision makers when planning production schedules.

The mathematical model is indicated as follows.

*Objective*

*Constraints*

Equation (2) is the objective function, which is aimed at minimizing the related makespan and ECC. Constraint (3) guarantees that each operation can be completed without interruption once it starts. Constraint (4) guarantees that only when the prior operation is finished can the next operation of the same job be processed. That is, it is the priority constraint of the operations. Constraint (5) ensures that each operation cannot be processed on more than one machine at one speed at the same time. Constraint (6) ensures that each machine can perform at most one operation at any time. Constraint (7) describes the relationship between decision variables and . Constraint (8) confines each machine to be operating at only one specific speed at each time period and requires that each machine must be in one of the three possible states at any time. Constraint (9) imposes the nonnegative restriction on the start time and finishing time of each operation. Constraints (10)–(13) specify that the decision variables , , , and are binary.

#### 4. MBBO Algorithm for FJSSP

In this section, the MBBO algorithm is proposed to solve the FJSSP considering both makespan and ECC. To enable the proposed approach adaptive to the investigated problem, a new three-vector string representation is employed. In addition, VNS is integrated to improve the performance of the basic BBO algorithm.

##### 4.1. Encoding and Decoding

The first step of solving FJSSP using the MBBO algorithm is to obtain a reasonable encoding and decoding scheme that can effectively represent the problem. Owing to the quality of the traditional FJSSP, the encoding scheme design includes two subtasks, namely, MA, which selects the machine where the operation will be processed, and OS, which determines the specific operation sequence of all the jobs. In this study, to meet the assumption that machine frequency is adjustable, a new subtask called speed selection (SS) is introduced. It must be noted that the item called habitat in the BBO algorithm is similar to that called chromosome in GA. Accordingly, we extended the representation employed in [32] to form a three-vector string. For a clear and intuitive understanding, we present Table 1 as a simple example to describe the encoding and decoding scheme.