Scientific Programming

Scientific Programming / 2020 / Article
Special Issue

Optimization Models and Algorithms for Services and Operations Management 2020

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 8860384 | https://doi.org/10.1155/2020/8860384

Huu Dang Quoc, Loc Nguyen The, Cuong Nguyen Doan, Toan Phan Thanh, Neal N. Xiong, "Intelligent Differential Evolution Scheme for Network Resources in IoT", Scientific Programming, vol. 2020, Article ID 8860384, 12 pages, 2020. https://doi.org/10.1155/2020/8860384

Intelligent Differential Evolution Scheme for Network Resources in IoT

Academic Editor: Tingsong Wang
Received05 Jun 2020
Accepted19 Jun 2020
Published14 Jul 2020

Abstract

Scheduling is a fundamental factor in managing the network resource of the Internet of things. For IoT systems such as production lines, to operate effectively, it is necessary to find an intelligent management scheme, i.e., schedule, for network resources. In this study, we focus on multiskill resource-constrained project scheduling problem (MS-RCPSP), a combinational optimization problem that has received extensive attention from the research community due to its advantages in network resource management. In recent years, many approaches have been utilized for solving this problem such as genetic algorithm and ant colony optimization. Although these approaches introduce various optimization techniques, premature convergence issue also occurs. Moreover, previous studies have only been verified on simulation data but not on real data of factories. This paper formulated the MS-RCPSP and then proposed a novel algorithm called DEM to solve it. The proposed algorithm was developed based on the differential evolution metaheuristic. Besides, we build the reallocate function to improve the solution quality so that the proposed algorithm converges rapidly to global extremum while also avoiding getting trapped in a local extremum. To evaluate the performance of the proposed algorithm, we conduct the experiment on iMOPSE, the simulated dataset used by previous authors in their research studies. In addition, DEM was verified on real dataset supported by a famous textile industry factory. Experimental results on both simulated data and real data show that the proposed algorithm not only finds a better schedule compared with related works but also can reduce the processing time of the production line currently used at the textile industry factory.

1. Introduction

Internet of things systems include not only sensors and automated machines like robots but also intelligent control and management algorithms [1]. Intelligent algorithms for managing network resource play a key role in controlling the operation of automated machines such as robots, enabling them to devote their full capacity to serve users.

Multiskill resource-constrained project scheduling problem (MS-RCPSP), a combinational optimization problem, has received extensive attention from the research community in recent years due to its advantages in network resources management and scheduling for Internet of things system. The main purpose of the MS-RCPSP is to find the way for assigning tasks to resources so that the execution time is minimal.

In recent years, MS-RCPSP is widely used, and this problem is found not only in critical civilian fields but also in the military setting. MS-RCPSP appears on a broad scale of military operations such as planning an onboard mission or routing an unmanned aerial vehicle (UAV) [2, 3].

The rest of the paper is organized as follows. In Section 2, we present some related works such as previous approaches to solving the MS-RCPSP. Specifically, the summary of the differential evolution method is included in this section. Differential evolution is a well-known evolutionary algorithm that has been applied effectively to solve NP-hard problems. In this paper, the proposed algorithm is inspired by the differential evolution method. In Section 3, we present the system model and the statement of MS-RCPSP. The formulation of MS-RCPSP is introduced in this section. Section 4 describes the proposed algorithm, which is built on differential evolution strategy [4]. In this section, we focus on presenting the reallocate function, which is the key component that creates the strength of the proposed algorithm. Besides, this section introduces our design for schedule representation and a new measurement model for measuring the difference between two schedules. The proposed algorithm (DEM) is also introduced in this section. Simulation results examining the proposed algorithms along with computational and comparative results are given and analyzed in Section 5. Simulations conducted by using the iMOPSE dataset [5] are presented in this section. In addition, this section describes experiments conducted on the factory’s dataset collected from TNG, a well-known national textile company. Finally, Section 6 concludes the paper and outlines some directions for future work.

2.1. Approximation Algorithms for MS-RCPSP

Scheduling problem arises in many practical situations [611]. In the general case, the scheduling is proved to be NP-hard. The authors of [5, 12, 13] show that MS-RCPSP, the problem we mentioned in this paper, is an NP-hard scheduling problem. In addition, MS-RCPSP is multimodal or discontinuous and very challenging to solve with traditional optimization methods. In the past few years, various approaches have been proposed to solve MS-RCPSP, thereby finding out intelligent algorithms to allocate and manage network resources in the Internet of things systems [1417]. Some popular metaheuristic algorithms are the genetic algorithm (GA), ant colony optimization (ACO), and particle swarm optimization [1822].

Su et al. [23] used a mixed-integer model and discrete constraints to solve the problems. Maghsoudlou et al. [24] and Bibiks et al. [25] applied the cuckoo Search algorithm to plan for the multirisk project with three distinct evaluation objectives. Lin et al. [26] proposed a new solution based on GA and some other heuristic algorithms.

