Research Article | Open Access

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

Ik Sun Lee, "A Scheduling Problem to Minimize Total Weighted Tardiness in the Two-Stage Assembly Flowshop", Mathematical Problems in Engineering, vol. 2020, Article ID 9723439, 10 pages, 2020. https://doi.org/10.1155/2020/9723439

# A Scheduling Problem to Minimize Total Weighted Tardiness in the Two-Stage Assembly Flowshop

Revised15 Jun 2020
Accepted10 Aug 2020
Published27 Sep 2020

#### Abstract

This paper considers an assembly scheduling problem which minimizes the total weighted tardiness. The system consists of multiple manufacturing machines in the first stage and an assembly machine in the second stage. A job is composed of multiple components, each component is manufactured on a dedicated machine in the first stage, and the assembly starts on the assembly machine after the completion of all the components’ manufacturing in the first stage. This paper is to try to find effective and efficient schedules to minimize the total weighted tardiness of jobs. This paper analyzes some solution properties to improve search effectiveness and efficiency. Moreover, some lower-bound procedures are derived, and then four constructive heuristics and a branch-and-bound algorithm are suggested. The performances of the derived algorithms are evaluated through numerical experiments. For the problems with no more than 15 jobs, the branch-and-bound algorithm finds the optimal solutions within 20 minutes, and the performance of the four heuristics seems to be similar. However, for problems with more than 30 jobs, HWT shows the best performance among the four heuristics except for two machines.

#### 1. Introduction

This study considers a two-stage assembly production system, which manufactures components in the first stage and then assembles all the components in the second stage. The assembly can begin after all the components belonging to the product are completed. In the system, two stages are physically located in a series as the components’ production stage and the assembly stage, but the components are produced on the manufacturing machines in parallel.

Framinan et al.  and Pan et al.  conducted a review paper on the assembly scheduling problem. From the review papers, there are two types of research depending on whether dedicated machines or flexible machines are used for manufacturing components in the first stage. First, the dedicated machine environment represents that each component is processed on a designated machine. Second, the flexible machine environment means that each machine can process all types of components.

This paper considers the first situation with dedicated machines in the first stage. Lee et al.  introduced a two-stage assembly scheduling problem with two dedicated manufacturing machines in the first stage and one assembly machine in the second stage. For the makespan minimization, Lee et al.  proved that the associated problem is NP-complete, proposed three heuristics based on Johnson's rule, and derived their error bounds. Potts et al.  conducted an extension of Lee et al.  by considering m machines in the first stage and proposed some heuristics using Johnson's rule. Hariri and Potts  proposed an effective branch-and-bound algorithm to find optimal solutions. Sun et al.  proposed 14 heuristics for the problem of Potts et al. . Tozkapan et al.  and Lee  considered the total weighted completion minimization in a two-stage assembly scheduling problem and suggested a branch-and-bound algorithm. Framinan and Perez-Gonzalez  proposed a metaheuristic algorithm for the total completion minimization in a two-stage assembly scheduling problem. Talens et al.  addressed the two-stage multimachine assembly scheduling problem under a situation with several dedicated parallel machines in the first stage and identical parallel machines in the second stage. They proposed two heuristic algorithms for the total completion time minimization.

Allahverdi and Al-Anzi  considered the maximum lateness objective in the two-stage assembly scheduling problem and derived a metaheuristic based on particle swarm optimization and tabu search. Allahverdi and Aydilek , Ha et al. , and Lee and Bang  considered the two-stage assembly scheduling problem to minimize the total tardiness with two dedicated machines in the first stage and one assembly machine in the second stage. Allahverdi and Aydilek  proposed metaheuristic algorithms, and Ha et al.  and Lee and Bang  derived a branch-and-bound algorithm. Kazemi et al.  investigated a two-stage assembly flowshop scheduling problem with a batched delivery system where there are m independent machines at the first stage and multiple assembly machines at the second stage. They proposed some heuristic algorithms so as to minimize the sum of tardiness penalty and delivery costs.

