A new AGV scheduling algorithm based on harmony search for

0 downloads 0 Views 2MB Size Report
Then, an improved harmony search algorithm is proposed, which includes an effective ..... The serial number of each. CNC material buffer .... Time. 420. 440. 480. 520. 580. 660. 750. 820. 860. 910. Table 3. Case 3. No. 8. 19. 28. 6. 12. 20. 24.
Research Article

A new AGV scheduling algorithm based on harmony search for material transfer in a real-world manufacturing system

Advances in Mechanical Engineering 2018, Vol. 10(3) 1–13 Ó The Author(s) 2018 DOI: 10.1177/1687814018765560 journals.sagepub.com/home/ade

Guomin Li1, Bing Zeng1,2, Wei Liao1, Xinyu Li1 and Liang Gao1

Abstract With the development of advanced manufacturing technology, more and more enterprises utilize automated guided vehicles to transfer materials for computer numerical control machines in the workshop. Because automated guided vehicle scheduling problem is non-deterministic polynomial–hard, highly efficient automated guided vehicle scheduling algorithms are necessary to improve productivity of the workshop. In order to solve this problem, a mathematical model with a new objective function is established firstly, which incorporates two indicators, that is, the standard deviation of the waiting time of computer numerical control material buffers and the total travel distance of automated guided vehicle, reflecting the overall capacity of computer numerical control machines and the energy efficiency of automated guided vehicles during the actual production, respectively. Then, an improved harmony search algorithm is proposed, which includes an effective discrete encoding scheme of harmony, a new initialization method for harmony memory based on opposition-based learning strategy, a dynamic harmony memory considering rate parameter and a local search strategy. The property of the proposed algorithm is evaluated through three cases from the real-world manufacturing system for producing back cover of smart phone. Compared with the competing algorithms, the results demonstrate the superiority of proposed algorithm for solving automated guided vehicle scheduling problem obviously. Keywords Automated guided vehicle, computer numerical control machines, harmony search, meta-heuristic, material transfer

Date received: 10 December 2017; accepted: 20 February 2018 Handling Editor: Jose Antonio Tenreiro Machado

Introduction Automated guided vehicles (AGVs) are widely used to transport commodities and materials between workstations and warehouse.1 AGVs are among the most efficient methods for material handling due to their supreme routing flexibility, space utilization, product quality, and safety. The delivery demands of computer numerical control (CNC) machines are unpredictable and highly dependent on the product. In general, an optimized transfer scheduling can effectively decrease the waiting time of CNC machines and improve the production efficiency. Thus, optimally addressing the problem of AGV scheduling for

transferring materials to CNC machines (AGVSP) is instrumental to developing an enterprise into a manufacturing giant.

1

The State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan, China 2 XEMC Windpower Company Limited, Xiangtan, China Corresponding author: Bing Zeng, The State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan 430074, China. Email: [email protected]

Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).

2 At present, for almost all the enterprises, the transfer AGV scheduling for CNC machines is based on ‘‘First Come First Served (FCFS),’’ that is, transfer materials to CNC machines according to their sequence of calling materials. However, this scheduling method may not be the optimal solution; sometimes, it tends to seriously compromise the production efficiency of enterprises. The count of feasible scheduling solutions increases exponentially with the number of CNC machines, so the AGVSP is a non-deterministic polynomial (NP)-hard problem. Therefore, it is practically impossible to use a brute-force approach for identifying the optimal solution to this problem. However, it is feasible to solve this problem with meta-heuristic algorithms. To achieve a significant improvement concerning the overall performance of AGVSP, the scheduling process must consider the interdependencies between CNC machines and AGVs. In this article, we research on a part of the real complex workshop where there has been single AGV feed for cells and an optimal methodology for solving AGVSP is proposed. Due to the discreteness of AGVSP, harmony search (HS) algorithm is one of the most appropriate meta-heuristics for settling the AGVSP. By a judicious AGV scheduling from three cases in a real-world manufacturing system for producing back cover of smart phone, the standard deviation of the waiting time of CNC material buffers can be reduced to a great extent, that is, the overall capacity of CNC machines can be improved significantly. Moreover, this effort can also potentially shorten an AGV’s travel routes, that is, improve the energy efficiency of AGV.

Related works Dispatching AGVs to transport the assigned jobs to destinations is a kind of AGV scheduling. An appropriate route should be searched when the AGV scheduling is planned. The AGV scheduling problems belong to a kind of vehicle routing problem (VRP) or sequencing problem.2,3 The research of the scheduling problem by French4 is classic. They not only introduced the jobshop scheduling but also used material to illustrate the ideas behind combinatorial optimization. A variety of literatures researched on AGV scheduling problems for different applications. Soylu et al.5 developed an artificial neural network algorithm to seek the shortest routing for the single, multi-capacitated AGV carrying out pick and deliver requests. Tsirimpas et al.6 proposed dynamic programming algorithms and illustrated three practical variations to find the optimal routing of the single, limited capacity vehicle transporting the assigned product to the appointed client. Minis and Tatarakis7 proposed a dynamic programming algorithm to minimize the routing cost of the single vehicle delivering random products and picking up items. Hu

