Abstract

Automatic production system scheduling problem under a just-in-time environment is researched in this paper. The automatic production system is composed of many tanks and one robotic, the tank of the researched problem is responsible for processing the job, and the robotic moves the job from one tank to the other tank. The difference between the researched problem and the classic shop scheduling problem is that the former must consider job scheduling and the robotic move sequence, but the latter considers only job scheduling. For optimizing simultaneously job scheduling and robotic move sequence in the proposed problem and minimizing total earliness/tardiness, an improved NEH (Nawaz-Enscore-Ham) and variable search (INEH-VNS) algorithm are developed. In the proposed method, firstly, to obtain initial solution, an improved NEH is shown. Secondly, for computing value of the objective function, the double procedure method is constructed. Thirdly, according to the properties of the proposed problem, three neighborhood structures, adjacent exchange, random insertion, and job exchange, are investigated. To test the performance of the INEH-VNS, 100 instances are randomly generated. When the run time is the same, compared with CPLEX 12.5, the INEH-VNS algorithm can find high-quality approximate optimal solution, a special big scale. Compared with the G-VNS algorithm, the average improvement rate of the approximate optimal solution is 45.9%, and the average stability rate of the INEH-VNS algorithm enhances 75.04%. That is to say, the INEH-VNS algorithm is outstanding and more effective.

1. Introduction

The competition between enterprises in the market is becoming increasingly intense. Therefore, the concept of just-in-time (JIT), which is paid attention to by many industries and has been widely adopted over the years to improve production efficiency, stems from the production field [1, 2]. One of the main concepts in JIT production systems is that the makespan is equal to the due date, avoiding unnecessary inventory and late delivery; that is to say, the makespan is not more than the due date, and the makespan is not less than the due date [3]. If the makespan is less than the due date, many questions, e.g., inventory cost, product damage, and economic depreciation, will be generated. If the makespan is more than the due date, a variety of adverse consequences, e.g., customer loss, contract default, and commercial reputation damage, will arise [4].

Automatic production system which is applied to printed circuit board (PCB) industry, wafer fabrication, automobile manufacturing industry, and steel industry [59] is an advanced intelligent manufacturing system. The products which are manufactured by the automatic production system are easy to depreciate, so JIT production mode is adopted. In published literatures, the classical shop scheduling problem considering due date is not researched in robotic move sequence [1013]. The automatic production system scheduling problem is studied by robotic move sequence but is not considered the due date [59]. To the best of our knowledge, the scheduling problem considering JIT environment and robotic move sequence is studied seldom. Due to the significance of JIT production in industrial environment to overcome variations from the demand providing customer satisfaction and widely applied of automatic production system in industrial environment, in this research, an improved NEH and variable neighborhood search (INEH-VNS) that optimizes the proposed problem is developed for improving production efficiency of the automatic production system and alleviating the contradiction between supply and demand.

In this paper, the automatic production system scheduling problem considering due date is investigated. The proposed automatic production system is composed of a number of tanks, including one input station, one output station, and one robotic which is controlled by computer. All jobs which will be processed are loaded into the input station and are processed in each tank in order; finally, all jobs which have been finished are loaded into the output station. Each job is moved from one tank to the other tank. The due date of each job is given.

To satisfy the due date of the job as much as possible, the INEH-VNS is designed for optimizing total earliness/tardiness of more than one tank’s automatic production system scheduling problem. Our contributions are listed in the following: (1)According to the properties of the proposed problem, the improved NEH is developed for generating initial solution(2)To compute the value of the objective function, the double course method is proposed(3)According to the solution properties of the proposed problem, three neighborhood constructs are shown for designing variable neighborhood search (VNS).(4)In the previous work, the automatic production system scheduling problem with one tank under JIT environment is studied. In this paper, the INEH-VNS is developed for solving the automatic production system scheduling problem with more than one tank under JIT environment and optimizing simultaneously the job scheduling and the robotic move sequence

The remainder of this paper is organized as follows. Section 2 describes the literature on scheduling problems minimizing total earliness/tardiness. Section 3 presents the INEH-VNS. Section 4 exhibits computational results. Finally, the conclusions and the challenge work are given in Section 5.

2. Literature Review