Within the author's knowledge, no studies have addressed the total weighted tardiness objective on an assembly scheduling problem. This paper is to find effective and efficient schedules to minimize the total weighted tardiness of jobs. Section 2 presents the problem description, and Sections 3 and 4 characterize some solution properties to improve search effectiveness and efficiency. Moreover, four constructive heuristics, HEDD, HWT, WSPTMAX, and WSPTAVER, are derived in Section 5, lower-bound procedures are derived in Section 6, and a branch-and-bound algorithm is derived in Section 7. The performance of the proposed algorithms is evaluated through numerical experiments in Section 8.

#### 2. Problem Description

This study has n jobs to process and a job consists of m components which are processed in the first stage, and the m manufactured components are assembled in the second stage. Component k (k = 1, 2, …, m) of each job is processed in dedicated machine k, and the processing time ajk is required for component k of job j (j = 1, 2, …, n, k = 1, 2, ..., m). In addition, the assembly of job j can start in the assembly machine after the completion of all the components, and the required assembly time is bj (j = 1, 2, ..., n). The assembly system in this paper is presented in Figure 1.

It is assumed that more than one job in a machine cannot be processed at the same time, and no preemption is allowed in any processing. The completion time of a job is the completion time of the assembly at the second stage, which is expressed as Cj. dj denotes the due date for job j, and then the tardiness value of job j is computed as . Thus, the total weighted tardiness of n jobs is , where denotes the weight for job j.

This paper introduces Property 1 which can reduce the solution space to search, where the permutation schedule means that the production order is the same in all the machines.

Property 1. An optimal schedule exists among permutation schedules.

Proof. It can be easily shown by a pairwise job interchange augment.

#### 3. Identical Weight Case

Now, we characterize some dominance properties for finding better schedules for the problem. These properties can be used in some heuristic algorithms and a branch-and-bound algorithm. This section considers the situation where all the job weights are identical; that is, .

For the property derivation, assume that jobs p and q are processed consecutively. Denote by S a schedule in which job p is processed immediately before job q, while R denotes a schedule which is identical to S but job q is processed immediately before job p. Moreover, π denotes a set of jobs processed before jobs p and q in S and R. Some solution properties are presented for jobs p and q as follows.

Property 2. If the following conditions hold, job q should be processed immediately before p:

Proof. We prove that schedule S is dominated by R under conditions (1) and (2). First, we can easily see that the completion time of job q in schedule R is smaller than that of job q in S, so that the relation holds. Condition (1) implies , so that the following relation holds:If condition (2) is satisfied, the earliest start time of a job after jobs p and q on the assembly machine in R is no later than the earliest start time of a job after p and q in S. Therefore, we can conclude that schedule S is dominated by schedule R.

Property 3. If the following conditions hold, job q should be processed immediately before p:

Proof. As we discussed earlier in Property 2, condition (4) implies that the following relation holds:In condition (5), the right-hand-side term means the possible earliest start time after jobs p and q on the assembly machine. Condition (5) represents that the q ⟶ p sequence in R does not delay the processing of any jobs after jobs p and q. Therefore, we can conclude that schedule S is dominated by schedule R.

Property 4. If the following conditions hold, job q should be processed immediately before p:

Proof. Remind that the q ⟶ p sequence is in schedule R, while the p ⟶ q sequence is in schedule S. From equations (7) and (8), we can see that the completion times of jobs q and p in R are earlier than the completion times of p and q in S, respectively. Now, the following cases are considered:(i)Case 1. and .In this case, both jobs p and q are completed before each due date, so the following relation holds:From the above equation, we can see that holds.(ii)Case 2. and .In this case, both jobs p and q are tardy, so the following relation holds:(due to conditions (7)–(9))(iii)Case 3. and .From the condition, (due to conditions (7) and (9)) .(iv)Case 4. and .From the condition, (due to condition (8)) .From cases 1, 2, 3, and 4, we can say that holds.