Advances in Mechanical Engineering et al.8 combined the static problem model and solution approach with a rolling-horizon approach to research the dynamic routing problem. In the matter of workshop layout, the multi-capacitated Rail Guided Vehicle (RGV) working on a linear track and the two-sided loading/unloading operations are similar to ours. There is also a resemblance between the objective of minimizing the energy consumption of RGV and the total travel distance of AGV. However, they regarded the routing problem as a mixed-integer linear program. Qudeiri et al.9 proposed an artificial neural network approach to predict the products’ routes’ group and then combined with genetic algorithms (GAs) to optimize proper routes for all product types in flexible manufacturing systems. Gultekin et al.10 developed a mathematical model for the scheduling problem of flexible robotic cell considering fixed process time and sequence for the CNC machines. They regarded the problem as a special traveling salesman problem. Nejad et al.11 developed a simulated annealing algorithm to study the problem of transferring materials to machines by robot, which is similar to this study. They considered the scheduling problem of a flexible robotic cell in which machines were identical and parallel, meanwhile the sequence of the actions was also determined, which was regarded as the traveling salesman problem. The objective is to minimize the cycle time of the robot. The above literatures research on AGV scheduling problem, but they are single objective optimization problems; however, the single-objective scheduling model cannot fully address the requirements of the realworld production system. To the best of our knowledge, few researches have been performed to solve the AGVSP using HS algorithm. In this article, we investigate a multi-objective AGVSP problem derived from a realistic manufacturing case and then propose an Improved Harmony Search–based Material Transfer Scheduling Algorithm (IHSMTSA) to solve AGVSP. The remainder of this article is organized as follows. The related works are elucidated in section ‘‘Related works.’’ Section ‘‘Mathematical model of AGVSP’’ describes the mathematical model of AGVSP. The proposed material transfer algorithm with sufficient detail is presented in section ‘‘Proposed algorithm of AGVSP.’’ In section ‘‘Experimental results and analysis,’’ an empirical evaluation through computational experiments is conducted. Section ‘‘Conclusion’’ presents the conclusions and outlines the potential future direction of this research.

Mathematical model of AGVSP Problem description The production line in the workshop illustrated in Figure 1 is popular at present. The workshop is

Li et al.

Figure 1. Schematic diagram of a CNC production line with cells and buffers.

composed of multiple CNC machine cells, and each cell contains multiple CNC machines. The production and processing operation is carried out in the form of CNC machine cells and each cell sets a buffer. Materials needed for CNC machine cells will be delivered by AGV. Then, when CNC machines are brought into production operations, materials will be obtained from the buffer. In this article, we research on one of the production lines in the workshop shown in Figure 1, which is briefly depicted in Figure 2. Here, we only study the optimization problem of AGV transferring materials to the buffer of each CNC machine cell, which aims to reduce the standard deviation of the waiting time of CNC material buffers and the total travel distance of AGV. With the purpose of avoiding wide variance between cells, the standard deviation of waiting time of CNC material buffers is considered. The risk of starved feeding will increase by the wide variance, which will influence capacity equalization of cells and then reduce overall productivity. Assuming a CNC production line is shown in Figure 2, and the CNC buffers’ sequence of calling

3 material is (1, 10, 11, 19, 2, 20), if an AGV feeds CNC material buffers in this order, it would need to travel repeatedly from one end of the track to the other. Consequently, it is clear that this transfer scheduling approach will make the total waiting time of the buffers greatly extended because the AGV spends too much time on traveling. More seriously, if a CNC material buffer waits too long, some CNC machines in this cell may run but lack of materials will reduce the productivity greatly. According to the specific situation of material calling depicted in Figure 2, a reasonable solution should make AGV feed the CNC material buffers at one end of the track followed by moving to the other end of the track to feed other CNC material buffers, a process being able to greatly reduce the total waiting time of the CNC material buffers so as to improve the productivity. However, during the actual machining, the CNC buffers’ sequence of calling material is dynamic and unpredictable. Therefore, based on the above introduction and marked signs, the problem that we will research in this article can be described as how to arrange scientifically the optimal sequence of tasks of having sent out demand signals for materials within record time slice (RTS) applied an effective modern heuristic algorithm and then produce the delivery sequence of AGV. RTS stands for a period of time. In other words, tasks of sending out demand signals for materials by buffers during this time period are researched in this article, and for the sake of research, RTS is assigned as 1000 in which there has been enough cells that sent out demand signals for materials. However, it must be specified that the method presented and validated in this article can also be applied to other situations when RTS is assigned other values.

IHSMTSA optimization model As mentioned above, the standard deviation of the waiting time of CNC material buffers and the total travel distance of the AGV are the two indicators should be considered when solving AGVSP. Therefore, we propose the mathematical model as follows

Figure 2. Schematic diagram of CNC material buffers calling material.

Min f (X ) = w1 f1 + w2 f2

ð1Þ

4

Advances in Mechanical Engineering sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi n 1X f1 = (wT(xi )  wTaver )2 n i=1 f2 =

n X

abs(pos(xi )  pos(xi1 ))

ð2Þ

ð3Þ

4. 5. 6.

Time needed for conducting a single loading of AGV to a buffer (s): Tsl = 30. The sum of the algorithm running time (s): rT. RTS = 1000 s.

i=1