The shop scheduling under JIT environment means the due date and the makespan are as consistent as possible, so the objective of the shop scheduling problem under JIT environment is to minimize total earliness/tardiness. In this section, the literatures which the objective is to optimize total earliness/tardiness are described.

For the single-machine scheduling problem under JIT environment, the mixed integer programming model is developed by Woodruff and Spearman when the due date is given [14]. Coleman researched the problem with job scheduling sequence-dependent setup times and investigated the mixed integer programming model [15]. To obtain the total weighted earliness/tardiness of job scheduling sequence-dependent switching cost, the nonlinear mathematic programming model is constructed [16]. For the problem with quadratic earliness/tardiness penalties, the mixed integer programming model with uncertainty due date is shown [17]. The exact and heuristic methods are proposed [18] and improved the result of literature [17]. According to literature [17], the mathematic programming with considering learning effect and release dates is exhibited [19]. Kayvanfar et al. researched the single machine scheduling problem with parallel machine under JIT environment and minimized total earliness/tardiness [11]. The single machine scheduling problem considering general due date is studied; if the objective function is to minimize the number of tardiness jobs, the problem is polynomially solvable; if the objective function is to optimize the number of weighted tardiness jobs, the problem is NP hard; if the problem is extended to allow job rejection, the problem is NP hard [20]. For multiscenario scheduling problem, to obtain the number of weighted JIT jobs, the problem is solvable in polynomial time [21]. Otherwise, if the objective function is a total weighted makespan, the problem is solvable in (is the number of jobs) [22]. When the due dates are assigned to each job depending on its order, and the objective is to minimize the total penalty for the earliness and tardiness of each job, the problem is NP hard [16]. For solving the single machine scheduling problem under JIT environment, the decomposition algorithm is proposed, so the upper bound and the lower bound are obtained, and the optimal solution is found by the proposed branch and bound; for optimizing total weighted tardiness, four heuristic methods are designed [23]. Mosheiov et al. introduced pseudopolynomial dynamic programming algorithms for minimizing the number of tardy jobs when the jobs allowed rejection [20]. Alidaee reviewed the published mathematic programming and then introduced future work [12].

For the two machines scheduling problem under JIT environment, the problem with simultaneous consideration of common due date assignment, convex resource allocation, and learning effect is solvable in polynomial time for minimizing the total resource consumption cost [24]. Shi and Wang researched the problem with common due window assignment, learning effect, and resource allocation, and they proofed the problem is solvable in polynomial time for optimizing the linear weighted sum of job earliness, tardiness, due window size, and resource cost [25]. The results are the same if the slack due date substitutes for common due date [26, 27]. Based on the results of literature [25, 27], for minimizing the scheduling cost and the resource consumption cost, the scheduling problem with different due window assignment and learning effects is investigated, and the polynomial algorithm is proposed [28]. The classical scheduling problem models considering due date assignment are extended and solvable in polynomial time [29].

For more than two machines scheduling problem under JIT environment, Jiang et al. addresses the problem considering position-dependent weights and discusses the common due date assignment model and slack due date assignment model. The polynomial time algorithm is developed for minimizing total cost [30]. Lv et al. proposed polynomial time algorithms for the problem in literature, the computational time complexity is reduced from O(n2logn) to O(nlogn) [31]. For the problem with two distinct due dates, Koulamas et al. show that the proportionate flow shop problem with variable machine speeds remains NP-hard; they then show that the no-wait problem with ordered jobs and a maximal last machine is solvable in time; They also show that the corresponding problem with the minimum number of tardy jobs is solvable in time [32].

In summary, for the flow shop scheduling problem under JIT environment, mathematics model, optimization algorithm, and computational complexity are studied, but the robot moving job between the operations is not considered. Thus, the results of the literatures [1432] cannot be applied to the proposed problem in this paper.