Property 5. If the following conditions hold, job q should be processed immediately before p:where tk represents the completion time of the last-positioned job before jobs p and q on machine k.
The indices x, u, and represent the machine indices satisfying the following relation:

Proof. At first, we try to show . and can be expressed as follows:where denotes the completion time of the job immediately before jobs p and q on the assembly machine. From condition (15), the following relation is true:Then, in order to prove the relation , we have to show that the holds as follows:From condition (15), the relation is true. Thus, the above equation is true; then,From equations (18) and (20), the following equation holds:From condition (14) and equation (21), we have and . In order to say that schedule R is better than schedule S, due date should be examined additionally as in the two following cases:(i)Case 1. . Both p and q are tardy from condition (13), so the following relation holds:From equations (14) and (21), we can see that holds.(ii)Case 2. . From (13) and (14), the following relation holds:From (14), (21), and (32), we can say that holds according to Property 4.
Property 6 is an extension of Schaller , dominance property which is applied in this paper.

Property 6. If the following conditions hold, job q should be processed immediately before p:where denotes the number of jobs in the set .

Proof. If condition (24) is satisfied, the sum of tardiness of jobs p and q in R is smaller than that of S. In condition (25), the right-hand-side term represents the upper limit of the potential tardiness increment of the jobs after p and q. From condition (25), we can say that schedule S is dominated by schedule R.

#### 4. Nonidentical Weight Case

Now, we present some dominance properties for the situation where all the job weights would be nonidentical. Some solution properties are characterized for jobs p and q as follows.

Property 7. If the following conditions hold, job q should be processed immediately before p:

Proof. The proof is the same as in Property 2.

Property 8. If the following conditions hold, job q should be processed immediately before p:

Proof. The proof is the same as in Property 3.

Property 9. If the following conditions hold, job q should be processed immediately before p:

Proof. The proof can be made in a similar manner in Property 6.

#### 5. Constructive Heuristic Algorithms

This paper proposes four constructive heuristic algorithms. Let π denote a partial sequence whose processing order is already assigned. The heuristic algorithms in this paper are based on a scheme which selects one job that has not been yet allocated and attaches it to the end position of the partial sequence π. Let denote the completion time of the last-positioned job in π.

This paper proposes an EDD (Earliest Due Date) type heuristic. A specific procedure is as follows, where denotes the possible completion time of job j when attaching job j to the end position of π, and thus, , where denotes the completion time of the last-positioned job in machine k in π.

##### 5.1. Heuristic HEDD

Step 0. , , , for k = 1, 2, …, 0Step 1. Select a job with the smallest due date value for j ∉ πStep 2. Add the job to the end position of the set πStep 3. If , then terminate. Otherwise, go to Step 1

It is true that EDD-type heuristics generally perform well in the scheduling problems associated with the due date objective measures. However, in the total weighted tardiness as in this study, it may be difficult to guarantee the performance of EDD-type heuristics because the schedule must be changed according to the job weight. Therefore, this paper presents another heuristic based on the job weight as follows.

##### 5.2. Heuristic HWT

Step 0. , , , for k = 1, 2, …, 0Step 1. Select a job with the largest weight value for j ∉ πSteps 2 and 3. The same as in the heuristic HEDD

Heuristic HEDD can be expected to show excellent performance in situations sensitive to the due date parameters, while heuristic HWT can be expected to show better performance in the weight-sensitive situations. However, it can be said that the two heuristics now have limitations in that they do not consider processing times on machines. Now, this study presents two additional heuristics, called WSPTMAX and WSPTAVER, which consider processing times. The first heuristic WSPTMAX calculates the maximum value of the processing times of each job in the first stage and selects a job with the smallest value of the value . The second heuristic WSPTAVER calculates the average value of the processing times of each job in the first stage and selects a job with the smallest value of the value . The two heuristics have an inherent limitation that they do not consider due dates, but they can be said to be very advantageous in that they consider processing times and weights together. If the completion time can be accelerated considering the processing times, it can be said that due date parameters are satisfied as a result. The specific procedures of the two heuristics are presented as follows.