Myszkowski et al. [5, 13] proposed the hybrid algorithm which is the combination of the differential evolution strategy and the greedy approach for arranging the human resources in product manufacturing projects and calculating time-of-use electricity tariffs and shift differential payments. The authors aimed to achieve the shortest makespan and low cost and gave a new benchmark dataset called iMOPSE [5], an artificially created dataset based on real-world instances for the MS-RCPSP.

The MS-RCPSP is studied and applied in the planning of many practical areas of the Internet of things. Some authors also study new variants of the MS-RCPSP and apply them in different fields. Mejia [3] developed a new variant of the MS-RCPSP to coordinate research activities in the nuclear laboratory. Hosseinian [27] proposed two new algorithms, P-GWO and MOFA, to solve the MS-RCPSP with the deterioration effect and financial constraints (a case study of a gas treating company). In another study, Hosseinian et al. introduced a new mixed-integer formulation for the time-dependent MS-RCPSP considering the learning effect [28].

Nemati-Lafmejani et al. [29] developed an integrated biobjective optimization model to deal with multimode RCPSP. The objective of the proposed model is to minimize both the costs and the makespan of the project. Tahrir and Yang [30] proposed a hybrid metaheuristic algorithm that combined ant colony optimization and variable neighbourhood search approaches for job scheduling in grid computing.

Several prior research has also been conducted to solve other subproblems of RCPSP, the general problem of MS-RCPSP, to improve the efficiency of network resource management in specific fields of the Internet of things. Ballestin and Leus [31] and Javanmard [32] and their colleagues studied specific cases of RCPSPs such as multiskill stochastic and preemptive problem to calculate project implementation time, i.e., makespan, and build mathematical models for the project resource investment.

2.2. Differential Evolution Method

Differential evolution (DE) is a population-based stochastic optimization algorithm introduced by Storn and Price [4]. The advantages of DE lies in its simplicity, efficiency, and speed thanks to the local search method. DE deals with the old generation of the original population by using the mutation operator to create better solutions in the new generation. Until now, DE is clearly recognized to be the approach that has the potential to solve NP-hard [4, 33] problems.

Tanabe et al. and Ghosh and Das [34, 35] proposed a DE-based method for scheduling problems in grid computing environment whose purpose is to minimize the execution time and add more parameter adjustment to get high effects.

Like other evolutionary approaches, DE performs the evolution of a population by using three kinds of the operator: crossover (recombination), mutation, and selection. The major difference between DE and genetic algorithm, a classical evolutionary algorithm, is due to the mutation operator. DE’s mutation operator uses orientation information to change solutions of the current population as follows.

Given the population that is composed of solutions, each of them consists of D components, and thus a particular solution is represented by a vector xi = (xi,2, xi,2,…, xi,D) where xi,j ∈ R, i = 1, 2, …, N; j = 1, 2, …, D.

To create the mutated solution of a given solution xi, DE picks out three random different solutions from the current population: xr1 ≠ xr2 ≠ xr3 ≠ xi, and then the mutation solution is determined as follows:where the value of the constant F, which plays the role of the mutant factor, is in the range [0, 1]. Because the mutant factor F is used to adjust the size of the directional vector, it is also called the directional hop length.

After executing the mutation step, the crossover operation is carried out by combining the parent solution xi and the mutation solution . The crossover operator is performed by selecting a random number CR (CR ∈ [0, 1]) as the probability of crossover. The result of the crossover step is represented by the vector (u1,1, u1,2, ... uN,D) where(i)i = 1, 2, …, N; j = 1, 2,…, D.(ii)randi,j: mutant factor.(iii)Irand: a random number varies within the range [1,D]. Thanks to Irand, the mutation solution is always different from the parent solution xi.

Finally, the next generation is created from the xi and the ui as follows:

Since 1997, various versions of DE have been developed by researchers. This paper proposes a variant of the DE algorithm which uses the following terminologies and conventions:(i)A solution represents a schedule, that is, a plan for performing given tasks(ii)The number of search space’s dimensions is equal to the number of tasks

3. System Model and Definitions

3.1. RCPSP

MS-RCPSP is just a special case of RCPSP. Before the outline of the multiskill RCPSP, the description of the classical RCPSP would be introduced as follows.(i)Assume a given project represented by a directed noncyclical graph G (V, E), where each node depicts a task and the arc represents finish-to-start precedence relationship between two tasks. The arc (u, ) ∈ E shows that the task u must be completed before task begins (Figure 1).(ii)Without any difference, two empty tasks were added to the project. The first one is placed at the beginning of the project as the predecessor of every other task, whereas the second empty task is placed at the end of the process as the successor of other tasks.(iii)The duration (also called execution time) of a certain task is calculated from the beginning time to the ending time of that task.(iv)The execution progress of every task cannot be stopped until it is done.(v)Each task needs an amount of resources for its execution.

The goal of the RCPSP is to find out a schedule of the task such that the precedence and resource constraints are satisfied and the makespan is minimized. RCPSP has been proved to be NP-hard [12]. Throughout the paper, we use the term “makespan” to refer to the execution time of the project, also called project duration. In the process of looking for the optimal schedule, both priority constraint between operations and resource constraint have to be met.

3.2. MS-RCPSP Formulation