For the automatic production system scheduling problem, Zhao and Guo [9] and Zhao and Guo [33] constructed optimization algorithm, respectively, for optimizing the production cycle. Wu et al. proposed a differential evolution algorithm for solving robotic cell scheduling problem with batch-processing machines [34]. Majumder et al. researched robotic cell scheduling problem with sequence-dependent setup times and addressed a bacterial foraging optimization algorithm [35]. Elmi et al. showed an integer programming for cyclic flow shop robotic cell scheduling problem with multiple part types [36]. The objective of the literatures [9, 3336] is to minimize the production cycle. Wang et al. studied the automatic production system scheduling problem under JIT environment, developed a mixed integer programming, and proposed a hybrid guided tabu search algorithm to optimize the makespan and total weighted tardiness [37]. The difference between the proposed problem and the problem of literature [37] is that the parallel-tank scheduling problem with hoist and group constraints is researched in literature [37], but more than one tank flow shop scheduling problem with hoist is investigated in this paper, so the method of the literature [37] cannot be applied in this paper. The method for solving the proposed problem is worth studying.

Since the automatic production system scheduling problem is NP-hard, all kinds of algorithms, e.g., chemical reaction optimization [9], branch and bound algorithm [33], differential evolution algorithm [34], bacterial foraging optimization algorithm [35], and hybrid guided tabu search algorithm [37], are proposed for obtaining the optimal solution. There are two shortcomings in these algorithms; one is it is difficult to find a solution for big scale problem, and the other is it is troublesome to find parameter values for giving problem. Variable neighborhood search that is proposed by Mladenović and Hansen [38] does not involve parameters and can find optimal solution of the problem. Variable neighborhood search is applied to solve various combinatorial optimization problem [3942]. In current, it is seldom found that variable neighborhood search is used to solve the automatic production system scheduling problem under just-in-time environment. In addition, for enhancing the performance of the proposed algorithm, the improved NEH that is developed generates initial solution. In this paper, an improved NEH and variable neighborhood search, namely, INEH-VNS, are investigated for minimizing total earliness/tardiness and finding the optimal job scheduling and robotic move sequence.

3. Problem Description

The automatic production system is composed of one input station , tanks , ,,, one output station , and single robot in this paper. The input station contains all jobs which will be processed, the output station holds all jobs which have been finished, and all tanks , , , process jobs. The robot moves the job from one tank to the other tank. jobs are processed by the automatic production system at the same time. Firstly, each job is unloaded from the input station . Then, the job is successively processed in tank , , , . Finally, the job is loaded into the output station , that means the job is finished. The processing time of each job is not all same in tank (). So, the proposed problem in this paper is a flow shop scheduling problem with single robot. It is difficult that job scheduling and robotic move sequence are simultaneously optimized in current researched problem. At the beginning of the period, the input station contains all unprocessed jobs, and all tanks and the output station are empty, which is different from published literatures automatic production system scheduling problem, because the tanks of published literatures automatic production system scheduling problem are not all empty. At any time, the robot can hold only one job at most, and the tank can contain only one job at most. Furthermore, processing interruptions and preemption are not considered. In this paper, three operations: (i) the robot unloads job from tank ; (ii) the robot carries the job from tank to tank ; (iii) the robot loads the job into tank (), are called robotic activity. The robot moves from tank to tank () is named void move. For satisfying the due date of job as much as possible, the objective of total earliness/tardiness is proposed in this paper. To optimize the objective value, three types of constraint, processing time constraints, robot-containing constraints, and tank capacity constraints, must be considered.

There is an obvious difference between the proposed problem and the classic shop scheduling problem in that job scheduling and the robotic move sequence of the proposed problem are simultaneously optimized. Since the flow shop scheduling problem for two machines is NP-hard to obtain total earliness/tardiness [43], the proposed problem is also NP-hard. The notations that are used in this paper are defined as follows:

Decision variables:

: the robotic activity start time (; ).

Intermediate variables:

: binary variables. When is true, , otherwise (; ).

: total earliness\tardiness time.

: earliness\tardiness time of the (). job.

: makespan of the (). job.

Parameters or notations:

: the processing time lower bound of the () job in tank ().

: the due date of the () job.

: the time of robotic activity(;).

: the time of void move().

Other notations:

: , the jobs permutation.

: , the robotic activities permutation.

: the robotic activity in feasible solution().

: the robotic activity corresponds to the subscript of the tank in feasible solution().

: the robotic activity corresponds to the job in feasible solution().

: the start time of robotic activity and ().

The mixed integer programming is shown.