##### 5.3. Heuristic WSPTMAX

Step 0. , , , for k = 1, 2, …, 0 Calculate the value for all job jStep 1. Select a job with the smallest value for j ∉ πSteps 2 and 3. The same as in the heuristic HEDD

##### 5.4. Heuristic WSPTAVER

Step 0. , , , for k = 1, 2, …, 0 Calculate the value for all job jSteps 1, 2, and 3. The same as in the heuristic WSPTMAX

#### 6. A Lower-Bound Procedure

This section proposes a lower-bound procedure to evaluate the performance of the derived heuristic algorithms, and this lower bound can also be used in the branch-and-bound algorithm developed in a later part of this paper. In the following lower bound LB1, [y] denotes the y-th smallest number in the parameters ajk’s, bj’s, and ’s.

##### 6.1. Lower-Bound LB1

Step 0. Calculate the Dj value as follows:for j = 1, 2, …, n, where Step 1. Calculate LB1 as , where

In the above lower bound, denotes a lower-bound value of the completion time of any possible j-th positioned job:

(k) denotes the k-th positioned job in any schedule, and then,

To derive a lower-bound value of equation (31), we sequence values in nondecreasing order:

Therefore, we can say that the lower-bound LB1 provides a lower bound for the original problem. The complexity order of the lower bound LB1 is . Two more lower bounds are derived on Property 10, which was verified by Chu  and Kim .

Property 10. For the Cj and dj values satisfying , then the following relation holds:Based on Property 10, this paper proposes two more lower bounds LB2 and LB3. The second lower bound LB2 is calculated as follows; we can begin with the fact that the completion time of the first job, C1, cannot be larger than the following value under the condition :In inequality (34), the value is the possible earliest starting time of processing, and is the smallest assembly time among all the jobs. From (34), an additional relation for the x-th smallest completion time Cx holds as follows:From (34) and (35), we can derive the second lower bound LB2 as follows:Now, we show the process of deriving the third lower bound LB3; we can derive another relation that can replace equation (34) as follows:The value is the possible earliest starting time of processing in the assembly machine. From (37), an additional relation for the x-th smallest completion time Cx holds as follows:From (37) and (38), we derive the second lower bound LB3 as follows:

#### 7. A Branch-and-Bound Algorithm

This paper derives a branch-and-bound method to find the optimal solution for medium-size problems. The proposed branch-and-bound algorithm uses depth-first search as a branching rule, which means selecting the lowest node (i.e., the node with the largest partial permutation) in the branching tree when selecting a node. In the case of a tie, the node having the lowest lower bound is selected.

In the B&B algorithm, the initial upper bound is obtained from the minimum value among the four heuristics H1, HEDD, WSPTMAX, and WSPTAVER proposed in Section 5:

This paper uses the largest of the three values LB1, LB2, and LB3 as the lower bound for the branch-and-bound algorithm as follows:

##### 7.1. Branch-and-Bound Algorithm

Step 0 (initialization). , , . Calculate the value .Step 1 (lower bound). Calculate .Step 2 (branching).If , then go to Step 4.If , then go to Step 5.Select k such as , andStep 3 (dominance properties).If then go to Step 5.If any of Properties 29. holds, then go to Step 5.Go to Step 1.Step 4 (update the upper bound).If the last level of the enumeration tree is reached, the total weighted tardiness WTT is calculated.If , then .Step 5. The current node no longer needs to be searched:

If , then go to Step 2. Otherwise terminate the B&B algorithm.

The optimal solution is UB.

#### 8. Computational Experiments

This section measures the performance of the branch-and-bound algorithms and the heuristic algorithms through various numerical tests. For the test, the generated experimental data is set as follows:(1)The parameters ajk, bj, and (j = 1,…,n, k = 1,…,m) are generated from U (1, 100), where U (x, y) represents a uniform distribution with the numbers x and y.(2)The due dates dj (j = 1,…,n) are generated from , where P denotes a possible lower-bound value of the makespan value:

We conducted computational experiments with the number of jobs n in {5, 7, 9, 11, 13, 15} and the number of machines m in {2, 4, 6}. The coefficients RDD (relative range of due date) and TF (tardiness factor) are from the values {0.2, 0.5, 0.8} for the due date generation. 20 problems were randomly generated for 162 combinations of jobs and machines and RDD and TF values, and then a total of 3,240 problems were generated.

For the performance evaluation of the proposed branch-and-bound algorithm, the computational experiments are summarized in Table 1. In the table, “the number of nodes” refers to the total number of the generated nodes, and “computational time” refers to the computing time of the branch-and-bound algorithm. From the table, we can see that the optimal solutions were found within 1,200 seconds (CPU) for instances with up to about 15 jobs, and the average number of nodes and computational time of the B&B algorithm do not increase significantly as the number of machines increases. From the fact that the number of nodes and computational time decrease in some cases as the number of machines increases, we can possibly say that the lower-bound mechanism of the B&B algorithm is effective. Moreover, this paper tries to evaluate the performance of the dominance properties in Sections 3 and 4. We can see the performance summary of the B&B algorithm without any dominance properties applied. From the table, the computational time and the average number of search nodes increased significantly when the dominance properties were not applied, and specifically optimal solutions were not found for problems with 15 jobs within 20 minutes. From the comparative experiments, we can say that the dominance properties contribute to the performance of the B&B algorithm.

 n m B&B with dominance properties B&B without dominance properties Computational time The number of nodes Computational time The number of nodes Average Max Average Max Average Max Average Max 5 2 0.0007 0.004 34.1 109 0.0008 0.004 47.09 176 4 0.0008 0.005 33.8 96 0.0007 0.001 47.5 106 6 0.0008 0.003 32.7 105 0.0008 0.003 47.6 121 7 2 0.002 0.007 456.5 2198 0.003 0.008 861.7 3093 4 0.002 0.01 460.2 3012 0.003 0.01 920.1 3215 6 0.002 0.01 399.9 1737 0.003 0.01 841.0 2471 9 2 0.05 0.2 11893.4 75986 0.09 0.4 26816.9 162938 4 0.05 0.3 12159.7 92172 0.1 0.5 27496.8 170525 6 0.04 0.2 8685.7 39096 0.1 0.3 24986.6 87934 11 2 2.04 33.5 383616.1 8290345 4.9 50.0 1231820.6 15696940 4 2.3 44.2 401569.1 9476135 5.6 53.2 1281161.9 16173866 6 2.0 13.7 287125.6 2230298 5.3 27.3 1060709.6 6218744 13 2 170.6 1147.5 26198107.8 257803159 352.1 1180.4 73275993.7 247993027 4 132.4 1140.9 18033822.7 185741409 376.7 1187.9 69507001.4 249138079 6 140.0 1120.2 16511705.6 155593682 378.2 1185.2 62695670.7 204011606 15 2 910.2 1194.7 113697715.3 119565585 4 913.8 1144.9 96860705.2 133536399 6 907.4 1171.9 86202170.5 87887442
means the problem instances that the B&B algorithm cannot find the optimal solutions.

This paper conducted some experiments to investigate the performance change according to the RDD and TF values. The computational tests are carried out with n = 9 and m = 2 and are presented in Table 2. For each combination of RDD and TF values from {0.2, 0.5, 0.8}, 20 problems were randomly generated. The performance of the proposed B&B algorithm is summarized in Table 2. From the table, the performance of the B&B algorithm does not seem to fluctuate significantly according to the RDD and TF values.

 RDD TF Computational time The number of nodes Average Max Average Max 0.2 0.2 0.029 0.175 11626.45 36520 0.5 0.028 0.071 4143.25 12049 0.8 0.039 0.15 5650.75 23505 0.5 0.2 0.049 0.179 9536.7 34911 0.5 0.045 0.122 7760.05 23526 0.8 0.051 0.146 8294.2 23623 0.8 0.2 0.031 0.127 5547.8 21737 0.5 0.061 0.164 11497.05 28191 0.8 0.076 0.205 14115.05 39096