In this paper, MS-RCPSP [5, 13, 36] was chosen because it is more practical than its general problem, RCPSP. MS-RCPSP uses the additional attribute, skill, in the problem formulation. Each task requires particular skills at the given levels to be executed, whereas each resource has some skills with given levels. Thus, the resource A is capable of performing the task B if A has skills required by B at the identical or higher level.

Before deepening the study of MS-RCPSP, we take a look at some notations as follows:(i)dj: the duration of task j; Pj: subset of the predecessors of task j(ii)qj: subset of skills required by task j(iii)Qk: the subset of skills owned by the resource k(iv)Q: the set of skills, Qk ⊆ Q(v)RS: the set of resources(vi)T: the set of tasks(vii)Tk: the set of tasks that resource k can perform, Tk ⊆ T(viii)RSk: the set of resources that can perform task k, RSk ⊆ RS(ix)Sj, Fj: start time and finish time of task j(x)Uj,kt: Uj,kt = 1 in case of resource k is assigned to task j at given time t; Uj,kt = 0 in other cases(xi)lq: the level of the given skill q; hq: type of the skill q(xii)qj: the skill required by task j(xiii)τ: duration of a schedule(xiv)PS: a feasible schedule(xv)PSall: the set of all feasible schedules(xvi)f(PS): makespan of schedule PS

Each schedule (PS) is responsible for coordinating n tasks and m resources.

Formally, MS-RCPSP could be stated as follows:where

4. Proposed Algorithm

4.1. Schedule Representation

In this study, every schedule or solution is represented by a vector that consists of elements; the number of elements (size of the vector) and the number of tasks are the same. The value of a certain element depicts the resource that can execute the corresponding task.

Example 1. Given a set of task T = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10} and a set of resource RS = {1, 2, 3}. The MS-RCPSP’s goal is to find out a schedule in which the makespan can be minimized while meeting priority constraints between tasks (Figure 1).
Consider a set of task T = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10} and set of resource RS = {1, 2, 3}. Priority constraints between tasks are depicted in Figure 1, whereby task 1 is the predecessor of task 2, and thus task 2 cannot be performed until task 1 is finished: S2 ≥ F1. There is no relation between task 3 and task 4; therefore, they are executed concurrently or sequentially.
Besides, in this example, we assume that(i)sk = sik, i ∈ RS(ii)Qk = Qk ∈ RS(iii)RSk = RS ∀k ∈ RSThe duration of the tasks is given in Table 1.
A feasible schedule for the project is demonstrated in Figure 2, where it can be seen that the makespan is 15.
The above schedule is demonstrated in Table 2, where resource 1 performed tasks number 1, 2, 6, 8, and 9; resource 2 is assigned to execute tasks number 3 and 7; resource 3 deals with the remaining tasks.


Task12345678910

Duration2435225342


Task12345678910

Resource1123312113

4.2. Measurement Model

The differential evolution algorithm was primarily designed for real-valued data. However, like most of the other scheduling problems, MS-RCPSP deals with discrete values, and thus there is a need for the discrete variant of DE. In addition, this discrete DE algorithm needs a model for representing the MS-RCPSP’s factors such as the difference between two schedules. In this paper, we build a new measurement model for measuring the difference between the two schedules as follows:(i)Unit vector P =(p1, p2, ..., pn); pi = 100/(ki − 1), where ki is the number of resources that can execute task i(ii)The difference between schedule X =(x1, x2, ..., xn) and schedule Y =(y1, y2, ..., yn) is denoted by the differential vector D = X – Y =(d1, d2, ..., dn)whereorder (xi): the position of the resource (xi) in the RSi(iii)The sum of the schedule Y and differential vector D is the schedule X =(x1, x2, ..., xn)whereposition (i): the resource corresponding to the position i

Example 2. Assume that there are 6 resources that could handle task 1: RS1 = {R1, R3, R4, R5, R9, R10}. Thus, k1 = 6; p1 = 100/(6–1) = 20. Assume that there are 5 resources that could handle task 2: RS2 = {R3, R5, R7, R8, R10 }. Thus, k2 = 5; p2 = 100/(5–1) = 25. The resource order in measurement model is presented in Table 3.
Consider 2 schedules: X =(3, 5); Y =(5, 10). The task assignment is shown in Table 4.
D = XY =(d1, d2), where d1 = p1. Abs (order (x1)order (y1)) = 20. Abs (1–3) = 40.
d2 = p2. Abs (order (x2)order (y2)) = 25. Abs (1–4) = 75. Thus, D =(40, 75). Consider D’=(35.71, 5.23). As shown in Table 5, we have Z = X + D’=(5, 8).


Order012345

ResourceR1R3R4R5R9R10
ResourceR3R5R7R8R10


ScheduleTask
12

XR3R5
YR5R10


Task12

XR3R5
D35.715.23
Pi. order (xi)2060
X + D55.7165.23
Z = X + D′R5R8

4.3. Proposed DEM Algorithm

The proposed algorithm (differential evolution multiskill, or DEM) is described in Algorithm 1:(i)Stop criterion: the loop while-end while would be stopped if the difference between makespan of some continuous generations is smaller than a threshold(ii)Size: number of solutions in the population(iii)F: parameter of the mutation operator(iv)CR: crossover probability(v)n: number of tasks