8 absðposðx1 Þ  posðx0 ÞÞ > > > + Tsl < RTS  calltimeðx1 Þ + rT + Vagv wTðxi Þ = > > calltimeðx Þ  calltimeðx Þ + wTðx Þ + absðposðxi Þ  posðxi1 ÞÞ + T > : i1 i i1 sl Vagv wTaver =

n 1X wT(xi ) n i=1

pos(x0 ) = pos(AGV )

ð5Þ ð6Þ

The objective function (1) is to minimize the standard deviation of the waiting time of CNC material buffers and the total travel distance of the AGV. X = (x1 , x2 , . . . , xi , . . . , xn ) denotes a solution, which represents a sequence of CNC material buffers fed by AGV; f1 denotes the standard deviation of the waiting time of CNC material buffers corresponding to the solution X ; f2 represents the total travel distance of the AGV; and w1 and w2 are the weights of f1 and f2 , respectively. The values of f1 and f2 are calculated in equations (2) and (3), respectively. Equation (4) presents the waiting time of CNC material buffers calling material. Where, calltime(xi ) denotes the time of the buffer xi calling material, pos(xi ) indicates the position of the buffer xi , and pos(AGV ) indicates the current location of the AGV. It should be noted that the location of the AGV needs to be updated in accordance with the transfer operation. It can be seen from equation (4) that the waiting time of the current buffer consists of five parts: The first part denotes the duration from the calling time to the end of RTS, the second part corresponds to the algorithm running time, the third part is the waiting time of the former buffer calling material, the fourth part is the time it takes for AGV to travel to the current buffer, and the fifth part is the loading time. Equation (5) is the average waiting time of CNC material buffers calling material. wTaver is the average waiting time of CNC material buffers calling material. Equation (6) presents that x0 is the initial location of AGV. The conditions associated with the AGVs-CNC system are specified in the following: 1. 2. 3.

The distance between the buffer centers of two neighboring CNC cells on the same track (m): 5.5. The initial position of AGV (m): iniPos. The average moving speed of AGV (m/s): Vagv = 0.45.

if i = 1 ð4Þ otherwise

Proposed algorithm of AGVSP Brief introduction of HS algorithm As a new meta-heuristics search algorithm, the HS algorithm12 is impressed by the searching process of a perfect harmony for music. HS employs a suite of newness stochastic derivatives employed to the discrete variables that leverages a musician’s experience to identify the searching direction so as to ensure the convergence of the algorithm. It does not need differential gradients and initial guesses to be applied to the discrete and continuous variables. The HS algorithm can also extract the optimal solution when it is in the local optima. It is able to overcome the shortcoming associated with the GA’s theory in building module and take the relationship into consideration definitely in an integration operation. HS algorithm is increasingly used for optimization of problem-solving schemes in many fields, such as a vehicle dispatching problem in the NC workshop which is equipped with the just-in-time (JIT) distribution,13 a flexible job-shop scheduling problem,14 the comprehensive production and transport scheduling problem associated with the make-to-order (MTO) manufacturing,15 the routing problem of wireless sensor network (WSN),16 the optimization of continuous problems,17 realization of two-dimensional (2D) nearest neighbor quantum circuits,18 discrete sizing optimization of truss structures,19 fine-tuning deep belief networks,20 the clustering problem of WSN,21 image reconstruction from projections,22 and single area unit commitment problem.23 HS first initializes the harmony memory, followed by populating the memory with randomly generated new harmony. Only when the newly generated harmony is better than the worst harmony in the original harmony memory, the new harmony is put in the memory and replaced with the worst harmony of the original memory. This replacement process repeats itself until the stop criterion is met. The detailed workflow of HS is outlined in the following:24

Li et al.

5

Step 1. Initialize the algorithm parameters and specify the optimization problem The size of the harmony memory size (HMS) refers to the number of solution vectors in HM and is usually set to 5. The algorithm parameters like the harmony memory considering rate (HMCR), pitch adjusting rate (PAR), bandwidth (bw), and the stopping criterion (maximum iteration number) are applied to explore better solution vectors and are all defined in Step 3. The optimization problem is defined as Minimize f ð X Þ subject to xi 2 Xi , i = 1, 2, . . . , n

Step 2. Initialize the harmony memory HM with random values drawn from vectors ½X L , X U  The HM, as shown in equation (8), is a memory space storing solution vectors, that is, harmony, which is generated randomly following the uniform distribution 6 6 6 6 6 HM = 6 6 6 6 6 4

X1 X2 .. . Xi .. .

3 7 7 7 7 7 7 7 7 7 7 5

XHMS x1, 1 6 6 x2, 1 6 6 . 6 .. 6 =6 6 xi, 1 6 6 . 6 .. 4 xHMS, 1 2

x1, 2 x2, 2 .. .

  .. .

x1, j x2, j .. .

  .. .

x1, n x2, n .. .

xi, 2 .. .

 .. . 

xi, j .. .

 .. . 

xi, n .. .

xHMS, 2

xHMS, j



x0 i 2 fx1, i , x2, i , . . . , xHMS, i g x0 i 2 Xi

3 7 7 7 7 7 7 7 7 7 7 7 5

xHMS, n ð8Þ

