Research Article  Open Access
A Mathematical Model and a Simulated Annealing Algorithm for Balancing Multimanned Assembly Line Problem with SequenceDependent Setup Time
Abstract
Multimanned assembly lines have been widely applied to the industrial production, especially for largesized products such as cars, buses, and trucks, in which more than one operator in the same station simultaneously performs different tasks in parallel. This study deals with a multimanned assembly line balancing problem by simultaneously considering the forward and backward sequencedependent setup time (MALBPS). A mixedinteger programming is established to characterize the problem. Besides, a simulated annealing algorithm is also proposed to solve it. To validate the performance of the proposed approaches, a set of benchmark instances are tested and the lower bound of the proposed problem is also given. The results demonstrated that the proposed algorithm is quite effective to solve the problem.
1. Introduction
As floworiented production systems, assembly lines have been widely used in the industrial production of high quantity standardized commodities since Ford developed such a line in 1913. With the usage and spread of assembly lines, a combinatory problem, named assembly line balancing problem (ALBP), has aroused great interest of researchers [1–3]. As the simplest problem of ALBP, the simple assembly line balancing problem (SALBP) is to assign tasks to an ordered sequence of stations such that the precedence constraint and the cycle time constraint are satisfied, and one or more objectives are optimized. In terms of the mathematical complexity, the SALBP is strongly NPhard, since it can be subsumed as a special case of the bin packing problem, which was being proved as a NPhard problem [4]. Thus, numerous researchers have been devoted to developing various approaches including exact algorithm, heuristics, and metaheuristics to solve the problem [5, 6].
However, when assembly lines are applied to produce largesized products such as cars, trains, or aircraft, there is enough space to assign two or more workers to each station; the SALBP is no longer suitable since it assumes that only one operator is allowed in one station. In such conditions, the multimanned assembly line balancing problem (MALBP) is proposed to bridge the gap by resuming that more than one worker is allowed in one station, as shown in Figure 1. Another very similar problems is twosided assembly line balancing problem (TALBP) firstly addressed by [7]. However, the main difference between TALBP and MALBP is that only two workers can be in the same multimanned workstation and there are preferred operation directions of tasks that restrict the assignments of them to workstations [8]. Thus, the TALBP can be regarded as a special case of MALBP.
(a)
(b)
As an extension of SALBP, the MALBP is also NPhard problem and an increasing number of approaches have been designed to solve the problem. To the best of authors’ knowledge, Dimitriadis [9] is the first researcher to address the MALBP by proposing a twolevel heuristic algorithm to study a real case from an automobile assembly plant. In terms of the characteristics of the objective functions, the MALBP can be divided into three versions: MABLPI to minimize the number of workers/stations, MALBPII to minimize the cycle time, and MALBPC to minimize the total production cost.
Regarding the MABLPI, Fattahi et al. [10] firstly developed a MIP model and an antcolonybased algorithm, and then Yilmaz and Yilmaz [11] corrected their mathematical model. Kellegöz and Toklu [12] modified the branchandbound algorithm to minimize the length of the line. Then, the same author developed a priorityrule based constructive heuristic approach for the same problem [13]. Roshani et al. [14] designed a simulated annealing algorithm to solve the problem by simultaneously optimizing the line efficiency, the line length, and the smoothness index. Recently, Kellegöz [15] modified the simulated annealing algorithm by performing on Gantt representations of solutions, and a new MIP model requiring less number of variables and constraints was also developed. Michels et al. [16] proposed a benders’ decomposition algorithm to minimize the number of workers as the primary objective. For other extension of MALBPI, Chen et al. [17] additionally considered the resource constraint. Şahin and Kellegöz [18] assumed that workers could walk between different stations. Lopes et al. [19] defined the flexible multimanned assembly lines with flexible station frontiers, in which multiple workers in the same station were allowed to start operating as possible as they can. Although the MALBPII is a major practice for reconfiguration of the installed assembly lines, research on the MALBPII is much less than that on the MALBPI. Roshani and Giglio [20] presented a mathematical model and two simulated annealingbased algorithms for solving the MALBPII: one solves it directly and another one solves it by repeatedly solving the MALBPI. Besides the above mentioned two timeoriented objectives, the costoriented objective gradually gets more attention since the production cost has been the key to win in competitions among factories. Regarding the MALBPC, Roshani and Giglio [21] presented a MIP mathematical model to describe the problem. Giglio et al. [22] assumed workers were skilled and a MIP model was also built. Şahin and Kellegöz [23] additionally considered the resource constraint in the MALBPC; a MIP model and a particle swarm optimization algorithm were both developed to solve it.
Furthermore, despite it being certainly an argument whose setups are ubiquitous to the realistic work environment and cannot be ignored for decisionmaking [24], most researchers still assumed setups as negligible in the literature. This phenomenon may be caused by the following: (i) the influence of setups has been reduced to a large degree by using some advanced manufacturing technologies, such as flexible manufacturing system; (ii) in some real cases, the setup time is ignored because it is too small; (iii) researchers simplified their problem by assuming the setup time as an assignment constraint (e.g., incompatible task). However, there are still many environments where setup time is significant, especially in the case when a station is operated at or near full capacity. The time required to perform a setup activity is called as setup time. Setup time can be classified as sequencedependent setup time and sequenceindependent setup time. In assembly lines, sequencedependent setup time occurs when the setup time of a task depends on which task was set up on the station prior to operate that task. With respect to the sequenceindependent setup time, the sequencedependent setup time needs more consideration (such as scheduling tasks intra each station). Hence, treating this kind of setup time separately from processing time allows operations to be performed simultaneously and hence improves performance. Furthermore, in multimanned assembly lines, setups usually cannot be ignored and should be considered more carefully with respect to the traditional singlemanned assembly lines. The reason of this is that the largesized products are often assembled in the former assembly lines. Thus, the setups tend to be large enough with respect to their task time, such as the travel time of operators caused by its moving around the work piece (largesized) to perform the assigned tasks [25].
Andres et al. [26] made the first attempt to solve the assembly line balancing problem with setups (ALBPS); they also proposed a binary linear program model and a greedy randomized adaptive search procedure. Since then, various metaheuristics were proposed for solving the problem [27–29]. Then, Scholl et al. [30] extended the problem by distinguishing the forward and backward setup time; a MIP model and several heuristics were also developed. For other extensions of ALBPS, Nazarian et al. [31] presented a MIP model for the multimodel ALBPS. Sahin and Kellegoz [32] defined the crossover setup time in ushaped assembly lines. Özcan [33] defined the lineswitching setup time in parallel assembly lines; a binary linear program model and a simulated annealing algorithm were proposed to solve it. Akpinar et al. [34] developed a hybrid antcolony optimization algorithm for solving the mixedmodel ALBPS; then, Akpinar and Baykasoğlu [35, 36] extended this problem by distinguishing the forward and backward setup time (mALBPS). Ozcan and Toklu [37] firstly considered the twosided assembly line balancing problem with setups (TALBPS) and developed a MIP model and a heuristic approach to solve it. Janardhanan et al. [38] also extended the TALBPS to the robotic twosided assembly lines by proposing a MIP model and a metaheuristic migrating bird optimization algorithm. Aghajani et al. [39] extended the TALBPS to mixedmodel twosided robotic assembly lines, and they proposed a MIP model and a simulated annealing algorithm. Furthermore, Esmaeilbeigi et al. [40] developed three formulations for the ALBPS and designed several possible improvements in the form of valid inequalities and preprocessing approaches. Akpinar et al. [41] improved the model for the ALBPS and mALBPS and an exact procedure and introduced the benders decomposition algorithm.
Although a great deal of research has been devoted to various approaches for solving MALBP, according to our best knowledge, no published paper on MALBP in the literature has simultaneously considered forward and backward setup time before.
In this study, the MALBP is extended by considering the sequencedependent setup time to minimize the number of workers as primary objective and minimize the number of stations as secondary objective (MALBPSI). A mixedinteger programming mathematical model is built to characterize the MALBPS, and a metaheuristic algorithm based on simulated annealing (SA) approach is also developed to solve the problem.
The rest of this article is organized as follows. In Section 2, the problem to solve is formalized and the MIP model is also presented. Besides, an example of the proposed problem is given. Section 3 is devoted to the description of the proposed SA algorithm. In Section 4, the design of experiment is presented and the results are discussed. The conclusion and future direction are given in Section 5.
2. The MALBPSI
In this section, the MALBPSI is described in detail and the problem assumptions are listed. Finally, before calculating the lower bound, the proposed MIP formulation is developed.
2.1. Problem Definition
A series of multimanned stations () are utilized on the paced straight assembly lines to produce single model products. A set of workers are assigned to each multimanned station. A set of tasks are being assigned to workers and stations to minimize the number of workers and the number of multimanned stations without violating the cycle time constraint and the precedence constraint. As depicted in Figure 2, worker 1, worker 2, and worker 3 are assigned to multimanned station 1; worker 4 and worker 5 are assigned to multimanned station 2. Besides, in such a multimanned assembly line, the sequencedependent idle time may occur. For example, task is delayed by its predecessor task , which is operated by different workers in the same multimanned station.
(a)
(b)
In assembly lines, setup time may occur in two ways: the forward and backward setup time. As we can see from Figure 2, when a task is immediately performed before another task operated by the same worker at the same multimanned station in the same cycle, then a forward setup occurs for the same work piece to perform task and a forward setup time is added to the finish time of task . Furthermore, when a task is the last task operated by a worker and, in the next cycle, task is the first task operated by the same worker at the same multimanned station, then a backward setup occurs for the next work piece to perform task , and a backward setup time is added to compute the global station time.
Moreover, the sequencedependent idle time can be used for dealing with setup operations. As depicted in Figure 2, the sequence idle time occurs in worker 1. Thus, the finish time of task is calculated as . Besides, the station idle time also can be used to deal with backward setup operations. As depicted in Figure 2, the station idle time occurs in worker 1 and a backward setup time also occurs; if , then constraint (8) is satisfied. Otherwise, the cycle time constraint is violated.
2.2. Problem Assumptions
The problem assumptions of the MALBPSI are listed as below. Task time, setup time, and precedence diagram are deterministic in nature and known in advance All stations are equally equipped and all workers are assumed having the same ability to perform any tasks More than one worker is allowed to be assigned to each station Forward and backward setup time may occur between two adjacency tasks The buffers or WIP are not allowed
2.3. Notations
The notation is given in Table 1.