(1)Begin
(2) Load and validate iMOPSE dataset
(3)t 0
(4) Sizesize of the population
(5)P(t)Initial population
(6)f(t)Calculate the fitness and makespan
(7) while (Stop criterion)
(8)  for (int i = 0; i < Size; i++)
(9)   xr1 xr2 xr3 xirand (1,Size)
(10)   F rand (0, 1)
(11)   (t) ⟵ xr1 + F ×(xr xr3)// parameter of the mutation operator
(12)   for (j = 0; j < n; j++)
(13)    
(14)   end for
(15)   if (f (ui (t)) ≤ f(xi(t)))
(16)    xi (t+1) = ui(t)
(17)   else
(18)    xi (t+1) = xi (t)
(19)   end if
(20)  end for
(21)  Calculate the fitness and bestnest
(22)  If (makespan > min (fitness))
(23)  makespan = min (fitness)
(24)  End if
(25)  bestnest ⟵ Reallocate (bestnest)
(26)  t ⟵ t + 1
(27) end while
(28)return makespan
(29)End

As shown in line 25, after classical DE steps such as mutation and selection are executed, function Reallocate() is called to reduce the makespan of the candidate schedule bestnest. Procedures of this function will be dealt with in more detail in the next section (Algorithm 2):

currentBest: the best schedule among the population of the current generation.Rb : According to the arrangement of schedule newbest (line 3), Rb is the last resource that finishes its job. In Figure 3, before the reallocate function was called, Rb is resource 1.maxResource(): the function which returns the last resource to finish (line 3).size(): the function which returns the number of elements of the given set or array (line 5).newmakespan: makespan of the new schedule (line 11), which is obtained from the old schedule by reassigning task i. The task i is moved from resource Rb (line 9) to resource Ri[j] (line 10).

Input: currentBest //the best schedule among the current population
Output: //the improved schedule
(1) Begin
(2) makespan = f(best)
(3) newbest = currentBest;
(4) Rb ⟵ maxResource (newbest) //the last resource to finish its job
(5) Tb ⟵ set of tasks is performed by resource Rb
(6) For i = 1 to size(Tb) // Consider each task in Tb, the set of tasks performed by resource Rb
(7)  Ti = Tb [i];
(8)  Ri ⟵ set of resource that are skilled enough to execute the task i except Rb
(9)  For j = 1 to size (Ri) // Consider each resource in turn
(10)   Ri[j] = Ri[j] + {i} // Reassign task i to resource Ri[j]
(11)   Rb = Rb – {i} // Remove task i from the task list of Rb
(12)   newmakespan = f(newbest) // The makespan of the new schedule
(13)   If newmakespan < makespan
(14)   makespan = newmakespan
(15)   Return newbest; // Successful Reallocate, return the new and better schedule
(16)   End if
(17)   newbest = currentBest; // Unsuccessful Reallocate, reuse the original schedule for the next iteration
(18)   End for
(19)  End for
(20) Return best
(21) End Function
4.4. Function Reallocate

Line 12 and 13 demonstrate that the new schedule (newbest) is better than the old schedule (currentBest) in terms of makespan, which means that the function Reallocate certainly returns a better schedule. Moreover, the function Reallocate guides DEM algorithm in the right direction by finding a new schedule (newbest) based on the best solution (currentBest); therefore, it inherits the advanced genomes of the population as well as the old schedule.

Example 3. Given set of task T = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10}, set of resource RS = {1, 2, 3}. Resource 1 can execute tasks 1, 2, 3, 4, 6, 8, 9, and 10; resource 2 can perform tasks 1, 3, 7, and 9; resource 3 can handle tasks 1, 4, 5, 8, 9, and 10. The relationship between tasks is depicted in Figure 1, and the duration of the tasks is given in Table 1.
Consider a schedule (PS) as shown in Table 6. From observation, Figure 4 shows that makespan f(PS) = 17.
New schedule, the result of the function Reallocate, is described in Table 7.
Figure 4 illustrates that function Reallocate converts the schedule PS into the new schedule by assigning task 9 to resource 2 instead of resource 1 as before. Thanks to this allocation, the new schedule achieves better makespan (15) compared to the old schedule (17).


Task12345678910

Resource1123312133


Task12345678910

Resource1123312113

5. Experiment

Experiments were conducted on the simulated dataset and the real dataset to evaluate the performance of proposed algorithms in comparison with previous algorithms such as GreedyDO and GA. We first utilize the simulated iMOPSE dataset, a collection of project instances that are artificially created based on a database obtained from the international company, to evaluate the performance of the proposed algorithm.

5.1. Simulation on Simulated Dataset
5.1.1. iMOPSE Dataset

The simulated dataset iMOPSE [5] has been used to investigate existing algorithms GreedyDO and GA [13, 36]; therefore, we conduct simulations on iMOPSE dataset to fairly compare between DEM and those algorithms.

As shown in Table 8, iMOPSE’s instances are classified into 15 categories according to the following attributes:(i)The number of tasks and resources(ii)The number of precedence relations between tasks(iii)The number of skills