Step 3. Improvise a new harmony vector A new harmony X 0 = ðx0 1 , x0 2 , . . . , x0 i , . . . , x0 n Þ is produced following two steps, that is, harmony memory considerations and pitch adjustments. The first step is shown as follows

if P1 \HMCR else

ð9Þ

where HMCR is a constant between 0 and 1, and P1 is a random number generated from the linear space spanning ½0, 1. Equation (9) indicates that x0i should be selected from the ith column of HM, that is, fx1, i , x2, i , . . . , xHMS, i g, if P1 is smaller than HMCR; otherwise, x0i can be a feasible value belonging to X l . We then need to determine whether or not to conduct pitch adjustment on these new harmonies selected from the HM. The criterion is shown as follows

ð7Þ

where X denotes the set of designed variables; f (X ) denotes the objective function; X l denotes the set of possible values of xi , which, in the case of continuous  with xLi optimization, follows X i = xi xLi  xi  xU i U being the lower bounds and xi being the upper bounds of xi respectively; and n denotes the number of designed variables.

2

x0i

Pitch adjusting decision for

x0i



Yes if ðP2 \ PARÞ No else ð10Þ

where PAR is a constant between 0 and 1 and P2 is a random number between 0 and 1. The formula above indicates that if PAR is greater than the random number P2 , the new harmony x0 must be pitch adjusted; otherwise, the original stays. In the case where a pitch adjustment is needed, x0 would be subject to the following process x0i

x0i 6 rand  bw

ð11Þ

where bw denotes the arbitrary distance bandwidth and rand is a random number generated between 0 and 1 following uniform distribution. Step 4. If the newly improvised harmony xi is superior to the worst harmony in HM, the former would replace the latter. Step 5. Terminate if the current iteration number exceeds MAXIMP. Otherwise, repeat Steps 3 and 4. At last, the optimal solution is derived from the best harmony in HM.

The proposed IHSMTSA for AGVSP The IHSMTSA for solving AGVSP is proposed in this section. Figure 3 describes the specific flow of IHSMTSA, which shows that the key steps of IHSMTSA include the initialization of HM, the process of the improvisation for a new harmony, and the local search. This section describes the encoding of harmony first and subsequently illustrates the initialization of HM and the improvisation of a new harmony. Finally, the local search scheme is presented. Encoding of harmony. Considering that AGVSP is a discrete optimization problem, it is more appropriate to use discrete coding method, which can greatly narrow down the search range in the solution space. We adopt

6

Advances in Mechanical Engineering

Figure 3. Flowchart of IHSMTSA.

the following vector for representing the ith harmony in the HM Xi = (xi, 1 , xi, 2 , . . . , xi, j , . . . , xi, n )

ð12Þ

where Xi denotes the ith harmony of HM, n denotes the length of the harmony, that is, the total number of CNC machines which need to be fed, and xi,j represents the jth pitch of harmony Xi . Since the data associated with the material calling differ from time slice to time slice, the following data structure is defined to facilitate the encoding of algorithm (Figure 4). A callMaterial type array named call is defined to record the information pertaining to the material called by CNC material buffers. The serial number of each CNC material buffer calling material and the corresponding time at which the material is called are saved

Figure 4. Structure type.

in the array call following the order in which the materials are called, as shown in Figure 5. Supposing the CNC buffers’ sequence of calling material and the corresponding time at which the material is called associated within the RTS are stored in the array call as shown in Figure 5. The tone in a harmony corresponds to the subscript number of the elements in the array call. For example, the ith harmony of HM is

Li et al.

7

Figure 5. Structure of the array call.

Xi = (2, 3, 5, 7, 8, 4, 0, 9, 6, 1), indicating that the order of AGV loading for CNC machines is (17, 11, 19, 13, 10, 6, 8, 14, 2, 4). Initialization of harmony. As a new method, oppositionbased learning (OBL) can improve the performance of machine learning at a large extent, which is proposed by Tizhoosh.25 OBL increases the probability of seeking the global optimal solution by extending the domain of the initial solution in the searching area. Building the initial HM by OBL can therefore assist in the identification of a more appropriate initial solution, which is illustrated in a previous study. Supposing P = (x1 , x2 , . . . , xi , . . . xD ) is a D-dimensional solution vector, and xi = ½xLi , xU i , i = 1, 2,    , D, the opposing solution vector OP = (x01 , x02 , . . . , x0i , . . . , x0D ) is defined as follows x0i = xLi + xU i  xi

ð13Þ

Furthermore, the harmony i is defined as Xi (xi, 1 , xi, 2 , . . . , xi, j , . . . , xi, n ), whose opposing solution vector is OXi = (x0i, 1 , x0i, 2 , . . . , x0i, j , . . . , x0i, n ). The element x0i, j in the opposing solution vector can be computed from equation (14) x0i, j = call:size  1  xi, j

ð14Þ

Figure 6. Strategy for dynamically adjusting HMCR.