2.4. The Mathematical Model
The mathematical model proposed by [10, 11, 37] are extended to develop a MIP model for the MALBPSI in this study. Task is assumed as a virtue node with zero task time, and it is a final node of the precedence graph; thus, the station, which task is assigned to, is the final station. The model is given as follows:
Objective function (1) minimizes the number of workers as the primary objective and minimizes the number of stations as the secondary objective. Constraint (2) ensures that each task is assigned to one position of one station . Constraint (3) ensures that at most one task will be assigned to one position of one station . Constraint (4) ensures that the position will be opened in increasing order. Constraint (5) ensures that all precedence relations among tasks are satisfied. Constraints (6)–(8) control the sequencedependent finish time of tasks. If task and its immediate predecessor task are assigned to the same station j, then constraint (6) becomes . Constraints (7) and (8) ensure when the forward or backward setup occurs; then, the forward or backward setup time must be considered. When two tasks are assigned to the successive position in the same cycle of a station, then constraint (7) becomes . When two tasks are assigned to the successive position in the next cycle of a station, then constraint (8) becomes . Constraint (9) ensures that if two tasks are assigned to two adjacent positions of station , then will be equal to one. Constraint (10) ensures that if task is the last task of station , then will be equal to one. Constraint (11) provides us to determine the backward setup between the last task and the first task of a station. Constraint (12) ensures that if any task has been assigned to station , then will be equal to one. Constraint (13) ensures that if workers are assigned to station j, then will be equal to one. Constraint (14) observes the sequence of workers’ index in a multimanned station. Constraint (15) ensures that if no task is assigned to station , then will be equal to zero. Constraint (16) observes the sequence of stations’ index in lines. Constraints (17) and (18) ensure that the range of the finish time of task is between its completion time and the cycle time. Constraints (19)–(23) are the internality constraints. The lower bound of the problem is given in Appendix A.
2.5. An Example to Illustrate the MIP Model
The Merten problem [42] with or without setup time are both solved optimally by using the MIP model. The setup time is generated in Section 4 and the detailed data is given in Table 2. The cycle time is set as seven and the is set as three. As shown in Figure 3, when considering setups, one more worker is needed with respect to the problem without setups.