Dataset instanceTasksResourcesPrecedence relationsSkills

100_5_48_91005489
100_5_64_1510056415
100_5_64_91005649
100_10_64_910010649
100_10_65_15100106515
100_20_65_15100206515
100_20_65_910020659
200_10_84_920010849
200_10_85_15200108515
200_20_55_920020559
200_20_97_15200209715
200_20_97_920020979
200_40_45_920040459
200_40_90_920040909
200_40_91_15200409115

All simulations were run on a PC with Intel Core i5-CPU 2.2 GHz, 6 GB RAM, and operating system Windows 10. All of iMOPSE instances used by our simulations are listed in Table 8.

5.1.2. Parameters and System Settings

The simulation is conducted by using the following:(i)Simulation environment: Matlab ver. 2014(ii)Simulation tool: the performance of previous algorithms such as GA is fairly evaluated by using GARunner [13]; the simulation tool was provided by the author of those algorithms(iii)Input data: 15 iMOPSE instances that are described in Table 8(iv)Number of solution in the population Np: 100(v)Number of generations : 50,000(vi)Since all of the considered algorithms are approximate, we run on each instance 35 times to get average value, the standard deviation value, and the best value

5.1.3. Simulation Results

Simulation results are shown in Tables 9 and 10, where the makespan of the best solutions found by GreedyDO, GA, and the proposed algorithm (DEM) is illustrated (in hours). Table 9 shows that the makespans acquired by GA are always better than those acquired by GreedyDO.


Dataset instanceGreedyDOGA

100_5_48_9779528
100_5_64_15640527
100_5_64_9597508
100_10_64_9533296
100_10_65_15426286
100_20_65_15310240
100_20_65_9408181
200_10_84_9999567
200_10_85_15706549
200_20_55_9999312
200_20_97_15680424
200_20_97_9816321
200_40_45_9821209
200_40_90_9963211
200_40_91_15519200


Dataset instanceGADEM
AvgBestStdAvgBestStd

100_5_48_95355289.74984920.50
100_5_64_155305272.55105041.50
100_5_64_95215089.94884851.00
100_10_64_93052966.62642571.50
100_10_65_152902865.02632580.50
100_20_65_152402400.02112100.50
100_20_65_91871814.51731581.00
200_10_84_958356711.45375273.00
200_10_85_155555494.94924900.50
200_20_55_93183124.23102900.50
200_20_97_154384249.73363300.00
200_20_97_93263216.23143063.00
200_40_45_92132092.92172071.50
200_40_90_92152113.12172061.00
200_40_91_152052003.42051913.50

Table 10 illustrates the average value, the best value, and the standard deviation of the makespan found by GA and the proposed DEM algorithm.

From the observation, the result in Tables 9 and 10 shows that(i)The solutions found by DEM are better than Greedy from 16% to 78% and GA from 0% to 21%.(ii)Regarding the average value (Table 9): DEM is better than Greedy from 15% to 77.5% and GA from 0.7% to 20.8%.(iii)Inspired by the differential evolution strategy, the proposed algorithm uses two parameters to guide the generation of the solution generated in the next generation: (i) the experience vector of both the population and (ii) the experience vector of the solution across different generations. DEM not only inherits those advantages of differential evolution approach but also utilizes the function reallocate, which improves its performance. Therefore, DEM has faster convergence and greater stability, which is shown by the standard deviation in Table 9 (column Std). The total standard deviation of DEM is 19.5 while the value of GA is 55.7; this result proves that the stability of DEM algorithm is better than GA.(iv)In addition, the use of a scale to represent solutions helps DEM calculate the deviations between solutions more accurately and create better solutions in the next generations.

5.2. Simulation on Real Dataset
5.2.1. TNG Dataset

Experiments conducted on the simulated dataset are sometimes inaccurate and unsatisfactory. To overcome this drawback, we collected the dataset of TNG Investment and Trading Joint Stock Company [37], a well-known national textile company. TNG dataset instances have the following characteristics and parameters:(i)The company has product sewing contracts under the business partners' contracts; each order places a sample of a product in large quantities(ii)Each product has a certain number of stages, each of which requires a certain time to complete(iii)Each order will be assigned to a group of workers(iv)Each group consists of many workers(v)Each worker's skills are rated based on his rank

The above parameters of the dataset have been converted to a suitable format as shown in Table 11.


Real factor or parameterThe role

OrderProject
Stage of the productTask
Time to execute the stageDuration of the given task
WorkerResource
The worker gradeSkill level of the resource
The order in which the tasks are executedThe priority of the tasks
Order execution timeThe total time of project implementation, i.e., the makespan

Table 11 demonstrates that the TNG dataset parameters will be represented as follows:(i)The order from the business partners’ is represented by a project(ii)The stage of the product is represented by a task(iii)A worker in the factory is represented by a resource in the problem formulation(iv)The worker grade is represented by the skill level of the resource(v)The order in which the tasks are executed is represented by the priority of the tasks(vi)Order execution time is represented by the total time of project implementation, i.e., the makespan

After the conversion, the dataset is shown in Table 12


Dataset instanceNumber of tasksNumber of resourcesPrecedence relationsSkill levelsProject time