dynamic adjustment of HMCR is conducted to prevent IHSMTSA from falling into the local optima in the early stage, giving rise to an enhanced global search ability. Meanwhile, the ability of local search in the late stage is also reinforced. An excessively small value of HMCR may compromise the local search ability of the HS algorithm, making it difficult to converge to the global optimal solution; an excessively large value of HMCR, however, may lead to an overly strong local search ability, rendering the algorithm prone to land on the local optimal solution instead of converging to the global optimal solution. Therefore, in this article, the HS algorithm conducts a dynamic update to the HMCR parameter. In the early stage, HMCR takes on a relatively small value to ensure a good global search ability. In the end, HMCR is assigned a relatively large value to improve the local search ability, which facilitates the convergence to the global optimal solution. The relationship between HMCR and the number of iterations is shown in equation (15)

 3ðHMCR HMCRðiterÞ =

max HMCRmin Þ itermax HMCRmax

where call:size denotes the number of elements in the array call. Improvisation of a new harmony. To apply the original classic HS optimization algorithm for solving a continuous optimization problem, one can use the Tone Tuning Bandwidth bw as a tuning parameter to adjust the harmonies selected from the HM. The AGVSP, however, is a discrete optimization problem, which may lead to tremendous amount of operational cost for adjusting harmonics based on the traditional HS algorithm. As a result, we choose to skip the step of pitch adjustment while a new harmony is improvised. Consequently, HMCR becomes the sole algorithmspecific parameter in the IHSMTSA. Specifically, the

 iter + HMCRmin ,

  if iter\ iter3max otherwise

ð15Þ