(a)
(b)
3. Proposed SA Algorithm for MALBPSI
As an extension of MALBP, the MALBPS is also strongly NPhard problem. Thus, it is necessary to develop a heuristic or metaheuristicbased algorithm to solve largesized problem. In this paper, a simulated annealing (SA) approach is proposed for solving the MALBPSI. Since the SA algorithm was introduced by Kirkpatrick et al. [43] as an iterative random search technique, it has widely been used to solve various combinatorial optimization problems including general assembly line balancing problem [14, 15, 20, 44, 45]. Basically, the SA algorithm is a local searchbased metaheuristic, which derives its acceptance mechanism from the annealing process to let the current solution escape from local optima. The detailed procedure of the proposed SA algorithm is given below.
3.1. Initial Solution
Considering that the number of stations is uncertain and it is essential to determine the task sequence in the proposed problem, a prioritybased coding method is adopted, in which solutions are constructed according to a priority list (PL) of tasks. Then, the initial solution is randomly generated as a sequence between 1 to by a uniform distribution (1, 2, 3, 4, 5, 6, 7), as shown in Figure 4. To obtain a feasible solution, the assignable task with the lowest priority value is being selected and then it is being assigned to a worker according to some given rules as Section 3.2. Then, the process continues until all tasks are assigned.
3.2. Building a Feasible Solution
A feasible balancing solution is to determine how to assign works to stations and how to assign tasks to workers without violating the precedence constraint and the cycle time constraint. The procedure to build a feasible solution is given as Algorithm 1 in Figure 5. An example is also illustrated in Table 3. The procedure to calculate the finish time () of a task is given as Algorithm 2 in Figure 6. The rules of accepting the task assignment to current multimanned station are defined as follows (if one of the following conditions (a or b) is fulfilled):

