Abstract
A parallel machine scheduling problem in plastic production is studied in this paper. In this problem, the processing time and arrival time are uncertain but lie in their respective intervals. In addition, each job must be processed together with a mold while jobs which belong to one family can share the same mold. Therefore, time changing mold is required for two consecutive jobs that belong to different families, which is known as sequence-dependent setup time. This paper aims to identify a robust schedule by min–max regret criterion. It is proved that the scenario incurring maximal regret for each feasible solution lies in finite extreme scenarios. A mixed integer linear programming formulation and an exact algorithm are proposed to solve the problem. Moreover, a modified artificial bee colony algorithm is developed to solve large-scale problems. The performance of the presented algorithm is evaluated through extensive computational experiments and the results show that the proposed algorithm surpasses the exact method in terms of objective value and computational time.
1. Introduction
Parallel machine systems are widely adopted in a variety of manufacturing environments, such as the semiconductor manufacturing industry [1] and the electronics industry [2]. Minimizing makespan is one of the commonly used objectives in manufacturing scheduling problems [3, 4]. Parallel machine scheduling problems have both theoretical and practical importance. Many literatures assume that the parameters of the problem (i.e., the processing time and job release time of jobs) are known in advance precisely before production process begins. However, the results of the solution derived under the deterministic assumption may deviate considerably in real situations [5]. In practice, it is difficult to get the exact parameters before production process begins due to some uncertainties, such as machine conditions, production environments, and jobs’ characteristics [6–9]. To characterize and overcome the impact of uncertainties, several robust scheduling approaches are proposed to enhance the quality and stability of the derived solution in real situations.
Unlike the typical parallel machine system, a parallel machine system for plastic product involves mold changes during the process of production, which is also known as setup requirements. In a plastic production system, jobs are allocated to one of parallel injection machines and the corresponding injection mold is required to be installed onto the injection machine before the injection process. Jobs belonging to one family can be processed with the same mold and the mold changing process costs a period of time. The production planner would try to arrange the jobs such that jobs that belong to the same family are processed together in order to avoid extra mold change process and improve efficiency as well. For two jobs that are scheduled to be processed consecutively but belong to different families, mold change must be conducted before the next job can be performed. Therefore, reducing mold changing time is crucial in plastic production scheduling.
This paper investigates the uniform parallel machine problem in the plastic production system, which involves uncertain processing time, job release time, and setup time for mold. It is an extended problem studied by [5]. The preliminary objective is to minimize the makespan. A scenario is firstly defined as a possible realization of processing times and job release time for all jobs. Due to the uncertainty in job processing time and release time, the makespan of a given solution can be different under various scenarios. It is proposed to identify an optimal solution with the strongest stability across all scenarios. The performance (or stability) of a given solution under its worst-case scenario is the major concern in robust scheduling. For generality, it is assumed that the processing times and job release time lie in their respective intervals with lower bound and upper bound value. The intervals of the corresponding parameters are considered to be known but not fall into any statistical distribution due to asymmetric information [5, 9]. And then the robust deviation criterion is adopted to evaluate the robustness of each candidate schedule. The robust deviation criterion is also known as min–max regret criterion, which aims to find out a solution with a minimum maximal deviation across all possible scenarios. This criterion has its advantages in the highly competitive environment for the reason that the robust decision performs well in any set of potential realizable scenarios [5]. The robust deviation criterion confines the magnitude of missed opportunities by identifying a schedule that possesses a performance close to that of the optimal (or near optimal) decision in any scenarios.
To the best knowledge of the authors, parallel machine scheduling problem with uncertain processing time, ready time, and mold change consideration has not been covered by the other researchers in the parallel machine scheduling problem. In this paper, a mixed integer linear programming (MILP) formulation is introduced to identify a robust schedule with minimum maximal regret across all scenarios. The concept of critical machine proposed in [5] is revised and adopted in the problem to eliminate the worst-case scenario into a finite number of extreme point scenarios in order to evaluate the robust deviation of given solution by calculating the maximal regret. Exact algorithm based iterative relaxation procedures are presented, and a modified artificial bee colony algorithm is proposed for the research problem. To demonstrate the effectiveness and efficiency of the proposed heuristic algorithm, a set of testing problem is carried out.
The contribution of this paper is shown as follows:(1)A parallel machine scheduling problem with consideration of uncertain processing time, job release time, and mold changing time is studied. For generality, the uncertain data is assumed to lie in intervals, which capture the situation in real environment.(2)Swarm intelligence algorithm is proposed to solve the NP-hard problem. The computational results demonstrate the stability and effectiveness of the proposed algorithm. To the best of our knowledge, this research is the first to adopt artificial bee colony algorithm in robust optimization for a parallel machine scheduling problem.The rest of this paper is organized as follows. Section 2 provides the literature review of works related to parallel machine as well as robust scheduling approaches. Section 3 presents the description and the mathematical model of the proposed problem. The exact algorithm based on iterative relaxation is provided in Section 4. The modified artificial bee colony algorithm for large-scale problem is proposed in Section 5. Results of extensive computational experiment and comparison are provided in Section 6. Section 7 provides the conclusion and direction for further works in related fields.
2. Literature Review
Stochastic approaches for tackling scheduling problem under uncertainties are available [10], but some of these approaches have their limitation due to the strict prerequisite and assumption [5]. For example, stochastic approaches require certain information on probability distribution of processing time or release time of each job, which can be inferred on the condition that a substantial amount of historical data is available [9, 11]. However, such amount of historical data is unavailable in highly uncertain environment and the only information is an educated guess of the lower bound and upper bound of some parameters, such as processing time and ready time [12]. For some one-time jobs, decision-makers are more interested in obtaining a robust schedule, which is against the worst-case performance across all scenarios, rather than obtaining an expected optimal performance under an expected situation. Robust deviation approach, which is known as min–max regret [13], is suitable for these circumstances to obtain such solution [14, 15]. The robust deviation approach is widely adopted in various combinational optimization problems when the input data are presented as intervals, such as the shortest path [16, 17], spanning tree [14, 17–19], and production problems [5, 20]. This approach yields a satisfactory result in an uncertain environment. Moreover, robust scheduling approaches are also adopted in other problems, such as maritime transportation problem [21, 22], routing problem [23], and scheduling problem in public health service department [24]. Most recently, a min–max regret makespan minimization in an identical parallel machine scheduling environment with interval data is studied [5]; in particular, they considered the processing time of jobs lies in respective intervals. To solve this nondeterministic polynomial-time hard (NP-hard) problem, the concept of critical machine and extreme point scenarios and two properties are proposed to avoid visiting an infinite number of possible scenarios and eliminating worst-case scenario into finite number of extreme point scenarios. A makespan minimization problem with interval job processing time on identical parallel machines is addressed in [5] for the first time in related field. Furthermore, the approach of robust deviation on a uniform parallel machine scheduling problem is adopted in [12] to minimize the total flow time with uncertain processing time which lies in respective intervals. The concept of worst-case scenario used to identify the maximal regret for a feasible solution is adopted in [5, 12] and then the exact and heuristic algorithms are proposed to find the robust schedule across all scenarios.
Scheduling problem with sequence-dependent setup time is a very active research area [25–28]. In the plastic manufacturing process, injection operation is a typical single-stage manufacturing process, which requires the plastic to be injected into specified mold to produce plastic products in different shape. Several jobs are manufactured by a single machine with their respective molds [29]. A setup of mold is required for a product to be produced in injection operation, while setup is a sequence-depending operation and parameters such as time for mold change and install time vary depending on two consecutive jobs [30]. If two consecutive jobs are processed by different mold, cost of setup is induced during the manual mold change operation. Therefore, the time for mold change as well as reinstall cannot be ignored in a mass production situation and planner would try to arrange those jobs by the same mold to be processed together in order to reduce the number of mold changes during production. The sequence-dependent family setup time represented by interval data on a single machine scheduling problem to minimize the total flow time is considered in [20].
Exact and heuristic algorithms have been proposed to solve min–max regret problem. A min–max regret problem of minimizing the total flow time on a single machine in [6] and then this problem are investigated further by the other researchers. In order to reduce the computational effort, different approaches have been proposed for the min–max regret model. It is proved that the optimal schedule under the mid-point scenario guarantees a 2-approximation of the optimal solution [15]. Heuristics method is an approximation algorithm to obtain near optimal result for min–max regret model. Typical example is shown in [5]. Some heuristic algorithms based on job swap moves and insert moves are presented by researchers. A hybrid tabu search algorithm for batching and sequencing decision-making in a single machine scheduling problem is proposed by Suppiah and Omar [31]. Job swapping and insertion approaches are applied to generate neighborhood solution and employed arcs in which the solution appears in the form of arcs in the tabu list to represent the sequence of job on a single machine; in addition, they also implemented a search depth strategy in the process of neighborhood generation to eliminate noneffective moves so as to reduce the computational burden while obtaining final schedule with outperformed quality. The work of Bilge et al. [32] points out that, for the situation where neighborhood of solution is in large number or its elements are expensive to evaluate, it is essential to restrict the number of solutions examined on a given iteration to screen the neighborhood so as to concentrate on promising moves at each iteration. Three candidate list strategies are proposed to confine the number of neighborhood solution for calculation efficiency.
Various studies on parallel machine problem have, respectively, considered uncertain processing time, arbitrary ready time, and mold changing time in the model development stage while there are no studies concerning these factors altogether in an integrated robust scheduling model. Therefore, this paper investigates these factors altogether and formulates a novel model in a plastic production environment. The presented paper is an extended version of the previous paper, containing more uncertain factors and features in a practical manufacturing system.
3. Problem Formulation and Mathematical Model
It is started by describing the problem and the definition of the maximal regret of a feasible schedule. Then the MILP is presented.
3.1. Problem Description
The problem under consideration deals with the scheduling of parallel machine over an assigned planning horizon in order to minimize the makespan. The job is processed by machine and mold and each job can be processed only once. Under the previous assumptions, the problem can be modelled as a mixed integer linear programming formulation as stated below.
3.2. MILP Formulation
Notations : job, , where is the total number of jobs, : machine, , where is the total number of machines, : mold installation time for job : mold removal time for job : job ’s processing time under scenario , : job ’s arrival time under scenario , : processing speed of machine : scenario which is a possible realization of the processing times and arrival time of jobs, : a large positive number
Decision Variables : 1 if job is scheduled to be processed on machine and 0 otherwise : 1 if job is scheduled to be processed after job on the same machine and 0 otherwise : the beginning time of job to be processed on machine in schedule under scenario : the completion time of job scheduled on machine in schedule under scenario : the completion time of machine in schedule under scenario A feasible solution should satisfy and , , . That is, each job should be processed exactly once in parallel machine system and there exists a processing sequence for jobs to be processed on the same machine. Let be the set of feasible schedules.
Completion time of job in schedule under scenario can be defined as follows:Calculation of completion time of machine under scenario :
The completion time of machine in schedule under scenario is equal to the completion time of the job that is scheduled to be processed by the last machineThat is,The starting of job on machine is equal to or greater than the arrival time of job under scenario :In addition, if job and its preceding job are processed by different mold, the time for mold change is involved before job starts to be processed:The makespan of schedule under scenario is equal to the completion time of the machine that finished last in parallel machine system:orWe define the minimum makespan under scenario as and the corresponding optimal schedule as :For a given schedule , its regret under scenario is defined asThe maximum regret of schedule isThe robust parallel machine scheduling with interval processing time and arrival time and mold change can be formulated asThe proposed problem is a generalization of classical schedule problem and the optimal solution derived from (11) corresponds to the robust schedule across all possible scenarios.
The parallel machine robust scheduling with uncertain processing time and arrival time (RS) can be formulated as follows:We note that the above formulation is nonlinear due to the two operators in the objective function. However, the original model of RS can be transformed into a mixed integer linear programming model:After linearization of the original model of RS problem, the formulation still cannot be solved directly due to an infinite number of possible scenarios. To tackle this problem, the primal goal is to identify the maximal regret of each feasible solution and then select the solution with the minimum maximal regret across all solutions. Therefore, two properties are proposed to identify and confine worst-case scenario to a finite number of extreme point scenarios, which means that we do not have to visit all possible scenarios to obtain maxima regret of a given solution. In order to capture the uncertainties in real manufacturing environment, we assume that the processing time and arrival time of each job are uncertain by lying in their respective intervals. Therefore, the number of possible scenarios is infinite under such assumption and it is impractical to enumerate all possible scenarios to get the worst-case scenario for a given solution. The derived properties are used to narrow down the worst-case scenario to several extreme point scenarios.
Definition 1. A machine is said to be critical in a schedule under scenario , if is the machine with maximum completion time in under :
Definition 2. An extreme point scenario for schedule is defined as follows:
Property 1. For any schedule , let be a worst-case scenario for in which machine is critical. Then a scenario exists for schedule such that(a)machine is also critical in under ;(b)scenario is a worst-case scenario for .Assume that a given solution has a worst-case scenario that satisfies conditions (a) and (b) of Property 1. Since machine is critical in under scenario , the makespan and maximal regret for solution can be calculated as follows.
Note that the set containing all the jobs on machine in schedule is denoted as , . If , then .
Proof. Worst-case scenario can be transformed into scenario by decreasing to and to for all and by increasing to and to for all , and is also the critical machine in under . Let denote the increase in the makespan under . Let be the optimal schedule for scenario . Because the makespan in cannot increase by more than under scenario in comparison with , we getFrom the abovementioned inequality, we have , indicating that the maximal regret cannot decrease if we replace with , which is also the worst-case scenario for .
Property 2. The maximal regret of a given solution can be expressed as follows:
Proof. Equation in Property 2 shows thatSuppose by contradiction that there exists a machine such thatSince , where ,
4. Exact Algorithm
RS is a min–max problem that can be solved using a general iterative relaxation (IR) procedure proposed by [33–35]. First, the set of all possible scenarios is replaced in RS with a finite set of scenarios , resulting in relaxed mixed integer program (RS-relaxed): Note. In this iterative relaxation approach, we replace the set of all possible scenarios () with , and the scenario in is the worst-case scenario for the feasible solution that is identified in the last iteration.
The minimal makespan can be solved by solving a mixed integer programming if the scenario is fixed.
We refer to the constraints and associated with as regret cuts.
Iterative Relaxation Algorithm. The IR procedures stop when , where is the maximum regret for current relaxed solution .
Step 0. Set ( is lower bound for regret value) and ( is upper bound for regret value); choose an initial solution .
Step 1. Identify solution ’s worst-case scenario and its respective through Property 1 and Property 2. If , then , . If , go to Step .
Step 2. Add regret cuts and ; to RS-relaxed scenario .
Step 3. and are identified by solving RS-relaxed. Set ; go to Step .
Step 4. Stop.
5. Swarm Intelligence Approach
5.1. Modified Artificial Bee Colony Algorithm
5.1.1. Description of Modified Artificial Bee Colony Algorithm
ABC algorithm has been well-studied in parallel machine scheduling [36–38]. ABC algorithm is a promising swarm intelligence algorithm for optimization problems, which simulates the foraging behavior of honeybees. Three types of bees serve different functions to contribute to the forage acquisition for the bee social structure in the hive. The collected nectar will be stored as a food supply during the shortage of nectar. Employed bees and onlooker bees respond to obtain nectar from the nearby plants, while the main force of scout bee is to gather nectar information to guide the foraging behavior of the employed bees and onlooker bees. The design of the ABC algorithm follows a recursive search procedure with four major phases, including initialization phase, employed bee phase, onlooker bee phase, and scout bee phase. The process flowchart is shown in Figure 1.
The size of the bee colony is equal to , which includes employed bees, onlooker bees, and one scout bee. Approximately half of honeybees are as employed bees, and the remaining are onlooker bees. The number of solutions in ABC algorithm is equal to , where . Each solution represents an occupied food source by one employed bee. Hence, the population of the solution is represented by . The role of employed bee is to gather nectar back to the hive; then the information of the food source will be shared among the hives. Onlooker bees will select one food source to further exploit based on the abundance level of the known nectar. Scout bee will seek for a new solution once a food source is exhausted.
5.1.2. Initialization Phase
For the initialization phase, the conventional approach is to generate a random assignment and sequence for an optimization problem. This prolongs the convergence to the global optimum and fails to achieve the global optimum, as the process of converging to the promising regions involves certain unsuccessful neighborhood searching under the large size solution space. A solution will be abandoned when there is no improvement after certain unsuccessful search. Therefore, a constructive heuristic for the initialization phase insists on developing the best possible solution at the initial stage in order to reduce the searching effort of the ABC algorithm in converging to the global optimum. The composition of the minimal makespan in PMS model aims to obtain a solution with the least total length of the parallel machine schedule to complete all the available jobs. The solution quality of an initial solution by sorting the ready time of each job in ascending order is expected to be the best alternative in timely fashion compared to random assignment. While maintaining a certain quality of initial solution, the diversity between solutions should be considered. Otherwise, certain solution space could not be visited. At the same time, it is difficult for the population-based metaheuristics to escape from the premature convergence [39]. One way of reducing the possibility of converging to an identical but not a global optimal solution for the population-based metaheuristics relies on randomness in the construction of initial solution. In our preliminary study, it is effective to develop a fair quality of solution by random machine assignment for a sorted sequence.
5.1.3. Employed Bee Phase
Each employed bee performs a neighborhood search on a solution to generate a neighborhood solution by three neighborhood operators to enhance the solution quality in terms of its objective value. A greedy selection is applied herein to continuously improve a solution by comparing the objective values between and . If the solution quality of neighborhood solution is better than the original one, the neighborhood solution will replace the previous solution . The neighborhood operators include swap operator, insert operator, and crossover operator as shown in Figures 2, 3, and 4, correspondingly. The neighborhood search is not restricted on either the same or another runway by randomly selecting two or more elements. The accumulative trial of unsuccessful improvement on a solution will be increased by 1 if the objective value of a neighborhood solution is worse than the previous solution. The crossover operator will select a single runway solution with the longest completion time and choose another solution by fitness selective probability to combine and develop an offspring solution.
5.1.4. Onlooker Bee Phase
Each onlooker bee will further act as neighborhood operators to a specified solution based on a probabilistic selection with regard to the fitness values. The fitness value of a solution is measured by (33). A larger number of fitness values imply a better quality of a solution . After the computation of fitness value for all the solutions, (34) is able to obtain the individual selective probability of each solution . A cumulative probability is computed by aggregating the preceding individual probabilities. The selection criterion follows the approach of roulette wheel selection. A random number will be generated for each onlooker bee for selecting a solution to further search by neighborhood operators:
5.1.5. Scout Bee Phase
During the neighborhood search, the cumulative number of unsuccessful update will be recorded as . A large number of imply a possibility that solution is tripped in a local optimum. The employment of scout bee phase is to abandon a solution once it reaches the maximum tolerance of searching. If one solution cannot be further improved by several iterations and reaching the maximum tolerance , the scout bee will be replaced by an exhausted solution by an initial solution.
5.1.6. Min–Max Regret Approach Using Artificial Bee Colony Algorithm
The design of the min–max regret approach is to obtain a robust PMS schedule to minimize the maximum regret value by considering a number of scenarios. Two main stages require optimization techniques to resolve the min–max regret model, which are optimal makespan under the worst-case scenario and makespan under robust schedule. The process workflows using the ABC algorithm in deriving the worst-case optimal makespan and robust makespan are the same, except the computation of objective value and fitness value. The objective function in optimal makespan under the worst-case scenario is , while the objective function in robust ALP makespan by considering all the worst-case scenarios follows (28) and (29).
6. Numerical Experiments and Results
In this section, the performance of the modified ABC algorithm is evaluated by random robust PMS instances and compared with the results using mixed integer programming. The MIP is computed by IBM ILOG CPLEX Optimization Studio 12.6.3 for comparison. All algorithms are written in C# language with visual studio 2015 on a computer with Intel Core i7 3.60 GHz CPU and 16.0 GB ram under the Window 7 Enterprise 64-bit operating environment. In order to evaluate the algorithm performance, two approaches have been considered in the numerical experiments, as shown in Table 1. Approach 1 using mixed integer programming in resolving RS and is considered as a baseline for comparison. In approach 2, the optimization technique for RS and is the proposed artificial bee colony algorithm. In this regard, the performance of convergence and computational speed of the proposed artificial bee colony algorithm can be measured by comparing the values obtained in exact method.
The test running is conducted in preliminary study to obtain the parameter setting. The maximum iteration MaxIter is equal to . The maximum tolerance of unsuccessful update is equal to . The colony size is set to be 80, and the number of solutions is 40. These results are also evaluated together with the optimal solution from MIP to measure the deviation between exact algorithm and swarm intelligence. The maximum computational time for MIP and proposed ABC algorithm is 7200 seconds.
Both algorithms are evaluated by test instances with jobs and machines. For each job, the lower bound and upper bound of ready time and processing time fall into a random uniform distribution integer with intervals of , , , and correspondingly, where value is considered to be (0.2, 0.4, 0.6, 0.8, 1.0). The total number of test instances is 45 with any combination of , and . The processing speed is equal to a random uniform distribution from 1 to 10, where . Assuming there are only three types of injection molds, the mold installation time and mold removal time of injection molds for each group are defined as (30, 60, 90) and (10, 20, 30). The number of jobs for each group is equally distributed.
The regret upper bound value and regret lower bound value have been recorded using MIP to indicate the performance of the iterative relaxation algorithm. If the regret gap is equal to 0%, the recorded robust optimal value implies a true optimum. Otherwise, the model has not been able to determine the optimal solution within 7200 seconds if the regret gap is larger than 0%. In order to compare between the performance of MIP algorithm and the performance of MABC algorithm, the deviation from the best-known solution indicates the deviation between average optimal and the optimal value using Cplex with limited computational time.
Table 2 shows the average computational performance for robust parallel machine scheduling using modified ABC algorithm in 10 run times. The results indicated that the computational time is exponentially increased when the problem size grows for both algorithms. For the problem size of , the robust optimal solution obtained by MIP performs better than the result using MABC algorithm in terms of optimal value and computational time. This implies the certain computational overhead for MABC algorithm in small size instance. When the problem size increases, MIP was not able to determine the optimal solution within a reasonable time, while MABC algorithm was able to obtain near optimal value. In order to measure the variability of algorithm performance, the number of best solutions denotes the number of optimal values (maximum number of 10) that surpasses the result obtained by MIP in Table 3. These results can provide evidence that the variability of MABC algorithm performance for the test instance with is in low level with respect to the quality of the solutions. The number of better solutions #best represents the number of solutions found by modified ABC algorithm in 10 run times which outperforms the solution obtained by MIP with time limit. As for the large size instance, the modified ABC algorithm is able to obtain better solution with lower computational effort than the exact method.
7. Concluding Remarks
Parallel machine model has been widely adopted in the manufacturing environment and its scheduling problem is extensively studied in literatures and minimizing makespan is the common objective in parallel machine scheduling problem. However, most of the literatures have the assumption that the parameters of job such as processing time and ready time are precisely known, which is counter to the practical situation for the set of jobs that are processed for the first time. Uncertain processing time and ready time in parallel machine scheduling are the common phenomenon in the production scheduling. In practice, the ready time of a job is uncertain and cannot be estimated by an exact probability distribution due to the effect of external factors, such as delivery variability, in-bound logistics delays, labour issues, and transportation schedules. The processing time of a job may deviate from the standard operation time, which is caused by corrupting operational events, corrective action, and reworking.
In order to characterize the uncertainties of parameters, it is assumed that the exact value of processing time and ready time of each job is unknown and the available information is the bound of each parameter. To resolve this robust version of problem, the min–max regret approach is introduced herein to obtain a robust parallel machine schedule that considers all possible scenarios. We first introduce a mixed integer linear programming formulation to characterize the problem and then propose two properties to eliminate worst-case scenario The proposed two properties allow the model obtain the maximal regret for each feasible schedule with a finite number of extreme point scenarios.
An exact algorithm based on general iterative relaxation algorithm is developed to solve the proposed problem. However, it is difficult to obtain a solution from large size instances by the exact algorithm since identifying worst-case scenario for each feasible is computationally expensive. Therefore, a modified artificial bee colony algorithm is proposed to solve the large-scale parallel machine scheduling with uncertainty. The experimental result indicates that a robust near optimal solution can be obtained with a reasonable time. The strength of artificial bee colony algorithm follows the three main features to obtain a near optimal solution by utilizing decentralization of honey bee, self-organizing in swarm performance, and collective behavior in algorithm structure. The proposed artificial bee colony algorithm is widely adopted in a static parallel machine scheduling and other related problems. In order to enhance the exploitation and exploration of the proposed algorithm to achieve a high level of optimization technique for robust modelling, certain modification on artificial bee colony algorithm has been proposed.
Future work can be concluded as follows: extension of parallel machine scheduling with realistic uncertainties, constraints, and objectives to retain feasible solution in actual operation and optimizing the algorithm structure and iteration process with other known optimization techniques to achieve optimum with less computational effort.
Notation of Artificial Bee Colony Algorithm
| : | The size of bee colony |
| : | The number of colony solutions |
| : | The maximum number of iterations |
| : | The dimension of an independent solution |
| The position of each solution in bee colony | |
| : | The objective value of solution |
| : | The fitness value of solution |
| : | The probability of an individual solution among the entire colony in terms of fitness value |
| : | The cumulative probability of an individual solution in ascending order among the entire colony in terms of fitness value |
| : | The neighbor solution of an individual solution |
| : | The accumulated trial value of an individual solution , which cannot enhance the quality of solution in terms of its objective value |
| : | The maximum tolerance |
| : | Random number, . |
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research is supported by The National Natural Science Foundation of China . 71201099, Innovation Program of Shanghai Municipal Education Commission . 14YZ111, Shanghai Young Eastern Scholar Programme QD2015041, Shanghai Pu Jiang Program (no. 13PJC066), Shanghai Youth Teacher Foundation (no. ZZshhs13021), and the Hong Kong Polytechnic University. The authors’ gratitude is also extended to the research committee and the Department of Industrial and Systems Engineering of the Hong Kong Polytechnic University for support of this project (RUF1 and RU8H).