The objective function that is to minimize total earliness\tardiness is shown by equality (1). Inequality (2) is a processing time constraint. Inequality (3) and Inequality (4) are robots containing constraints. Equality (5) is the tank capacity constraint. Makespan is computed by equality (6). Inequality (7) and Inequality (8) are triangle inequality constraints. Inequality (9) and Inequality (10) are nonnegative constraints and binary variables, respectively.

4. INEH-VNS

Variable neighborhood search that improves local search algorithm is composed of initial solution, neighborhood structure, and stop condition. Next initial solution and neighborhood structure of the proposed algorithm that are designed are introduced.

4.1. INEH

Since the job scheduling and the robot move sequence of the proposed problem are optimized, if the job scheduling is discussed first of all and, secondly, the robot move sequence is optimized, it is hard to find the optimal solution of the proposed problem, vice versa. In addition, since the job scheduling and the robot move sequence are not one-to-one, it is not easy to find the optimal solution to the current researched problem if the job scheduling or the robot move sequence represents the initial solution. In this paper, to obtain initial solution, the robotic activity sequence [12] that is composed of the job scheduling and the robot move sequence is applied. It has been proven that the NEH heuristic is an effective heuristic\method to address flow-shop scheduling problems [44]. The improved NEH is developed for constructing the initial solution, and the detailed process is as follows:

Step 1. Since the time which the job moves from tank to tank is considered. In order to apply NEH method, the processing time of the job, called , is revised. The detailed process is shown in equation (3).

Step 2. Compute the sum of the processing time of the () job in all tanks, named . It means equation (4) holds

Step 3. All () are arranged in descending order; then, the job sequence set is obtained.

Step 4. The first two robotic activities of the first job of set are inserted in the first two places of the initial solution, namely, . The last two robotic activities of the last job of set are inserted in the last two places of initial solution .

Step 5. The first robotic activity of each job, except the first job and the last job, in set is inserted in the place of initial solution in turn.

Step 6. According to property 1 (the appendix), the remainder robotic activities of each job in set are inserted in the initial solution ; then, the initial solution is obtained.

Example 1. The number of jobs is four, and the number of tanks, including the input station and output station, is five. The processing time of jobs is listed in Table 1. The robotic activities are shown in Table 2. The time of the void move is equal to (). The time of robotic activity is 6 (,). The due date of jobs is .

According to Step 1, Step 2, and Step 3, the set . Figures 1(a) and 1(b) are obtained according to Step 4 and Step 5, respectively. Figures 1(c) and 1(b) are obtained according to Step 6. Finally, the initial solution , Figure 1(d), is given.

4.2. Double Procedure Method

In order to calculate the makespan (), the double procedure method is developed. Firstly, the relation between the makespan () and the start time of robotic activity is built. Secondly, the start time of robotic activity is computed that is separated into two procedures: one is the complete time of job on the tank, and the other is the time of void move. If the complete time of job on the tank is more than the time of the void move, the start time of the robotic activity is equal to the complete time of job on the tank, vice versa. If is true, the makespan holds, where the is obtained according equation (5), equation (14), and equation (15).

Example 2. The relation between the robotic activities in feasible solution is exhibited in Figure 2. The start time of each robotic activity is computed as follows. So, the makespan is equal to 62.
, , , , , , .

After the double procedure method is applied, , , and are true.

4.3. Neighborhood Structure

The feasible solution properties of the proposed problem are applied, and three neighborhood structures, adjacent exchange, random insertion, and job exchange, are shown.

4.3.1. Adjacent Exchange

The represents a feasible solution of the current researched problem, and and are two robotic activities. If the inequality holds, the robotic activities, and , are exchanged. New solution is obtained and is feasible. For improving convergence, when the robotic activities are exchanged, the job of smaller processing time is moved from right to left. According to the following steps, the neighborhood structure of adjacent exchange is realized.

Step 1. The place, called , is randomly selected in feasible solution. There is a robotic activity at place .

Step 2. If the inequality holds, Step 3 is implemented. Otherwise, Step 4 is implemented.

Step 3. If the inequality is true, the robotic activities, and , are exchanged. Otherwise, Step 4 is implemented.

Step 4. If the inequality holds, Step 5 is implemented. Otherwise, Step 1 is implemented.

Step 5. If the inequality is true, the robotic activities, and , are exchanged. Otherwise, Step 1 is implemented.