3.3. Neighborhood Generation
Neighborhood structures are main approaches to produce new solution in the SA algorithm by moving from a solution to its neighborhood one. For the proposed SA algorithm, the neighborhood structures are designed the same as [14] swap and insert based operators. An example is explained in Figure 7.
3.4. Objective Function
In this paper, three factors are taken into consideration including the line efficiency (), the line length (), and smoothing index () for the MALBPSI. Moreover, these criterions are combined together to build a single objective function by using a minimum deviation method [14]. Let , , and , which are obtained from a initial solution, be the least desirable objective value of ,, and . The objective function of the proposed algorithm is formulated as follows:where , , and , are, respectively, the most desirable objective value of ,, and . For the problem of MALBPSI, the value of the is set as 100, is set as , and is set as 0.
3.5. Description of Proposed SA Algorithm
The proposed SA algorithm starts with an initial solution (Section 3.1) in a predetermined temperature (initialized as ) and navigates around it through some predetermined neighborhood structures (Section 3.3) to find better neighborhood solution. With the searching process of the algorithm, the temperature decreases during the iteration as , where is the predetermined cooling rate. In each temperature, the remains the same for a period of time (, the Marko chain length), and then it falls through the algorithm, the probability decreases, and the result converges to the best solution. The neighborhood solution will be accepted as the probability , where is the difference of fitness between the current solution and its neighbourhood solution. The algorithm stops when the is lower than a predetermined value . Then, output the best found solution. The procedure is given in Figure 8.
4. Computational Results
In this section, there are mainly two folds. Firstly, the test instances are solved by using the MIP model. The model is solved by IBM ILOG CPLEX 12.6.3. Secondly, the test instances in the literature are solved by the proposed SA algorithm, which is coded by MATLAB 2016a software. All experiments are executed on PC with Inter® Core™ i3, 3.4 GHz processor, and 4.0 GB of RAM.
The test instances are a group of wellknown Talbot data set including 64 problems from [42]. The forward setup time is randomly generated as for problems with lowlevel setups and for problems with highlevel setups. The backward setup time is randomly generated as for problems with problems with lowlevel setups and for problems with highlevel setups. The number of maximum allowed workers is set as two. All the parameters of the proposed SA algorithm are determined through preliminary experiments by using the Taguchi method [46]. To determine the values of parameters including , three test instances are randomly selected as Mansoor (), Kilbridge, and Arcus problem(). Three levels of , three levels of , and three levels of are all tested. The value of is fixed to . According to the number of levels and factors, the Taguchi method L9 is used for the adjustment of the parameters. Each test problem is solved ten times, and a performance measure is defined aswhere and are the best solution obtained for a given instance and the objective value obtained for a trial, respectively. According to the results of the statistical analysis, the value of are determined as , and , respectively.
4.1. The Result of the MIP Model
The results of the MIP model for solving the MALBPSI with different level setups are given in Table 4. The model is coded by OPL language and each test stops until the running time is larger than 7200 seconds. As we can observe from the table, only very limited smallsized (n ≤ 20) test instances can be solved optimally in the time limit. The best found number of workers/stations () and the CPU time are also reported. The bold text in the table shows that compared with the problem without setups, the same problem with setups need more stations or workers, and the higher the setups are, the more stations or workers are needed. means that not optimal results, but the best feasible results found by the model in time limit are reported.