where itermax denotes the maximum number of iterations, iter represents the current iteration number, and HMCRmin and HMCRmax denote the minimum and maximum values of HMCR, respectively. In general, the values of HMCRmin and HMCRmax are set to be 0.2 and 0.99, respectively. The relationship between HMCR and the iteration number is shown in Figure 6. Supposing X 0 = (x01 , x02 , . . . , x0i , . . . , x0n ) is a new harmony, the element x0i of X 0 can be generated using equation (16) ( x0i

x0i !fx1, i , x2, i , . . . , xHMS, i g

if (P1 \HMCR)

x0i

Otherwise

rd

!X rd

i

ð16Þ

8

Advances in Mechanical Engineering

where ! denotes the process of random selection.

crossover, two bit-swap mutation, and local search of inserting a gene into another position; and HS algorithms. The CNC production line used in our experiments in shown in Figure 8, which includes 30 CNC machine cells, that is, 30 material buffers and 1 AGV. To make the comparison comprehensive, three cases from the real-world manufacturing system for producing back cover of smart phone are tested here involving the aforementioned five algorithms.

Equation (16) dictates that x0i is randomly selected from fx1, i , x2, i , . . . , xHMS, i g in the case where the random number P1 is smaller than HMCR; otherwise, x0i is chosen from X i randomly.

Case 1. There are 15 CNC material buffers calling material within the RTS, whose information is shown in Table 1.

Figure 7. Schematic diagram of local search.

rd

Local search. As a heuristic algorithm for solving optimization problems, local search is performed within the vicinity of the current solution. Specifically, it searches new solution within the neighborhood of the current solution. When the new solution is superior to the current solution, the former will substitute the latter. In view of above features of the AGVSP, the local search in IHSMTSA is implemented by moving a randomly selected tone to a new position in the harmony, as shown in Figure 7. If a new harmony with an improved fitness value is identified, it will replace the existing one. With this local search strategy, the quality of the optimal solution will be improved greatly.

Experimental results and analysis Experimental settings The experimental results are presented in this section. The algorithms compared here are all implemented using C + + programming language by Microsoft visual studio 2015 and executed on a personal computer with 2.5 GHz Intel Core i7 3610QM processor, 8 GB RAM and Microsoft Windows 10 operating system. In order to verify the superiority of the proposed IHSMTSA, the comparisons are conducted by incorporating the particle swarm optimization (PSO) with inertia weight; GA with elite selection; single-point crossover and two bit-swap mutation strategies; memetic algorithm (MA) with elite selection, single-point

Case 2. There are 20 CNC material buffers calling material within the RTS, whose information is shown in Table 2. Case 3. There are 25 CNC material buffers calling material within the RTS, whose information is shown in Table 3.

Parameter settings The parameter value of AGVs-CNC system and the experiments are listed in Table 4. It is worth noting that the experimental results were obtained by taking an average of 25 independent runs for fairness. The parameters of IHSMTSA and the comparison of algorithms are set as shown in Table 5.

Results and analysis To test the efficiency of the present algorithm, the IHSMTSA is compared with other meta-heuristicbased transfer scheduling algorithms in terms of the convergence rate for the cases 1–3, respectively. Figures 9–11 show the convergence behaviors of all the algorithms, plotting the average fitness value and the optimal fitness value as a function of evaluation number for each algorithm. Tables 6–8 show the standard deviation of the waiting time of CNC material buffers,

Figure 8. Schematic diagram of the CNC production line used in the present experiments.

Li et al.

9

Table 1. Case 1. No.

1

16

28

21

2

8

3

5

15

30

18

19

12

9

6

Time

10

150

160

190

250

280

450

560

650

760

820

880

950

980

990

Table 2. Case 2. No.

2

14

30

16

7

19

29

5

20

11

Time No. Time

6 1 420

50 21 440

80 8 480

100 26 520

150 3 580

200 6 660

260 17 750

280 21 820

320 22 860

380 15 910

Table 3. Case 3. No.

8

19

28

6

12

20

24

2

30

25

13

3

1

Time No. Time

10 29 451

48 22 482

82 14 522

110 7 546

146 26 587

180 10 622

212 23 678

250 16 755

278 5 787

325 15 825

356 9 866

388 18 970

421

Table 4. Parameter settings of experiments.

Table 5. Parameter settings of algorithms.

Items

Value

Algorithms

Parameters

Distance between two buffer centers iniPos Vagv Tsl rT RTS Independent runs FEs w1 w2

5.5 (m) 0 (m) 0.45 (m/s) 30 (s) 10 (s) 1000 (s) 25 10,000 0.7 0.3

IHSMTSA HS PSO

HMS = 6, HMCRmin = 0.2, HMCRmax = 0.99 HMS = 6, HMCR = 0.95, PAR = 0.05 PopSize = 100, vmax = 0.9, vmin = 0.1, c1 = 0.9, c2 = 0.1 PopSize = 100, Pc = 0.9, Pm = 0.01 PopSize = 100, Pc = 0.9, Pm = 0.01

RTS: record time slice.

the average waiting time of CNC material buffers, and the total travel distance of AGV. Case 1: the experimental results of 15 buffers calling material. Figure 9(a) shows the average fitness value versus FEs associated with five algorithms in case 1. The solution found by IHSMTSA is much better than those found by other algorithms including PSO, MA, and GA. Although the HS converges faster than IHSMTSA in early iterations as HS tends to get trapped in local optima, the HS converges slower than IHSMTSA in late iterations. Compared to the rate of convergence, it is of more importance to examine the quality of the solution. We can get the summary that IHSMTSA yields the highest probability of searching

GA MA

IHSMTSA: Improved Harmony Search based Material Transfer Scheduling Algorithm; HMS: harmony memory size; HMCR: harmony memory considering rate; PAR: pitch adjusting rate; HS: harmony search; PSO: particle swarm optimization; GA: genetic algorithm; MA: memetic algorithm; PopSize: population size; vmax: upper limit of inertia weight; vmin: lower limit of inertia weight; c1, c2: acceleration factors; Pc: crossover probability; Pm: mutation probability.

the best solution among the algorithms assessed here. Figure 9(b) shows the convergence curve of the best fitness value. The first population fitness value of IHSMTSA is worse than those of other algorithms due to its limited population capacity which makes it yield a lower probability of obtaining good solutions. However, the IHSMTSA can eventually converge to a more optimal solution than other algorithms after an increased number of iterations due to the fact that the dynamic changed HMCR can avoid falling into the local optima in early iterations, giving rise to an enhanced global search ability, and improve the ability of local search in later iterations. Besides, the proposed local search strategy has a great contribution on

10

Advances in Mechanical Engineering

Figure 9. Convergence behaviors of various algorithms in case 1: (a) average fitness value versus FEs and (b) best fitness value versus FEs.

Table 6. Comparison of algorithms in case 1. Methods

Standard deviation of the waiting time of CNC material buffers

Average waiting time of CNC material buffers (s)

Total travel distance of AGV (m)

IHSMTSA HS PSO MA GA FCFS

149.712 156.735 228.282 154.56 150.959 150.773

1146.89 1155.43 1155.56 1206.63 1174.07 1404.37

313.06 316.58 346.28 317.46 326.26 555.5

AGV: automated guided vehicle; CNC: computer numerical control; IHSMTSA: Improved Harmony Search based Material Transfer Scheduling Algorithm; HS: harmony search; PSO: particle swarm optimization; GA: genetic algorithm; MA: memetic algorithm; FCFS: First Come First Served. Note: The significance of the values given in boldface is reflected in the experimental analysis paragraph of case 1.

improving the quality of the optimal solution. Based on case 1 results, it can be concluded that IHSMTSA is superior to other four algorithms in terms of efficiency and accuracy. In this case of AGVSP, the comparison involving IHSMTSA, HS, PSO, MA, GA algorithms and FCFS has also been conducted. Table 6 illustrates the statistical results associated with the five algorithms and FCFS in case 1. The data used in FCFS are collected from engineering project in order to examine the science and validity of the proposed algorithm. Meanwhile, the data used in other cases are also credible. As we can see from Table 6, the proposed IHSMTSA performs best in terms of all the indicators, that is, the standard deviation of the waiting time of CNC material buffers, the average waiting time of CNC material buffers, and the total travel distance of AGV. Especially, the average waiting time of CNC material buffers obtained by IHSMTSA has been decreased by 18.33% compared with those obtained by FCFS, which means that IHSMTSA performs far better than FCFS in terms of the overall capacity. Moreover, the total travel distance of AGV obtained by IHSMTSA has been decreased by 43.64% compared with those obtained by FCFS, which

means that IHSMTSA can improve the energy efficiency of AGV significantly. Case 2: the experimental results of 20 buffers calling material. As we can see from Figure 10, IHSMTSA performs far better than other algorithms compared in terms of the rate of convergence and the accuracy of the best solution, similar to case 1. It can be seen from Table 7 that the standard deviation of the waiting time of CNC material buffers, the average waiting time of CNC material buffers, and the total travel distance of AGV associated with IHSMTSA are all shorter than those of the other four algorithms and FCFS. Especially, the values of the three indicators obtained by IHSMTSA have been decreased by 36.79%, 33.55%, and 76.7% compared with those obtained by FCFS, indicating that the advantages of IHSMTSA can stand out when solving larger scale AGVSP. Case 3: the experimental results of 25 buffers calling material. As we can see from Figure 11, IHSMTSA also performs far better than other algorithms compared in terms of the rate of convergence and the accuracy of the best solution, similar to case 1 and case 2, which

Li et al.

11

Figure 10. Convergence behaviors of various algorithms in case 2: (a) average fitness value versus FEs and (b) best fitness value versus FEs.

Figure 11. Convergence behaviors of various algorithms in case 3: (a) average fitness value versus FEs and (b) best fitness value versus FEs.

Table 7. Comparison of algorithms in case 2. Methods

Standard deviation of the waiting time of CNC material buffers

Average waiting time of CNC material buffers (s)

Total travel distance of AGV (m)

IHSMTSA HS PSO MA GA FCFS

161.959 186.565 326.832 215.986 194.398 256.223

1061.03 1143.38 1292.59 1179.49 1158.83 1596.76

138.38 180.18 314.38 208.34 216.7 594

AGV: automated guided vehicle; CNC: computer numerical control; IHSMTSA: Improved Harmony Search based Material Transfer Scheduling Algorithm; HS: harmony search; PSO: particle swarm optimization; GA: genetic algorithm; MA: memetic algorithm; FCFS: First Come First Served. Note: The significance of the values given in boldface is reflected in the experimental analysis paragraph of case 2.

further confirms the advantages of IHSMTSA on the rate of convergence and accuracy. It can be seen from Table 8 that the values of the three indicators associated with IHSMTSA are far smaller than those of the other four algorithms and FCFS. Especially, the values of the three indicators obtained by IHSMTSA have

been decreased by 58.54%, 28.24%, and 67.55% compared with those obtained by FCFS, which further the advantages of IHSMTSA on solving larger scale AGVSP. By running the program 25 times in three cases, separately, the optimal sequences of CNC material buffers

12

Advances in Mechanical Engineering

Table 8. Comparison of algorithms in case 3. Methods

Standard deviation of the waiting time of CNC material buffers

Average waiting time of CNC material buffers (s)

Total travel distance of AGV (m)

IHSMTSA HS PSO MA GA FCFS

160.854 228.171 439.593 250.121 236.749 388.008

1277.79 1323.59 1491.77 1325.58 1308.04 1780.56

223.08 253 436.7 284.24 272.58 687.5

AGV: automated guided vehicle; CNC: computer numerical control; IHSMTSA: Improved Harmony Search based Material Transfer Scheduling Algorithm; HS: harmony search; PSO: particle swarm optimization; GA: genetic algorithm; MA: memetic algorithm; FCFS: First Come First Served. Note: The significance of the values given in boldface is reflected in the experimental analysis paragraph of case 3.

Table 9. The optimization results. Item

Best sequences

Case 1 Case 2 Case 3

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

fed by AGV can be obtained by applying the method of IHSMTSA, as shown in Table 9. From the three cases mentioned above, it is obvious that the result obtained by IHSMTSA is much better than ‘‘FCFS’’ in real-world production. And the proposed IHSMTSA has better performance than other algorithms compared in terms of both the quality of solutions and the rate of convergence, when solving the AGVSP described in this article. Especially, the advantages of IHSMTSA will stand out when dealing with larger scale AGVSP such as cases 2 and 3. Moreover, the performance of IHSMTSA is much better than other algorithms and FCFS in terms of the three indicators, that is, the standard deviation of the waiting time of CNC material buffers, the average waiting time of CNC material buffers, and the total travel distance of AGV, which means that IHSMTSA has a better performance than other algorithms and FCFS compared to increasing the overall capacity of CNC machines, as well as improving the energy efficiency of AGVs. Therefore, the IHSMTSA algorithm can greatly improve the profitability of enterprises.

heuristic method based on HS algorithm, that is, IHSMTSA, is applied to solve the AGVSP. The current design of IHSMTSA incorporates four major improvements to the traditional HS algorithm: (1) Considering the characteristics of the problem of interest here, the discrete coding method is adopted which can greatly reduce the search range in the solution space. (2) The OBL approach is adopted here to initialize the HM, which can expand the initial solution space in the searching area, thereby enhancing the probability of landing on the global optimum. (3) The parameter HMCR is dynamically changed, which can prevent IHSMTSA from falling into the local optima in early iterations and can enhance the local search ability in later iterations. (4) A local search strategy is recommended to improve the quality of the optimal solution. As demonstrated in the simulation study, our IHSMTSA is notably superior to other algorithms in terms of the performance in solving AGVSP. The contributions of this article are summarized as follows: 1.

Conclusion In this article, we first proposed an effective objective function model for AGVSP, which incorporates two indicators, that is, the standard deviation of the waiting time of CNC material buffers and the total travel distance of AGV, reflecting the overall capacity of CNC machines and the energy efficiency of AGVs during the actual production, respectively. Then, a novel meta-

2.

An improved HS algorithm, IHSMTSA, is developed to address the AGVSP, which incorporates a multitude of major improvements to the traditional HS algorithm. The standard deviation of the waiting time of CNC material buffers and the total travel distance of AGV are reduced greatly. The workshop capacity is also increased. As a result, the production efficiency is increased and the energy consumption is reduced to a great extent.

Li et al. On a going-forward basis, attempts will be made to implement the IHSMTSA in a real-world intelligent production workshop so as to further assess its effectiveness. Then, this algorithm can be applied to more production workshop. On this basis, we can continue to study and improve the algorithm so that the algorithm can be applied to more areas.

13

12.

13.

14.

Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

15.

Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is supported by National Natural Science Foundation of China (NSFC) under grant no. 51421062.

16.

References

18.

1. Miyamoto T and Inoue K. Local and random searches for dispatch and conflict-free routing problem of capacitated AGV systems. Comput Ind Eng 2016; 91: 1–9. 2. Qiu L, Hsu W-J, Huang SY, et al. Scheduling and routing algorithms for AGVs: a survey. Int J Prod Res 2002; 40: 745–760. 3. Clarke LW, Hane CA, Johnson EL, et al. Maintenance and crew considerations in fleet assignment. Transp Sci 1996; 30: 249–260. 4. French S. Sequencing and scheduling: an introduction to the mathematics of the job-shop. Hoboken, NJ: John Wiley & Sons, 1982. 5. Soylu M, Ozdemirel NE and Kayaligil S. A selforganizing neural network approach for the single AGV routing problem. Eur J Oper Res 2000; 121: 124–137. 6. Tsirimpas P, Tatarakis A, Minis I, et al. Single vehicle routing with a predefined customer sequence and multiple depot returns. Eur J Oper Res 2008; 187: 483–495. 7. Minis I and Tatarakis A. Stochastic single vehicle routing problem with delivery and pick up and a predefined customer sequence. Eur J Oper Res 2011; 213: 37–51. 8. Hu WH, Mao JF and Wei KJ. Energy-efficient rail guided vehicle routing for two-sided loading/unloading automated freight handling system. Eur J Oper Res 2016; 258: 1–30. 9. Qudeiri JEA, Umer U, Abidi MH, et al. On multistage approach for flexible routing in flexible manufacturing systems. Adv Mech Eng 2016; 8: 1–11. 10. Gultekin H, Karasan OE and Akturk MS. Pure cycles in flexible robotic cells. Comput Oper Res 2009; 36: 329–343. 11. Nejad MG, Gu¨den H, Vizva´ri B, et al. A mathematical model and simulated annealing algorithm for solving the

17.

19.

20.

21.

22.

23.

24.

25.

cyclic scheduling problem of a flexible robotic cell. Adv Mech Eng 2018; 10: 1–12. Geem ZW, Kim JH and Loganathan G. A new heuristic optimization algorithm: harmony search. Simul 2001; 76: 60–68. Zhang D, Zhang C, Xu F, et al. A hybrid genetic algorithm used in vehicle dispatching for JIT distribution in NC workshop. IFAC 2015; 48: 898–903. Maroosi A, Muniyandi RC, Sundararajan E, et al. A parallel membrane inspired harmony search for optimization problems: a case study based on a flexible job shop scheduling problem. Appl Soft Comput 2016; 49: 120–136. Guo Z, Shi L, Chen L, et al. A harmony search-based memetic optimization model for integrated production and transportation scheduling in MTO manufacturing. Omega 2017; 66: 327–343. Zeng B and Dong Y. An improved harmony search based energy-efficient routing algorithm for wireless sensor networks. Appl Soft Comput 2016; 41: 135–147. Abedinpourshotorban H, Hasan S, Shamsuddin SM, et al. A differential-based harmony search algorithm for the optimization of continuous problems. Expert Syst Appl 2016; 62: 317–332. Alfailakawi MG, Ahmad I and Hamdan S. Harmonysearch algorithm for 2D nearest neighbor quantum circuits realization. Expert Syst Appl 2016; 61: 16–27. Cheng MY, Prayogo D, Wu YW, et al. A hybrid harmony search algorithm for discrete sizing optimization of truss structure. Automat Constr 2016; 69: 21–33. Papa JP, Scheirer W and Cox DD. Fine-tuning deep belief networks using harmony search. Appl Soft Comput 2016; 46: 875–885. Zeng B, Dong Y, Li X, et al. IHSCR: energy-efficient clustering and routing for wireless sensor networks based on harmony search algorithm. Int J Distrib Sens N 2017; 13: 1550147717741103. Ouaddah A and Boughaci D. Harmony search algorithm for image reconstruction from projections. Appl Soft Comput 2016; 46: 924–935. Kamboj VK, Bath S and Dhillon J. Implementation of hybrid harmony search/random search algorithm for single area unit commitment problem. Int J Elec Power 2016; 77: 228–249. Lee KS and Geem ZW. A new meta-heuristic algorithm for continuous engineering optimization: harmony search theory and practice. Comput Method Appl M 2005; 194: 3902–3933. Tizhoosh HR. Opposition-based learning: a new scheme for machine intelligence. In: Proceedings of the 2005 international conference on computational intelligence for modelling, control and automation, and international conference on intelligent agents, web technologies and internet commerce, Vienna, 28–30 November 2005, pp.695–701. New York: IEEE.