TNG1713710266409
TNG2713910266325
TNG3714110266296
TNG4714510266392
TNG513737189461174
TNG613739189461052
TNG71374118946871
TNG81374518946996

5.2.2. Parameters and System Settings

Experiment setting:(i)Input data: 8 instances that are listed in Table 12(ii)Number of the solution in the population Np: 100(iii)Number of generations : 50,000(iv)Like simulation on the simulated dataset, we run on each TNG dataset instance 35 times to get the best value

5.2.3. Simulation Results

To prove the effectiveness of the proposed algorithm (DEM), we carry out experiments by using the real-world dataset collected from the TNG company. DEM will be compared with other algorithms from previous research studies, GreedyDO and GA. The performance of these algorithms is recorded and observed by using GARunner, the tool created by Myszkowski et al. [13].

Table 13 contains the makespan (measured in hours) of the best solution achieved by GreedyDO, GA, and DEM. In addition, the performing intervals of the actual projects at the TNG factory are shown in column TNG.

Table 13 proves that schedules found by DEM are better than those acquired by GA and GreedyDO. Table 13 and Figure 5 also highlight that the schedules found by evolutionary algorithms are always shorter than the execution time of the actual projects at the TNG factory.


Dataset instanceTNGGreedyDOGADEM

TNG1409236201165
TNG2325243198166
TNG3296258212166
TNG4392248176165
TNG51174972751712
TNG61052963791712
TNG7871834810728
TNG8996906720675

Compared to the execution time at the TNG factory, the GreedyDO and GA algorithms reduce the makespan by 4%–42% and 7%–51%, respectively, while DEM has the best results of 16%–61%. This result is understandable because so far all projects have been manually scheduled based on the scheduler's experience.

This statistical analysis of the experimental result proves that DEM achieved better solutions than the previous algorithms. Thanks to the directed nature, DEM not only ensures fast convergence but also averts getting trapped in local extrema. In experiments conducted on the TNG dataset, DEM reduces the makespan from 16% to 61% compared to the current factory schedule.

6. Conclusion and Future Work

This paper proposes to investigate the MS-RCPSP, one of the combinatorial optimization problems that have a wide range of practical applications in science, technology, and life. The MS-RCPSP is formulated and presented as a mathematical model, and then a novel algorithm called DEM which is based on differential evolution strategy is proposed.

The MS-RCPSP, which is the subject of the investigation in this paper, has significant meaning for network resource management in the Internet of things. The MS-RCPSP we study takes into account the resource’s skill level, an important factor of the practical projects. Inspired by differential evolution metaheuristic, we have devised a new evolutionary algorithm that has the ability to outperform all the previous approaches. We build a function called reallocate, which is the key component that creates the strength of the proposed algorithm. This function avoids the idle time intervals in a given schedule, improving the solution quality.

To analyze the efficiency of the proposed algorithm, various simulations are conducted and compared with the previous algorithms. Besides the simulated dataset, the real dataset supported by the TNG company has also been utilized for the assessment. The simulation result shows that our developed algorithm is more effective than existing algorithms. DEM achieves not only better quality solution but also faster convergence to global extremum in comparison with the prior algorithms.

In the next steps, we are planning to improve DEM algorithm by using a blended random walk and the Gauss statistical function. We will implement local search operations to further improve the quality of the solutions. At the same time, we will integrate the Gaussian probability distribution to improve the search direction modification in the solution space. Thanks to that improvement, we expect that the proposed algorithm will not fall into the local extreme trap.

In the future, we will utilize the DE metaheuristic with the integration of the Gauss distribution and Gauss statistical function [38, 39] to improve the effectiveness of the searching for the optimum solution and avoid the local extreme trap.

Data Availability

The paper uses the standard iMOPSE dataset to test the efficiency of the algorithm. This dataset is publicly available at http://imopse.ii.pwr.wroc.pl/ and is free of charge. In addition, we also tested the algorithm with TNG's garment manufacturing dataset. We have obtained permission from TNG to use these data.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

We are grateful to the TNG Investment and Trading Joint Stock Company for giving us the information regarding their projects as well as providing us with their factory dataset.