The performance of four heuristics HEDD, HWT, WSPTMAX, and WSPTAVER are given in Tables 3 and 4. In the table, the term “error” means the average error of the associated heuristic compared to the lower-bound value and is defined as follows:where WT (SH) is the total weighted tardiness of an obtained solution from a heuristic and LB is a lower-bound value in Section 6; thus, , and WT (SWorst) is the largest value among solutions obtained from four heuristics HEDD, HWT, WSPTMAX, and WSPTAVER. The smaller the value error, the better the performance of the associated heuristic. 180 instances are solved for each combination of jobs and machines. In the table, NO denotes the number of problems that the associated heuristic found the best solution among four heuristic values. NO value can reach up to 180, but since multiple heuristics can produce the same value, the sum of the NO values of the four heuristics can exceed 180.

 n m HEDD HWT WSPTMAX WSPTAVER Aver. error NO Aver. error NO Aver. error NO Aver. error NO 5 2 0.881 32 0.650 80 0.639 89 0.648 94 4 0.875 35 0.656 83 0.655 94 0.695 71 6 0.882 30 0.655 98 0.676 92 0.696 76 7 2 0.867 33 0.633 68 0.639 72 0.644 58 4 0.868 26 0.655 73 0.660 68 0.698 48 6 0.898 25 0.667 83 0.695 64 0.703 56 9 2 0.832 43 0.678 61 0.685 63 0.688 53 4 0.846 37 0.650 80 0.685 56 0.700 35 6 0.851 36 0.659 83 0.691 45 0.709 46 11 2 0.821 42 0.703 52 0.675 58 0.699 48 4 0.829 44 0.674 64 0.695 46 0.721 42 6 0.837 40 0.647 86 0.712 37 0.743 28 13 2 0.831 41 0.680 45 0.637 49 0.661 53 4 0.843 41 0.653 77 0.683 37 0.707 29 6 0.826 45 0.675 82 0.707 29 0.715 26 15 2 0.809 46 0.673 52 0.665 50 0.680 40 4 0.807 49 0.662 73 0.714 34 0.733 26 6 0.830 44 0.656 73 0.696 36 0.703 30
 RDD TF HEDD HWT WSPTMAX WSPTAVER Aver. error NO Aver. error NO Aver. error NO Aver. error NO 0.2 0.2 0.868 4 0.501 6 0.568 10 0.561 10 0.5 0.982 0 0.593 10 0.601 8 0.579 5 0.8 0.997 0 0.677 10 0.674 8 0.705 5 0.5 0.2 0.624 11 0.676 6 0.814 2 0.8 3 0.5 0.931 5 0.627 6 0.626 8 0.652 6 0.8 0.977 1 0.769 6 0.72 9 0.725 6 0.8 0.2 0.221 16 0.856 3 0.796 2 0.802 0 0.5 0.905 5 0.74 7 0.723 6 0.728 7 0.8 0.985 1 0.667 7 0.646 10 0.639 11

From the test result, it seems that HEDD does not show good performance and HWT shows relatively excellent results in the experiments. It may be because algorithms are more sensitive to the job weight fluctuations. Therefore, it can be interpreted to be more advantageous to process high-weight jobs preferentially.

Table 4 shows the performance change of the four heuristics according to the RDD and TF values. The experiment was carried out with the instances of 9 jobs and 2 machines, and 20 problems are solved for combinations of RDD and TF. In the table, NO value can reach up to 20, but the sum of NO values can exceed 20. From Table 4, it can be seen that the performance of the heuristics varies depending on the RDD and TF values. When the value of RDD increases and the value of TF decreases, the performance of HEDD tends to improve. The larger the RDD value, the greater the deviation of the due dates, and the smaller the TF value, the greater the chance of tardy jobs. It can be seen that HEDD heuristic performs better where the problem is more sensitive to due dates.