Example 3. A feasible solution is shown in Figure 3(a). The place, , is selected, and there are two robotic activities , at place and , respectively. The inequality is true, and the inequality holds. After the robotic activities and are exchanged, new solution which is shown in Figure 3(b) is obtained and is feasible.

4.3.2. Random Insertion

For exploring field, the neighborhood of random insertion is developed. Since there is only one tank between successive two tanks (), or there is only one tank between successive two tanks (), the neighborhood structure of random insertion which is constructed is listed as follows:

Step 1. The place, denoted , is randomly selected in feasible solution.

Step 2. The tanks, and , which are located in the right-most at the left of the place are found. The tanks and correspond to the places and , respectively.

Step 3. The tanks, and , which are located in the left-most at the right of the place are found. The tanks and correspond to the places and , respectively.

Step 4. Let and .

Step 5. If is true, go to Step 1. Otherwise, Step 6 is implemented.

Step 6. A random number is obtained between and , and holds.

Step 7. The robotic activity is inserted at the place . The relative position of the remainder robot activities is unchanged.

Example 4. The value of is 9, and . The result is shown in Figure 4(a). According to Step 2, , , and , are shown in Figure 4(b). According to Step 3, , , , and are shown in Figure 4(c). Figure 4(d) shows the results of and , and is true. Finally, the place is selected that is shown in Figure 4(e). Figure 4(f) shows that the robotic activity is inserted at the place .

4.3.3. Job Exchange

Two neighborhood structures, neighborhood exchange and random insertion, change the robotic move sequence, but the job scheduling does not change. Therefore, to find a much better solution, the neighborhood of job exchange is investigated for changing job sequence. After the robotic activities of the job and the robotic activities of the job have corresponded to swap in feasible solution, new solution is obtained and is feasible. To enhance algorithm efficiency, the job that the due date is smaller is prior processed. The following steps show the process that job exchange is constructed.

Step 1. The two different jobs, and , are randomly selected. Without loss of generality, the job is processed before the job .

Step 2. The due dates of job and are and , respectively.

Step 3. If holds, implement Step 4; otherwise, implement Step 1.

Step 4. The robotic activities of the job and the robotic activities of the job correspond to swap.

Example 5. Figure 5(a) shows the randomly selected two jobs, (the shadow is the slash) and (the shadow is the square), respectively. The due dates of job and are and , respectively. Since is true, the robotic activities of the job and the robotic activities of the job corresponded to swap. New feasible solution is obtained and is shown in Figure 5(b).

As mentioned above, INEH-VNS flowchart is shown in Figure 6. The detailed procedure of INEH-VNS is illustrated as algorithm in Algorithm 1.

Algorithm INEH-VNS
Input: processing time of job, the time of robotic activity, the time of void move, the number of jobs , the number of tanks , the due date of job, the time of void move, neighborhood structures, (Adjacent exchange), (Random insertion), and (Job exchange). Let .
Output: optimal solution , and optimal objective value .
//First stage: generate initiation solution by INEH.
For=0 to
  For =0 to
    Record the amend processing time of job as
  End for
End for
For=0 to
  For =0 to
    Calculate the sum of processing time of each job in all tanks as
  End for
End for
For=0 to
  For =0 to
    Descending order , obtain the job sequence set .
  End for
End for
For=0 to
  For = to
    According to property 1, initiation solution is obtained.
  End for
End for
The double procedure method is proposed, the objective value is computed. Let and .
//Second stage: find optimal solution by VNS.
Repeat
  Forto 3
    Find the best neighbor of in ;
    When the solution is , the objective value is ;
    If , then and ;
    Otherwise ;
  End for
Until Stopping criteria.

5. Computational Result

To test the performance of the proposed method, the computational experiments are shown in this section. Extensive experiments are conducted on a set of problems. All experiments are implemented using Microsoft Visual Studio 2010 and are run on an Intel(R) Core(TM)i5-4344 [email protected] GHz and RAM 8.0 GB PC.

In order to evaluate the performance of the proposed method, we develop the other algorithms, namely, G-VNS and simulated annealing (SA), respectively. There are two differences between the G-VNS and the proposed method. One is all the due dates () are arranged in descending order, and then the job sequence set is obtained for generating the initial solution of the G-VNS; the other is the which is replaced by in Step 3 of Section 4.3.3, where is computed by equation (4). The initial solution of SA is generated by the INEH method, and the proposed neighborhood structures in this study are applied in the SA.