4.2. The Result of the Proposed Algorithm
In this section, the 62 test instances are solved by using the proposed SA algorithm. Each instance is solved ten times and the results are reported in Table 5. The and mean the best and mean found number of stations/workers (), respectively. The is the optimal value of the ALBPI from [42], which can be regarded as the lower bound of the number of workers for the MALBPI. The is the optimal result of number of stations/workers () for MALBPI obtained by [10]. The is the best reported number of stations/workers () in the literature for the MALBPI without setups [14, 15]. The is the lower bound of the number of workers given in this paper for the MALBPSI.

In order to evaluate the effectiveness of the proposed algorithm, the index represents the gap between the best found number of worker and the lower bound of MALBPSI and can be calculated as . Another index represents the gap between the best found number of worker and the and can be calculated as . Moreover, represents the difference between the best and average found number of workers and represents the variance of number of workers found by the SA algorithm in ten times’ runs.
The results in Table 5 show that the proposed SA algorithm is able to find good solutions to MABLPSI in a reasonable time. Besides, a summary of computational results is given in Table 6. For all test instances, the average values of from for problems with lowlevel and highlevel setups are 1.23% and 0.98%, respectively. These results show that the proposed algorithm produces very close results to the lower bound. Besides, among 50 test problems in Table 5, the proposed SA algorithm obtained 46 and 47 optimal solutions equal to the lower bound for problems with lowlevel and highlevel setups, respectively. In addition to this, the average value of the and the average value of best found for problem with highlevel setups are both larger than the value for problems with lowlevel setups. The results show that the setups may increase the number of utilized workers, and the number of workers will be increased with the level of setups. Furthermore, the average value of and for all test problems are 0.007 and 0.006, respectively, which show that the proposed SA algorithm has a high rate of convergence and stability.