We conducted computational experiments with large-size jobs from 30 to 270 and the number of machines in {2, 4, 6}. 180 problems were randomly generated for 9 combinations of RDD and TF values, and then a total of 4,860 problems were generated. Table 5 shows the performance summary of the four heuristics. NO value can reach up to 180, but the sum of NO values can exceed 180. From the table, HWT shows the best performance among the four heuristics, while when the number of machines is 2, HEDD performs well.

 n m HEDD HWT WSPTMAX WSPTAVER Aver. error NO Aver. error NO Aver. error NO Aver. error NO 30 2 0.776 50 0.681 51 0.690 48 0.698 31 4 0.760 51 0.658 74 0.712 35 0.728 20 6 0.761 53 0.679 74 0.728 27 0.756 26 60 2 0.743 59 0.694 46 0.693 53 0.710 22 4 0.743 60 0.686 78 0.743 27 0.774 15 6 0.747 56 0.658 94 0.737 21 0.759 9 90 2 0.740 58 0.696 53 0.703 44 0.718 25 4 0.731 60 0.671 89 0.739 21 0.760 10 6 0.734 60 0.675 104 0.753 11 0.782 5 120 2 0.733 60 0.713 51 0.731 49 0.752 20 4 0.738 59 0.676 98 0.745 15 0.773 8 6 0.740 58 0.672 111 0.761 8 0.789 3 150 2 0.729 61 0.707 43 0.728 50 0.746 26 4 0.732 60 0.676 99 0.748 16 0.774 5 6 0.729 60 0.676 104 0.771 12 0.802 4 180 2 0.746 57 0.705 52 0.716 45 0.735 26 4 0.731 60 0.676 106 0.761 9 0.788 5 6 0.724 60 0.659 115 0.765 4 0.796 1 210 2 0.744 59 0.711 49 0.719 55 0.741 17 4 0.736 60 0.674 105 0.764 10 0.788 5 6 0.737 59 0.669 115 0.774 6 0.803 0 240 2 0.743 59 0.707 56 0.723 43 0.743 22 4 0.728 60 0.676 110 0.763 8 0.793 2 6 0.721 60 0.668 116 0.776 3 0.809 1 270 2 0.730 60 0.697 52 0.711 50 0.730 18 4 0.726 61 0.673 107 0.756 10 0.782 2 6 0.727 59 0.658 120 0.765 1 0.802 0

#### 9. Concluding Remarks

This paper considers a two-stage assembly scheduling problem where there are m dedicated machines in the first stage and an assembly machine in the second stage. This paper proposes some solution approaches to establish effective and efficient schedules to minimize the total weighted tardiness measure. In this paper, we develop some dominance properties and lower-bound procedures, derive a branch-and-bound algorithm for finding optimal solutions, and derive four heuristics. Computational tests are performed to evaluate the solution approaches in this paper. For problems with no more than 15 jobs, the B&B algorithm finds optimal solutions within 20 minutes, and the performance of the four heuristics seems to be similar. However, for problems with more than 30 jobs, HWT shows the best performance among the four heuristics, while when the number of machines is two, the performance of HEDD is better.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This study was supported by research funds from Dong-A University.