References

  1. Z. Ning, G. Xu, N. Xiong et al., “TAW: cost-effective threshold authentication with weights for Internet of Things,” IEEE Access, vol. 7, pp. 30112–30125, 2019. View at: Publisher Site | Google Scholar
  2. M. Alirezaei, S. T. A. Niaki, and S. A. A. Niaki, “A bi-objective hybrid optimization algorithm to reduce noise and data dimension in diabetes diagnosis using support vector machines,” Expert Systems with Applications, vol. 127, pp. 47–57, 2019. View at: Publisher Site | Google Scholar
  3. O. P. Mejia, “A new RCPSP variant to schedule research activities in a nuclear laboratory,” in Proceedings of the 47th International Conference on Computers and Industrial Engineering (CIE47), Lisbon, Portugal, October 2017. View at: Google Scholar
  4. K. Price, R. Storn, and J. Lampinen, Differential Evolution-A Practical Approach to Global Optimization, Springer, Berlin, Germany, 2005.
  5. Y. Myszkowski, B. Paweł, and E. S. Marek, “A new benchmark dataset for multi-skill resource-constrained project scheduling problem,” in Proceedings of the 2015 Federated Conference on Computer Science and Information Systems (FedCSIS), IEEE, Sofia, Bulgaria, September 2015. View at: Google Scholar
  6. Y. Zeng, N. Xiong, and T. H Kim, “Channel assignment and scheduling in multichannel wireless sensor networks,” in Proceedings of the 2008 33rd IEEE Conference On Local Computer Networks (LCN), pp. 512-513, Montreal, Canada, October 2008. View at: Publisher Site | Google Scholar
  7. H. Cheng, N. Xiong, and L. T. Yang, “Distributed access scheduling algorithms in wireless mesh networks,” in Proceedings of the 22nd International Conference On Advanced Information Networking And Applications (Aina 2008), pp. 509–516, Okinawa, Japan, November 2008. View at: Publisher Site | Google Scholar
  8. G. Wei, A. V. Vasilakos, and N. Xiong, “Scheduling parallel cloud computing services: an evolutional game,” in Proceedings of the 2009 First International Conference on Information Science and Engineering, pp. 376–379, Nanjing, China, December 2009. View at: Publisher Site | Google Scholar
  9. W. Nie, N. Xiong, J. H. Park, and S. Yeo, “A fair-oriented two-level scheduling scheme for downlink traffic in wimax network,” in Proceedings of the 2010 2nd International Conference on Information Technology Convergence and Services, pp. 1–6, Cebu, Phillippines, January 2010. View at: Publisher Site | Google Scholar
  10. W. Zheng, N. Xiong, N. Ghani, M. Peng, A. V. Vasilakos, and L. Zhou, “Adaptive scheduling for wireless video transmission in high-speed networks,” in Proceedings of the 2011 IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS), pp. 180–185, Shanghai, China, January 2011. View at: Publisher Site | Google Scholar
  11. P. Dong, J. Xie, W. Tang, N. Xiong, H. Zhong, and A. V. Vasilakos, “Performance evaluation of multipath TCP scheduling algorithms,” IEEE Access, vol. 7, pp. 29818–29825, 2019. View at: Publisher Site | Google Scholar
  12. R. Klein, Scheduling of Resource Constrained Project, Springer Science Business Media, NewYork, NY, USA, 2000.
  13. P. B. Myszkowski, M. Laszczyk, I. Nikulin, and E. Skowronski, ““iMOPSE: a library for bicriteria optimization in multi–skill resource–constrained project scheduling problem,” Soft Computing, vol. 23, no. 10, pp. 3397–3410, 2019. View at: Publisher Site | Google Scholar
  14. H. Cheng, N. Xiong, X. Huang, and L. T. Yang, “An efficient scheduling model for broadcasting in wireless sensor networks,” in Proceedings of the 2013 IEEE International Symposium on Parallel & Distributed Processing, pp. 1417–1428, IEEE, Cambridge, MA, USA, May 2013. View at: Publisher Site | Google Scholar
  15. B. Lin, W. Guo, N. Xiong, G. Chen, A. V. Vasilakos, and H. Zhang, “A pretreatment workflow scheduling approach for big data applications in multicloud environments,” IEEE Transactions on Network and Service Management, vol. 13, no. 3, pp. 581–594, 2016. View at: Publisher Site | Google Scholar
  16. B. Lin, W. Guo, G. Chen, N. Xiong, and R. Li, “Cost-driven scheduling for deadline-constrained workflow on multi-clouds,” in Proceedings of the 2015 IEEE International Parallel and Distributed Processing Symposium Workshop, pp. 1191–1198, IEEE, Hyderabad, India, May 2015. View at: Publisher Site | Google Scholar
  17. L. Tan, Z. Zhu, F. Ge, and N. Xiong, “Utility maximization resource allocation in wireless networks: methods and algorithms,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, vol. 45, no. 7, pp. 1018–1034, 2015. View at: Publisher Site | Google Scholar
  18. P. Guo, W. Cheng, and Y. Wang, “A general variable neighborhood search for single-machine total tardiness scheduling problem with step-deteriorating jobs,” Journal of Industrial and Management Optimization, vol. 10, no. 4, pp. 1071–1090, 2014. View at: Publisher Site | Google Scholar
  19. S. Kavitha and P. Venkumar, “A vibrant crossbreed social spider optimization with genetic algorithm tactic for flexible job shop scheduling problem,” Measurement and Control, vol. 53, no. 1-2, Article ID 0020294019889085, 2020. View at: Google Scholar
  20. A. P. Agrawal, A. Choudhary, and A. Kaur, “An effective regression test case selection using hybrid whale optimization algorithm,” International Journal of Distributed Systems and Technologies, vol. 11, no. 1, pp. 53–67, 2020. View at: Publisher Site | Google Scholar
  21. W. Guo, J. H. Park, L. T. Yang, A. V. Vasilakos, N. Xiong, and G. Chen, “Design and analysis of a MST-based topology control scheme with PSO for wireless sensor networks,” in Proceedings of the 2011 IEEE Asia-Pacific Services Computing Conference, pp. 360–367, IEEE, Jeju, Korea, November 2011. View at: Publisher Site | Google Scholar
  22. X. Zhuang, H. Cheng, N. Xiong, and L. T. Yang, “Channel assignment in multi-radio wireless networks based on pso algorithm,” in Proceedings of the 2010 5th International Conference on Future Information Technology, pp. 1–6, Busan, Korea, November 2010. View at: Publisher Site | Google Scholar
  23. C. T. Su, M. C. Santoro, and A. Mendes, “Constructive heuristics for project scheduling resource availability cost problem with tardiness,” Journal of Construction Engineering and Management, vol. 144, Article ID 04018074, p. 8, 2018. View at: Publisher Site | Google Scholar
  24. H. Maghsoudlou, B. Afshar-Nadjafi, S. T. Akhavan Niaki, and S. Niaki, “Multi-skilled project scheduling with level-dependent rework risk; three multi-objective mechanisms based on cuckoo search,” Applied Soft Computing, vol. 54, pp. 46–61, 2017. View at: Publisher Site | Google Scholar
  25. K. Bibiks, Y.-F. Hu, J.-P. Li, P. Pillai, and A. Smith, “Improved discrete cuckoo search for the resource-constrained project scheduling problem,” Applied Soft Computing, vol. 69, pp. 493–503, 2018. View at: Publisher Site | Google Scholar
  26. J. Lin, L. Zhu, and K. Gao, “A genetic programming hyper-heuristic approach for the multi-skill resource constrained project scheduling problem,” Expert Systems with Applications, vol. 140, p. 112915, 2020. View at: Publisher Site | Google Scholar
  27. A. H. Hosseinian and V. Baradaran, “P-GWO and MOFA: two new algorithms for the MSRCPSP with the deterioration effect and financial constraints (case study of a gas treating company),” Applied Intelligence, vol. 50, pp. 2151–2176, 2020. View at: Google Scholar
  28. A. H. Hosseinian, V. Baradaran, and M. Bashiri, “Modeling of the time-dependent multi-skilled RCPSP considering learning effect,” Journal of Modelling in Management, vol. 10, 2019. View at: Google Scholar
  29. R. Nemati-Lafmejani, H. Davari-Ardakani, and H. Najafzad, “Multi-mode resource constrained project scheduling and contractor selection: mathematical formulation and metaheuristic algorithms,” Applied Soft Computing, vol. 81, Article ID 105533, 2019. View at: Publisher Site | Google Scholar
  30. Y. Tahrir and S. Yang, “Hybrid meta-heuristic algorithms for independent job scheduling in grid computing,” Applied Soft Computing, vol. 72, pp. 498–517, 2018. View at: Google Scholar
  31. F. Ballestín and R. Leus, “Resource-Constrained project scheduling for timely project completion with stochastic activity durations,” Production and Operations Management, vol. 18, no. 4, pp. 459–474, 2009. View at: Publisher Site | Google Scholar
  32. S. Javanmard, B. Afshar-Nadjafi, S. T. Akhavan Niaki, and S.T.K. Niaki, “Preemptive multi-skilled resource investment project scheduling problem: mathematical modelling and solution approaches,” Computers & Chemical Engineering, vol. 96, pp. 55–68, 2017. View at: Publisher Site | Google Scholar
  33. V. Feoktistov, Differential Evolution, Springer, New York, NY, USA, 2006.
  34. S. Tanabe, F. Ryoji, and A. Fukunaga, “Success-history based parameter adaptation for differential evolution,” in Proceedings of the 2013 IEEE congress on evolutionary computation, IEEE, Cancun, Mexico, November 2013. View at: Google Scholar
  35. T. K. Ghosh and S. Das, “A novel hybrid algorithm based on firefly algorithm and differential evolution for job scheduling in computational grid,” International Journal of Distributed Systems and Technologies, vol. 9, no. 2, pp. 1–15, 2018. View at: Publisher Site | Google Scholar
  36. H. Najafzad, H. Davari-Ardakani, and R. Nemati-Lafmejani, “Multi-skill project scheduling problem under time-of-use electricity tariffs and shift differential payments,” Energy, vol. 168, pp. 619–636, 2019. View at: Publisher Site | Google Scholar
  37. http://www.tng.vn TNG Investment and Trading Joint Stock Company, 434/1 bac kan street, Thai nguyen city, Viet Nam,.
  38. L. Kang, R.-S. Chen, N. Xiong, Y. C. Chen, Y.-X. Hu, and C.-M. Chen, “Selecting hyper-parameters of Gaussian process regression based on non-inertial particle Swarm optimization in internet of things,” IEEE Access, vol. 7, pp. 59504–59513, 2019. View at: Publisher Site | Google Scholar
  39. M. Wu, N. Xiong, and L. Tan, “Adaptive range-based target localization using diffusion gauss–newton method in industrial environments,” IEEE Transactions on Industrial Informatics, vol. 15, no. 11, pp. 5919–5930, 2019. View at: Publisher Site | Google Scholar

Copyright © 2020 Huu Dang Quoc 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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views211
Downloads265
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.