5. Conclusion
In this paper, the typeI problem of balancing multimanned assembly lines with sequencedependent setup time is introduced and characterized. A MIP and a simulated annealing algorithm are proposed to solve the problem. To demonstrate the efficiency of SA algorithm and MIP model, a computational experiment is conducted. The results show that the proposed approach is effective and successful for the problem. Furthermore, according to our best knowledge, this paper is the first study on the multimanned assembly line balancing problem with forward and backward setup time; thus, this paper gives a good starting point for future studies. The following directions may be studied in future: (1) more recent developed metaheuristics and exact algorithm should be developed; (2) other extension and assignment constraints, such as ushaped lines and zoning constraints should be considered according to the need of real line; (3) more effective procedure to build a feasible solution for the proposed should be developed; (4) tighter lower bound should be given; (5) dynamic job rotation tools can be combined with the proposed problem [47]; (6) multipleuser defined criteria should be considered and a real case should be studied [47].
Appendix
A. The Procedure to Calculate the Lower Bound of MALBPSI
The optimal solution for ALBPI from [42] is set as . Andres, Miralles, and Pastor [26] designed a lower bound for ALBPSI; then, this method is being modified to develop the lower bound of MALBPSI () and the detail of the procedure is given below. Step 1: calculate as If then go to Step 5; otherwise, = Step 2: calculate minimum forward setup times as . The method is to add lowest forward setup times between tasks using logical operator OR. Then, calculate the minimum backward setup times as . The method is to add lowest backward setup times between tasks. Then, add all possible setup times as . Then, calculate as Step 3: Step 4: if and go to Step 5; otherwise, go to Step 2. Step 5: output .
Data Availability
The benchmark datasets that support the findings of this study are previously reported by Scholl et al. [30] and are available at https://assemblylinebalancing.de/salbp/benchmarkdatasets1993/. Other 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 research was supported by the National Natural Science Foundation of China (Grants no.51675450) and Sichuan Science and Technology Program (Grant nos. 2019YFG0300 and 2019YFG0285).
References
 C. Becker and A. Scholl, “A survey on problems and methods in generalized assembly line balancing,” European Journal of Operational Research, vol. 168, no. 3, pp. 694–715, 2006. View at: Publisher Site  Google Scholar
 E. Erel and S. C. Sarin, “A survey of the assembly line balancing procedures,” Production Planning & Control, vol. 9, no. 5, pp. 414–434, 1998. View at: Publisher Site  Google Scholar
 A. Scholl and C. Becker, “Stateoftheart exact and heuristic solution procedures for simple assembly line balancing,” European Journal of Operational Research, vol. 168, no. 3, pp. 666–693, 2006. View at: Publisher Site  Google Scholar
 M. R. Garey and D. S. Johnson, “Computers and intractability; a guide to the theory of NPcompleteness,” 1990. View at: Google Scholar
 O. Battaïa and A. Dolgui, “A taxonomy of line balancing problems and their solutionapproaches,” International Journal of Production Economics, vol. 142, no. 2, pp. 259–277, 2013. View at: Publisher Site  Google Scholar
 P. Sivasankaran and P. Shahabudeen, “Literature review of assembly line balancing problems,” The International Journal of Advanced Manufacturing Technology, vol. 73, no. 9–12, pp. 1665–1694, 2014. View at: Publisher Site  Google Scholar
 J. J. Bartholdi, “Balancing twosided assembly lines: a case study,” International Journal of Production Research, vol. 31, no. 10, pp. 2447–2461, 1993. View at: Publisher Site  Google Scholar
 A. Roshani, P. Fattahi, A. Roshani, M. Salehi, and A. Roshani, “Costoriented twosided assembly line balancing problem: a simulated annealing approach,” International Journal of Computer Integrated Manufacturing, vol. 25, no. 8, pp. 689–715, 2012. View at: Publisher Site  Google Scholar
 S. G. Dimitriadis, “Assembly line balancing and group working: a heuristic procedure for workers' groups operating on the same product and workstation,” Computers & Operations Research, vol. 33, no. 9, pp. 2757–2774, 2006. View at: Publisher Site  Google Scholar
 P. Fattahi, A. Roshani, and A. Roshani, “A mathematical model and ant colony algorithm for multimanned assembly line balancing problem,” The International Journal of Advanced Manufacturing Technology, vol. 53, no. 1–4, pp. 363–378, 2011. View at: Publisher Site  Google Scholar
 H. Yilmaz and M. Yilmaz, “Erratum to: note to: a mathematical model and ant colony algorithm for multimanned assembly line balancing problem,” The International Journal of Advanced Manufacturing Technology, vol. 89, no. 5–8, p. 1941, 2017. View at: Publisher Site  Google Scholar
 T. Kellegöz and B. Toklu, “An efficient branch and bound algorithm for assembly line balancing problems with parallel multimanned workstations,” Computers & Operations Research, vol. 39, no. 12, pp. 3344–3360, 2012. View at: Publisher Site  Google Scholar
 T. Kellegöz and B. Toklu, “A priority rulebased constructive heuristic and an improvement method for balancing assembly lines with parallel multimanned workstations,” International Journal of Production Research, vol. 53, no. 3, pp. 736–756, 2015. View at: Publisher Site  Google Scholar
 A. Roshani, A. Roshani, A. Roshani, M. Salehi, and A. Esfandyari, “A simulated annealing algorithm for multimanned assembly line balancing problem,” Journal of Manufacturing Systems, vol. 32, no. 1, pp. 238–247, 2013. View at: Publisher Site  Google Scholar
 T. Kellegöz, “Assembly line balancing problems with multimanned stations: a new mathematical formulation and Gantt based heuristic method,” Annals of Operations Research, vol. 253, no. 1, pp. 377–404, 2017. View at: Publisher Site  Google Scholar
 C. Michels, A. Sato, and L. Magatão, “A Benders’ decomposition algorithm with combinatorial cuts for the multimanned assembly line balancing problem,” European Journal of Operational Research, vol. 22, 2019. View at: Publisher Site  Google Scholar
 Y.Y. Chen, C.Y. Cheng, and J.Y. Li, “Resourceconstrained assembly line balancing problems with multimanned workstations,” Journal of Manufacturing Systems, vol. 48, pp. 107–119, 2018. View at: Publisher Site  Google Scholar
 M. Şahin and T. Kellegöz, “Balancing multimanned assembly lines with walking workers: problem definition, mathematical formulation, and an electromagnetic field optimisation algorithm,” International Journal of Production Research, vol. 19, 2019. View at: Publisher Site  Google Scholar
 T. C. Lopes, G. Vidal Pastre, A. S. Michels, and L. Magatão, “Flexible multimanned assembly line balancing problem: model, heuristic procedure, and lower bounds for line length minimization,” Omega, vol. 34, 2019. View at: Publisher Site  Google Scholar
 A. Roshani and D. Giglio, “Simulated annealing algorithms for the multimanned assembly line balancing problem: minimising cycle time,” International Journal of Production Research, vol. 55, no. 10, pp. 2731–2751, 2017. View at: Publisher Site  Google Scholar
 A. Roshani and D. Giglio, “A mathematical programming formulation for costoriented multimanned assembly line balancing problem,” IFACPapersOnLine, vol. 48, no. 3, pp. 2293–2298, 2015. View at: Publisher Site  Google Scholar
 D. Giglio, M. Paolucci, A. Roshani, and F. Tonelli, “Multimanned assembly line balancing problem with skilled workers: a new mathematical formulation,” Spectrum, vol. 50, 2017. View at: Google Scholar
 M. Şahin and T. Kellegöz, “A new mixedinteger linear programming formulation and particle swarm optimization based hybrid heuristic for the problem of resource investment and balancing of the assembly line with multimanned workstations,” Computers & Industrial Engineering, vol. 133, pp. 107–120, 2019. View at: Google Scholar
 A. Allahverdi, “The third comprehensive survey on scheduling problems with setup times/costs,” European Journal of Operational Research, vol. 246, no. 2, pp. 345–378, 2015. View at: Publisher Site  Google Scholar
 G. Michalos, “Dynamic job rotation for workload balancing in human based assembly systems,” CIRP Journal of Manufacturing Science and Technology, vol. 2, no. 3, pp. 153–160, 2010. View at: Publisher Site  Google Scholar
 C. Andrés, C. Miralles, and R. Pastor, “Balancing and scheduling tasks in assembly lines with sequencedependent setup times,” European Journal of Operational Research, vol. 187, no. 3, pp. 1212–1223, 2008. View at: Publisher Site  Google Scholar
 N. Hamta, S. M. T. Ghomi, F. Jolai, and M. Akbarpour Shirazi, “A hybrid PSO algorithm for a multiobjective assembly line balancing problem with flexible operation times, sequencedependent setup times and learning effect,” International Journal of Production Economics, vol. 141, no. 1, pp. 99–111, 2013. View at: Publisher Site  Google Scholar
 S. A. SeyedAlagheband, S. M. T. F. Ghomi, and M. Zandieh, “A simulated annealing algorithm for balancing the assembly line type II problem with sequencedependent setup times between tasks,” International Journal of Production Research, vol. 49, no. 3, pp. 805–825, 2011. View at: Publisher Site  Google Scholar
 A. Yolmeh and F. Kianfar, “An efficient hybrid genetic algorithm to solve assembly line balancing problem with sequencedependent setup times,” Computers & Industrial Engineering, vol. 62, no. 4, pp. 936–945, 2012. View at: Publisher Site  Google Scholar
 A. Scholl, N. Boysen, and M. Fliedner, “The assembly line balancing and scheduling problem with sequencedependent setup times: problem extension, model formulation and efficient heuristics,” Spectrum, vol. 35, no. 1, pp. 291–320, 2013. View at: Publisher Site  Google Scholar
 E. Nazarian, J. Ko, and H. Wang, “Design of multiproduct manufacturing lines with the consideration of product change dependent intertask times, reduced changeover and machine flexibility,” Journal of Manufacturing Systems, vol. 29, no. 1, pp. 35–46, 2010. View at: Publisher Site  Google Scholar
 M. Sahin and T. Kellegoz, “Increasing production rate in Utype assembly lines with sequencedependent setup times,” Engineering Optimization, vol. 49, no. 8, pp. 1401–1419, 2017. View at: Publisher Site  Google Scholar
 U. Özcan, “Balancing and scheduling tasks in parallel assembly lines with sequencedependent setup times,” International Journal of Production Economics, vol. 213, pp. 81–96, 2019. View at: Publisher Site  Google Scholar
 S. Akpinar, G. M. Bayhan, and A. Baykasoglu, “Hybridizing ant colony optimization via genetic algorithm for mixedmodel assembly line balancing problem with sequence dependent setup times between tasks,” Applied Soft Computing, vol. 13, no. 1, pp. 574–589, 2013. View at: Publisher Site  Google Scholar
 Ş. Akpinar and A. Baykasoğlu, “Modeling and solving mixedmodel assembly line balancing problem with setups. Part I: a mixed integer linear programming model,” Journal of Manufacturing Systems, vol. 33, no. 1, pp. 177–187, 2014. View at: Publisher Site  Google Scholar
 Ş. Akpinar and A. Baykasoğlu, “Modeling and solving mixedmodel assembly line balancing problem with setups. Part II: a multiple colony hybrid bees algorithm,” Journal of Manufacturing Systems, vol. 33, no. 4, pp. 445–461, 2014. View at: Publisher Site  Google Scholar
 U. Özcan and B. Toklu, “Balancing twosided assembly lines with sequencedependent setup times,” International Journal of Production Research, vol. 48, no. 18, pp. 5363–5383, 2010. View at: Publisher Site  Google Scholar
 M. N. Janardhanan, Z. Li, G. Bocewicz, Z. Banaszak, and P. Nielsen, “Metaheuristic algorithms for balancing robotic assembly lines with sequencedependent robot setup times,” Applied Mathematical Modelling, vol. 65, pp. 256–270, 2019. View at: Publisher Site  Google Scholar
 M. Aghajani, R. Ghodsi, and B. Javadi, “Balancing of robotic mixedmodel twosided assembly line with robot setup times,” International Journal of Advanced Manufacturing Technology, vol. 74, no. 5–8, pp. 1005–1016, 2014. View at: Publisher Site  Google Scholar
 R. Esmaeilbeigi, B. Naderi, and P. Charkhgard, “New formulations for the setup assembly line balancing and scheduling problem,” Spectrum, vol. 38, no. 2, pp. 493–518, 2016. View at: Publisher Site  Google Scholar
 S. Akpinar, A. Elmi, and T. Bektaş, “Combinatorial Benders cuts for assembly line balancing problems with setups,” European Journal of Operational Research, vol. 259, no. 2, pp. 527–537, 2017. View at: Publisher Site  Google Scholar
 F. B. Talbot and J. H. Patterson, “An integer programming algorithm with network cuts for solving the assembly line balancing problem,” Management Science, vol. 30, no. 1, pp. 85–99, 1984. View at: Publisher Site  Google Scholar
 S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, “Optimization by simulated annealing,” Science, vol. 220, no. 4598, pp. 671–680, 1983. View at: Publisher Site  Google Scholar
 Y.Y. Chen, “A hybrid algorithm for allocating tasks, operators, and workstations in multimanned assembly lines,” Journal of Manufacturing Systems, vol. 42, pp. 196–209, 2017. View at: Publisher Site  Google Scholar
 A. Roshani, “Mixedmodel multimanned assembly line balancing problem: a mathematical model and a simulated annealing approach,” International Journal of Production Research, vol. 37, 2017. View at: Google Scholar
 K.L. Tsui, “AN overview OF taguchi method and newly developed statistical methods for robust design,” IIE Transactions, vol. 24, no. 5, pp. 44–57, 1992. View at: Publisher Site  Google Scholar
 G. Michalos, A. Fysikopoulos, S. Makris, D. Mourtzis, and G. Chryssolouris, “Multi criteria assembly line design and configuration  an automotive case study,” CIRP Journal of Manufacturing Science and Technology, vol. 9, pp. 69–87, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Wucheng Yang and Wenming Cheng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.