For single machine scheduling problem under JIT environment, the due dates obey the uniform distribution of , where is the lower bound, is the tardiness factor, and is the due date range [45]. The computational method of the lower bound is developed for machine scheduling problem under JIT environment by Taillard [46]. The computational method which is proposed by Taillard does not consider the time of robotic activity. In this paper, to obtain a lower bound , equation (3) is used. Then, the tardiness factor and the due date range are 0.2 and 0.6, respectively [47].

Since bench marks are not found for the proposed problem, one hundred instances are generated according to literature [48]. is an integer and . The robotic move time for all ;. The void move time for all . The number of tanks is an even number between 2 and 20. The number of jobs varies between 5 and 50 and is an integer multiple of 5.

On the test of the instances, the proposed algorithm, called INEH-VNS, compares with the G-VNS, CPLEX12.5, and SA. Termination condition is the maximum computational time does not exceed 600 seconds. The other parameters of CPLEX 12.5 are default. Each instance is run 10 times, and the average of each instance, namely, (), is recorded. The performance measure is defined by equation (8), equation (9), and equation (10).

5.1. Comparison Approximate Optimal Solution

The proposed algorithm compares with CPLEX 12.5, the G-VNS, and the SA, and the results are displayed in Tables 37. The results indicate high-quality solutions that can be found by the INEH-VNS.

The INEH-VNS algorithm compares with CPLEX 12.5, and the INEH-VNS algorithm is more effective for large-scale problems. The approximate optimal solutions of all instances are found by the INEH-VNS. Otherwise, 27% of the instances of the approximate optimal solutions are solved by CPLEX 12.5. Furthermore, in the instances in which the approximate optimal solutions are discovered by CPLEX 12.5, the approximate optimal solutions of eleven instances that are solved by INEH-VNS are better than the optimal solutions obtained by CPLEX 12.5. The best improvement rate is 77.96%, the worst improvement rate is 0.44%, and the average of is 34.55%. The best improvement rate and the worst improvement rate are shown in bold.

The INEH-VNS algorithm compares with the SA algorithm, and the INEH-VNS algorithm is more effective than the G-VNS algorithm. That is to say, if the run time is the same, better approximate solutions are found by the INEH-VNS algorithm in all instances. The best improvement rate is 77.78%, the worst improvement rate is 7.09%, and the average of is 54.18%. The best improvement rate and the worst improvement rate are exhibited by bold and shadow. The G-VNS algorithm compares with the SA algorithm, The G-VNS algorithm is better than the SA algorithm. 85 out of 100 instances of the solution quality which is found by the G-VNS are better than the solution quality which is solved by the SA algorithm. We find that if the number of tank is more than 4, the G-VNS algorithm is more competitive than the SA algorithm. The INEH-VNS algorithm and the G-VNS algorithm are more effective than the SA algorithm.

The INEH-VNS algorithm compares with the G-VNS algorithm, and the INEH-VNS algorithm is more competitive than the G-VNS algorithm. The approximate optimal solutions of all instances are obtained by the INEH-VNS algorithm and the G-VNS algorithm, but 94 out of 100 instances of the solution quality which is found by the INEH-VNS are better than the solution quality which is solved by the G-VNS. The biggest is 69.86%, the smallest is 0.59%, and the average of is 45.9%. The biggest and the smallest are shown by shadow. The results show job exchange that is designed according to descending order due date is more effective.

5.2. Stability Comparison

In this subsection, the stability of the INEH-VNS algorithm, the G-VNS algorithm, and the SA algorithm is discussed. For evaluating algorithm stability, the standard deviation is used. Since CPLEX 12.5 can find 27% of the instances of the approximate optimal solutions, CPLEX 12.5 is not considered in this subsection. The results are shown in Tables 711. The stability rate, called and , is computed according to equation (17) and equation (18), where STSA, STG-VNS, and STINEH-VNS represent the standard deviation of the SA algorithm, the G-VNS algorithm, and the INEH-VNS algorithm, respectively.