1. J. M. Framinan, P. Perez-Gonzalez, and V. Fernandez-Viagas, “Deterministic assembly scheduling problems: a review and classification of concurrent-type scheduling models and solution procedures,” European Journal of Operational Research, vol. 273, no. 2, pp. 401–417, 2019. View at: Publisher Site | Google Scholar
2. Q.-K. Pan, L. Gao, L. Xin-Yu, and F. M. Jose, “Effective constructive heuristics and meta-heuristics for the distributed assembly permutation flowshop scheduling problem,” Applied Soft Computing, vol. 81, p. 105492, 2019. View at: Publisher Site | Google Scholar
3. C.-Y. Lee, T. C. E. Cheng, and B. M. T. Lin, “Minimizing the makespan in the 3-machine assembly-type flowshop scheduling problem,” Management Science, vol. 39, no. 5, pp. 616–625, 1993. View at: Publisher Site | Google Scholar
4. C. N. Potts, S. V. Sevast'janov, V. A. Strusevich, L. N. van Wassenhove, and C. M. Zwaneveld, “The two-stage assembly scheduling problem: complexity and approximation,” Operations Research, vol. 43, no. 2, pp. 346–355, 1995. View at: Publisher Site | Google Scholar
5. A. M. A. Hariri and C. N. Potts, “A branch and bound algorithm for the two-stage assembly scheduling problem,” European Journal of Operational Research, vol. 103, no. 3, pp. 547–556, 1997. View at: Publisher Site | Google Scholar
6. X. Sun, K. Morizawa, and H. Nagasawa, “Powerful heuristics to minimize makespan in fixed, 3-machine, assembly-type flowshop scheduling,” European Journal of Operational Research, vol. 146, no. 3, pp. 498–896, 2003. View at: Publisher Site | Google Scholar
7. A. Tozkapan, Ö. Kırca, and C.-S. Chung, “A branch and bound algorithm to minimize the total weighted flowtime for the two-stage assembly scheduling problem,” Computers & Operations Research, vol. 30, no. 2, pp. 309–320, 2003. View at: Publisher Site | Google Scholar
8. I. S. Lee, “Minimizing total completion time in the assembly scheduling problem,” Computers & Industrial Engineering, vol. 122, pp. 211–218, 2018. View at: Publisher Site | Google Scholar
9. J. M. Framinan and P. Perez-Gonzalez, “The 2-stage assembly flowshop scheduling problem with total completion time: efficient constructive heuristic and metaheuristic,” Computers & Operations Research, vol. 88, pp. 237–246, 2017. View at: Publisher Site | Google Scholar
10. C. Talens, V. Fernandez-Viagas, P. Perez-Gonzalez, and J. M. Framinan, “New efficient constructive heuristics for the two-stage multi-machine assembly scheduling problem,” Computers & Industrial Engineering, vol. 140, p. 106223, 2020. View at: Publisher Site | Google Scholar
11. A. Allahverdi and f. S. Al-Anzi, “A PSO and a Tabu search heuristics for assembly scheduling problem of the two-stage distributed database application,” Computers and Operations Research, vol. 33, no. 4, pp. 1056–1080, 2006. View at: Publisher Site | Google Scholar
12. A. Allahverdi and H. Aydilek, “The two stage assembly flowshop scheduling problem to minimize total tardiness,” Journal of Intelligent Manufacturing, vol. 26, no. 2, pp. 225–237, 2015. View at: Publisher Site | Google Scholar
13. G. R. Ha, I. S. Lee, and S. H. Yoon, “A scheduling problem to minimize total tardiness in the two-stage assembly-type flowshop,” Journal of the Korean Operations Research and Management Science Society, vol. 33, no. 3, pp. 1–16, 2008. View at: Google Scholar
14. J. Y. Lee and J. Y. Bang, “A two-stage assembly-type flowshop scheduling problem for minimizing total tardiness,” Mathematical Problems in Engineering, vol. 2016, Article ID 6409321, 10 pages, 2016. View at: Google Scholar
15. H. Kazemi, M. M. Mazdeh, and M. Rostami, “The two stage assembly flow-shop scheduling problem with batching and delivery,” Engineering Applications of Artificial Intelligence, vol. 63, pp. 98–107, 2017. View at: Publisher Site | Google Scholar
16. J. Schaller, “Note on minimizing total tardiness in a two-machine flowshop,” Computers & Operations Research, vol. 32, no. 12, pp. 3273–3281, 2005. View at: Publisher Site | Google Scholar
17. C. Chu, “A branch-and-bound algorithm to minimize total tardiness with different release dates,” Naval Research Logistics, vol. 39, no. 2, pp. 265–283, 1992. View at: Publisher Site | Google Scholar
18. Y.-D. Kim, “A new branch and bound algorithm for minimizing mean tardiness in two-machine flowshops,” Computers & Operations Research, vol. 20, no. 4, pp. 391–401, 1993. View at: Publisher Site | Google Scholar