In Tables 812, the stability of the INEH-VNS algorithm and the G-VNS algorithm is much better than the stability of the SA algorithm. The stability of the INEH-VNS algorithm is much better than the stability of the G-VNS algorithm except for six instances. In the remainder instances, the standard deviation of the G-VNS algorithm is more than the standard deviation of the INEH-VNS. When and or , the same solution is obtained by the INEH-VNS algorithm; thus, the stability of the INEH-VNS is the best. The results are shown in bold. The standard deviation of the INEH-VNS algorithm, compared with the standard deviation of the G-VNS algorithm, is the worst improvement.

The rate is 3.46%, and the result is shown by shadow. That is to say, when and , there is little difference for the stability of the INEH-VNS algorithm and the G-VNS algorithm. In these instances where the standard deviation of the INEH-VNS algorithm is more than the standard deviation of the G-VNS, the average of is 75.04%. These results show that job exchange that is designed according to descending order due date is more stable; that is to say, the INEH-VNS is better convergent.

5.3. Convergence

Because CPLEX 12.5 can find 27% of the instances of the approximate optimal solutions and the performance of the SA algorithm is poor, in this subsection, the convergence of the INEH-VNS algorithm and the G-VNS algorithm is described. When the number of tanks is 50, and the number of jobs is 10, 12, 14, 16, 18, and 20, respectively, the convergence of the INEH-VNS algorithm and the G-VNS algorithm is shown in Figures 7(a)7(f). In Figure 7, the INEH-VNS algorithm and the G-VNS algorithm cause a quick decrease in the objective function value at the beginning and continue reducing this function with lower speed until it converges to its best solution with a reduced speed of approximately zero. But the best solution which is found by the INEH-VNS algorithm is much better than the best solution which is found by the G-VNS algorithm. That is to say, the proposed algorithm is effective.

6. Conclusions

In this article, the automatic production system scheduling problem under just-in-time environment is investigated, and the INEH-VNS algorithm is developed for minimizing total earliness/tardiness time.

In the proposed algorithm, the improved NEH robotic activity method is investigated to obtain initial solution. The double procedure method is designed to compute the value of the objective function. According to the property of the current researched problem, three neighborhood structures, adjacent exchange, random insertion, and job exchange, are discussed for finding approximate optimal solution. Compared with CPLEX 12.5, the G-VNS algorithm, and the SA algorithm based on solving randomly generating instances, the proposed algorithm is outstanding. Next, we will research multiobjective automatic production system scheduling problems under just-in-time environment.

Appendix

Property 1. If is a solution of the proposed problem, , , , and , , all hold.

Proof. Because at the beginning of the period, the input device is the full exclusion of all tanks and the output device. The job is moved to the tank ; that is to say, the tank is occupied. According to condition (iii) of definition 1, the other job cannot be loaded into the tank . The robot stays in the tank . After the job is processed completely on the tank , the robot moved the job from tank to tank . So, , , and are all true. All jobs are done, the output device is full, and all tanks and the input device are empty. The last two robotic activities of the solution must implement the job from tank to tank , and the robot stays in the tank ; when the job is done on the tank , the robot moves the job from tank to tank . That is to say, , , and all hold.

Data Availability

The data used to support the findings of this study are included within the article.

Disclosure

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

L.Q. and Z.X. were responsible for the conceptualization. H.Y. was responsible for the methodology. L.Q. was responsible for the software. Z.X. and H.Y. were responsible for the validation. Z.X. was responsible for the formal analysis. L.Q. was responsible for the investigation. H.Y. and Y. S. were responsible for the resources. H.Y. and Y. S. were responsible for the data curation. L.Q. and Y. S. were responsible for the writing—original draft preparation. Z.X. and H.Y. were responsible for the writing—review and editing. Z.X. and Y. S. were responsible for the visualization. Z.X. was responsible for the supervision. Z.X. and L.Q. were responsible for the project administration. Z. X and L.Q. were responsible for the funding acquisition. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This research was funded by the Humanities and Social Sciences Project of the Ministry of Education of the People’s Republic of China under Grant No. 18YJCZH26. It was also supported by the Social Science Planning PhD Project of Chongqing under Grant No. 2019BS072, and the Science and Technology Research Program of Chongqing Municipal Education Commission under Grant Nos. KJQN201901337, KJQN202001316, and KJQN202101324.