A multi-objective genetic algorithm for a mixed-model assembly U-line balancing type-I problem considering human-related issues, training, and learning

Rabbani, Masoud; Montazeri, Mona; Farrokhi-Asl, Hamed; Rafiei, Hamed Article A multi-objective genetic algorithm for a mixed-model assembly U-line balancing type-I problem considering human-related issues, training, and learning Journal of Industrial Engineering International Provided in Cooperation with: Islamic Azad University (IAU), Tehran Suggested Citation: Rabbani, Masoud; Montazeri, Mona; Farrokhi-Asl, Hamed; Rafiei, Hamed (2016) : A multi-objective genetic algorithm for a mixed-model assembly U-line balancing type-I problem considering human-related issues, training, and learning, Journal of Industrial Engineering International, ISSN 2251-712X, Springer, Heidelberg, Vol. 12, pp. 485-497, https://doi.org/10.1007/s40092-016-0158-6 This Version is available at: https://hdl.handle.net/10419/157515 Standard-Nutzungsbedingungen: Die Dokumente auf EconStor dürfen zu eigenen wissenschaftlichen Zwecken und zum Privatgebrauch gespeichert und kopiert werden. Sie dürfen die Dokumente nicht für öffentliche oder kommerzielle Zwecke vervielfältigen, öffentlich ausstellen, öffentlich zugänglich machen, vertreiben oder anderweitig nutzen. Sofern die Verfasser die Dokumente unter Open-Content-Lizenzen (insbesondere CC-Lizenzen) zur Verfügung gestellt haben sollten, gelten abweichend von diesen Nutzungsbedingungen die in der dort genannten Lizenz gewährten Nutzungsrechte. Terms of use: Documents in EconStor may be saved and copied for your personal and scholarly purposes. You are not to copy documents for public or commercial purposes, to exhibit the documents publicly, to make them publicly available on the internet, or to distribute or otherwise use the documents in public. If the documents have been made available under an Open Content Licence (especially Creative Commons Licences), you may exercise further usage rights as specified in the indicated licence. http://creativecommons.org/licenses/by/4.0/ ORIGINAL RESEARCH A multi-objective genetic algorithm for a mixed-model assembly U-line balancing type-I problem considering human-related issues, training, and learning Masoud Rabbani 1 •Mona Montazeri 1 •Hamed Farrokhi-Asl 2 •Hamed Rafiei 1 Received: 29 March 2015 / Accepted: 7 June 2016 / Published online: 1 July 2016 The Author(s) 2016. This article is published with open access at Springerlink.com Abstract Mixed-model assembly lines are increasingly accepted in many industrial environments to meet the growing trend of greater product variability, diversification of customer demands, and shorter life cycles. In this research, a new mathematical model is presented considering balancing a mixed-model U-line and human-related issues, simultaneously. The objective function consists of two separate components. The first part of the objective function is related to balance problem. In this part, objective functions are minimizing the cycle time, minimizing the number of workstations, and maximizing the line efficiencies. The second part is related to human issues and consists of hiring cost, firing cost, training cost, and salary. To solve the presented model, two well-known multi-ob- jective evolutionary algorithms, namely non-dominated sorting genetic algorithm and multi-objective particle swarm optimization, have been used. A simple solution representation is provided in this paper to encode the solutions. Finally, the computational results are compared and analyzed. Keywords Mixed-model assembly lines U-shaped assembly lines Learning and training effect  Human-related issues Multi-objective Introduction An assembly line is a group of successive workstations, joined by a material handling system. In each workstation, a set of tasks are carried out using a predefined assembly process, in which the time required to carry out each task and a set of priority relations which determines the order of the tasks are defined. The current market is severely competitive and consumer-centric with high variety in demands. As a result of high cost to establish and maintain an assembly line, the manufacturers produce one model with various features or several different models on a single assembly line. In situations like this, the mixed-model assembly line balancing problem arises to smooth the production and decreases the cost. Mixed-model Assembly Line (MMAL) is a kind of production line, where a set of similar models of a product are assembled to respond to the diversity of customer’s demands. There are two types of assembly line balancing problems. The purpose of type-I problems are minimizing the number of workstations. In this problem, the required production rate, assembly tasks, tasks times, and precedence requirements will be given. In type-II problems, the goal is to minimize the cycle time and maximize the production rate with fixed number of workstations or production employees. This study is mainly focused on the type-I problem, which wants to minimize the number of workstations. U-type line balancing was first invented by Miltenburg and Wijngaard (1994). The U-type assembly line is an attractive substitute for assembly production systems from the time operators became multi-skilled by performing tasks defined on different parts of assembly line (Go ¨kc¸en et al. 2005). The advantage of the U-type assembly line is the flexibility that it offers to choose an appropriate number of operators to satisfy demand changes (Aigbedo and Monden 1997). &Masoud Rabbani [email protected] 1 School of Industrial Engineering, College of Engineering, University of Tehran, Tehran, Iran 2 School of Industrial Engineering, Iran University of Science and Technology, Tehran, Iran 123 J Ind Eng Int (2016) 12:485–497 DOI 10.1007/s40092-016-0158-6 Learning effect is another important factor at assembly lines in the time of the new product lunch, or start of production (Baloff 1971). The length of the learning stage has becomean important performance indicator for a firm because of some common topics, such as shortened product life cycles, high innovation rates, and, therefore, more frequent product launches. Learning effect has to be considered in firms, because shorter learning stages enable firms to increase sales and, as a result, achieve more profits with the highest revenues, by the time, the new product reaches the market. Learning effects may occur by a highly repetitive execution of certain tasks. ‘‘A worker learns as he works; and the more often, he repeats an operation’’. Andress (1954) mentioned, learning effects at assembly lines and overall for repetitive operations. According to aircraft construction, Wright (2012)described learning effects at assembly lines and overall for repetitive operations. He figured out that by making the cumulated output double, average construction costs per unit sunk by about 20 %. This observation was formalized as an inversely proportional relationship between unit costs and cumulated output called learning curve. After that, for assembly lines in different industries, the presence of significant learning effects was confirmed. Basically, in mixed-model U-shaped assembly lines, workers are capable of operating several tasks. As Park (1991) said, training, the process by which workers become multi-skilled, has been recognized as a tool for boosting production flexibility. The minimum introduction of worker cross-training has the most significant improvement from no cross-training, and the subsequent increase of the cross-training has a diminishing return. In this research for the first time, a new model is presented considering both line balancing and worker assignment simultaneously, considering human-related issues. Two meta-heuristic algorithms [i.e., multi-objective particle swarm optimization (MOPSO) and non-dominated sorting genetic algorithm (NSGA-II)] are used to solve the proposed bi-objective problem, and a simple method is applied to represent solutions. The rest of the paper is organized as follows: in ‘‘Literature review’’, the relevant literature is reviewed. In ‘‘Problem description’’, the bi-objective problem, the objective function, and a mathematical model are presented. The methodology is described in ‘‘Methodology’’, and the illustrative examples are presented in this section. In ‘‘Parameters tuning’’, comparisons and discussion are brought. The study is finally ended by conclusions and future research in ‘‘Conclusion’’. Literature review The existing competitive and consumer-centric market and the observed trend of diversification of customer demands and high fluctuations is an important subject that is worth studying. Firms should improve their performance for dealing with these pressures to meet the customers demand within a short delivery time and with the lowest possible cost. Mixed-model assembly lines are one of the most relevant production environments that deal with these problem. The assembly line balancing problem encompasses assigning tasks to an ordered sequence of stations, such that precedence relations among tasks should not be violated (Erel and Sarin 1998). A mixed-model assembly line is assembly line, in which some similar product type with some insignificant difference is assembled. Many attempts have been made to solve the assembly line balancing (ALB) problems using the exact solution methods, heuristics, and meta-heuristic approaches. Some comprehensive reviews of such studies have been done (Becker and Scholl 2006; Erel and Sarin 1998). Some researches solved the assembly line balancing problem using a ranked positional weight technique (Helgeson and Birnie 1961). Monden (2011) was concerned with the sequencing of assembly lines, such as considering the stability of parts usage rates. Kim et al. (2009) presented a mathematical formulation and a genetic algorithm for the ALB-II problem. Some practitioners presented a formal ALB-I problem, and they also developed a branch-and-bound algorithm to solve the problem (Wu et al. 2008). Erel and Gokcen (1999) proposed a study that was concerned with minimizing the task time for different models considering precedence constraints using shortest route formulation. A binary integer formulation for the mixed-model assembly line balancing problem is developed by Go ¨kcen and Erel (1998). In another work, Gokcen and Erel (1997) extended a goal programming approach which was previously developed by Thomopoulos (1967), using a combined precedence diagram. Vilarinho and Simaria (2002) develop a two-stage heuristic method for balancing mixed-model assembly lines. The application of genetic algorithms (GA) for assembly line balancing has widely been considered in many studies. A genetic algorithm for type-II problems was presented by Anderson and Ferris (1994), and Leu et al. (1994) presented a GA-based approach to solve type-I problems with multiple objectives. Kim et al. (1996) presented a genetic algorithm for work load smoothing. In another study, a hybrid genetic algorithm approach to the assembly line planning problem was developed (Chen et al. 2002). There are only a few studies which use more than one meta-heuristic approach to solve their problem, but in this study, two meta-heuristic algorithms (i.e., MOPSO and NSGA-II) are used to solve the proposed bi-objective problem. Many practitioners studied the mixed-model straight line assembly line balancing problem which has been reported in the literature (Erel and Gokcen 1999; McMullen and Frazier 1998; Simaria and Vilarinho 2004,2009; 486 J Ind Eng Int (2016) 12:485–497 123 Thomopoulos 1967; Vilarinho and Simaria 2002). Simaria and Vilarinho (2009) proposed a mathematical programming model to formally describe the MMALB problem presenting an ant colony optimization algorithm. One of the effective factors for realizing the objectives of lean manufacturing is workforce planning. Several options of alternative production planning that can be applied for dealing with changing demand patterns, considering use of variable workforce, overtime, seasonal inventory, and planned backlogs have been developed by Hax and Candea (1984). Several classical LP models combining the production, manpower, and inventory-related trade-offs in each of the options mentioned above have been presented (Bhatnagar et al. 2003). Just-in-time (JIT) is able to adjust to changes in the external environment of the firm, because of several reasons, including efficient facility layouts and multi-functional workers (Monden 2011; Schonberger 1983). Japanese companies are operating with very low level of inventory and recognizing a high level of productivity using the just-in-time (JIT) manufacturing system which has the goal of continuously reducing and ultimately removing all forms of wastes (O ¯no 1988). The replacement of the traditional straight lines with U-shaped production lines is one of the most important changes resulting from JIT implementation (Chiang and Urban 2006). Reducing the work in process inventory and wasted operator’s movement, labor productivity improvement, material handling improvement, zero-defects campaign’s implementation, and higher flexibility in workforce planning in the face of changing demand patterns (Monden 2011) are the main benefits of the U-line as compared to a straight line. (In some reference, it is shown that one of the best applicable types of line is U-shape line and they illustrate that the benefits are impressive. The main characteristics in a U-shaped line are (Miltenburg and Wijngaard 1994): the U-line arranges machines around a U-shaped line in an operators work inside the U-line; U-lines are rebalanced periodically when production requirements change; the operators must be multi-skilled and versatile to do several different processes; it requires operators to walking, when setup times are negligible; U-lines are operated as mixedmodel lines, where each station is able to produce any product in any cycle; when setup times are larger, multiple U-lines are formed and dedicated to different products. Miltenburg and Wijngaard (1994) have a comprehensive article in the subject of U-shaped production line. In his article, the benefits of U-shape line were mentioned, and by some statistic information, they are proved for all). There are several studies on line balancing problems. Most of them assumed that the time of tasks for repetition tasks is independent from learning of workers. A few researchers have examined the learning effect on assembly line balancing problems (Chakravarty and Shtub 1988; Cohen and Ezey Dar-El 1998; Cohen et al. 2006). Learning can play a considerable role in manufacturing environments and there are many empirical studies that have proven learning effects (Cochran 1960; Yelle 1979). Learning occurs on the part of workers directly involved into manufacturing of the product (Andress 1954). The first model of Wright (2012) describes the learning rate as a relative decrease in average costs per product unit over the whole history of production. The second learning curve model, called Crawford or Stanford model (Yelle 1979), introduces the learning rate as a relative decline in the marginal costs, i.e., costs required to produce the last product unit. It is being observed that learning is present only in the initial production state, i.e., after a while task times converge to steady-state task (Table 1). A brief review of the related literature and contributions of this study is presented in Table 2. Problem description In this study, the focus is on minimizing the number of stations to achieve an optimum balance; therefore, the idle time should be minimized and the efficiency of the line should be enhanced. These goals may be achieved by smoothing the amount of workload and maximizing the equalization of the workload among stations. It was assumed that training, which is done to promote workers to upper levels, is performed between periods and it takes zero time. Workers are classified into four types based on their skill levels. The level of each work station indicates types of workers allowed to work at that station. Each worker has exactly one skill and exactly belongs to one skill level. Workers with skill level 4 can work on task levels 1, 2, 3, and 4. Workers with skill level 3 can work on task levels 1, 2, 3, and so on. In each period, workers can be trained to improve their working abilities to operate Table 1 Worker skills promotion possibility Skill level 1 Skill level 2 Skill level 3 Skill level 4 Skill level 1 –*** Skill level 2 ––** Skill level 3 –––* Skill level 4 –––– J Ind Eng Int (2016) 12:485–497 487 123 other task levels. The initial number of workers with skill level O in the beginning of the planning horizon is known. Levels of tasks are known, and the level of each station is equal to the maximum level of tasks which are assigned to it. Assumptions •Parallel stations are not allowed. •Operator walking time is ignored. •All parameters in the model are assumed to be deterministic. •There is no uncertainty. •Each task must be assigned to exactly one station. •All predecessors or successors of a task have already been assigned to a station (the precedence constraint. •The total time of the tasks assigned to each station, (i.e., the station time), may not exceed the cycle time (the cycle time constraint). •Salary is merely dependent on worker’s skill level and not depending on machine levels. •All of the machine types which need the same skill levels assumed to be similar in worker assignment. •Cost of hiring and firing are given, and they merely depend on skill levels. •Each task needs just one worker. •Training, which is done to promote workers to upper levels, is performed between periods and it takes zero time. •The productivity of experienced workers is assumed to be equal to 100 %. •The productivity of newly trained workers is assumed to be fewer than that of experienced ones, and it depends on the skill level to which they are trained. •Productivity of newly hired workers is assumed to be fewer than that of experienced ones, and it depends on the skill level for which they are hired. •Cost of training from one skill level to another is given, and it depends on both skill levels. Objective functions Minimizing the number of stations which is equivalent to the minimization of the idle time related to the line is one of the most important objectives in this article. Each model’s idle time is multiplied by the corresponding proportion (q0j). Computation of total weighted idle time (WIT) is shown below (Manavizadeh et al. 2013; Simaria and Vilarinho 2009). Minimize WIT ¼X M j¼1 q0 jX R r¼1 CX I i¼1 xir tir ! :ð1Þ Table 2 Overview of the related literature and contributions of this study Study Human–related issues (hiring, firing, and salary) Training and learning ALB-I ALB-II Parallel stations assembly line balancing U-l type assembly line balancing Simple assembly line balancing Two-sided assembly line balancing Mixed-model assembly line balancing Method This study * * * * * NSGA-II MOPSO Yuan et al. (2015) * * HBMO and SA Kucukkoc and Zhang (2014) * * * Agent-based ACO enhanced heuristics and model sequencing agent Manavizadeh et al. (2013) ** * * SA O ¨zcan and Toklu (2010) **GA Simaria and Vilarinho (2009) ** *SA Kim et al. (2009)* *GA Wu et al. (2008)* * Branch-and-bound Aryanezhad et al. (2009) ** Model numerical examples Simaria and Vilarinho (2004) * * * * Genetic algorithm 488 J Ind Eng Int (2016) 12:485–497 123 By balancing the workloads between stations, the idle time will be distributed across the workstations as equally as possible for each model. The workload balance between workstations will be computed by function B b . Therefore, the objective would be minimization of B b , as shown below (Simaria and Vilarinho 2009): Minimize Bb¼v v1X R r¼1 idr WIT 1 v  ð2Þ where id r is the idle time of workstation r: idr¼X M j¼1 q0 jidrj:ð3Þ The value of function B b is within the value range of [0,1]. In worst case, where the average idle time of the line is equal to the idle time of one of the workstations, the value equals 1, and in optimal case, equals zero when it is equally distributed among all workstations in the line actually it. By minimizing B w , the optimal value for B w is calculated as shown below (Simaria and Vilarinho 2009). Minimization Bw¼M2 vM 21ðÞ XX q0 jidrj idrj 1 M2  : ð4Þ The value of B w is within the value range of [0,1]. In worst case, when only model attributes to the idle time of each workstation, it equals 1, and when all models attribute equally to the idle time at each workstation it equals zero (Simaria and Vilarinho 2009). The value of WIT is different from one problem to another due to their dependence on the cycle time and task processing of each specific problem, whereas the function B b and B w are always within the value range of [0,1]. An alternative measurement, which is always within a fix range of values, is the weighted line efficiency (WE) (Simaria and Vilarinho 2009). This value varies between 0 and 1 is a direct indication of the efficiency of the line; 1 being the optimal value which indicates no idle time is found. The WE in an objective function computed as follows (Simaria and Vilarinho 2009): WE ¼Xq0 jPI i¼1tij vc ! :ð5Þ Another important objective is ti distribute tasks among workstations in a balanced fashion based on the job processing time. To achieve this goal, the difference between processing time of each model in each station and the average processing time for each model should be minimized. The formula is given below: WI ¼X R r¼1X M j¼1X I i¼1 tij xij meanj ð6Þ where t ij is the processing time of task irelated to model j, and x ir is equal 1 if task iassigned to station rand mean j is the average processing time workload needed for model j(Simaria and Vilarinho 2009). meanj¼P I i¼1 tij I:ð7Þ In Z 2 , we want to minimize all costs related to operators considering: Hiring cost: X T s¼1X 4 o¼1X MS k¼1 ho;suo;s;kð8Þ Training cost: X T s¼1X 4 o¼1X o0X MS k¼1X MS k0¼1 co;o0 ;sUXo0 ;o;s;k0 ;k ð9Þ Salary: X T s¼1X 4 o¼1X MS k¼1 so;sEo;s;kð10Þ Firing cost: X T s¼1X 4 o¼1X MS k¼1 fo;sWo;s;k:ð11Þ Mathematical model Parameters i,bIndex of task RMaximum number of stations r, r0Index of station JModel (product) {1,…,M} sIndex of period OWork skills category {1, 2, 3, 4} k,k0Index for station levels {1, …,MS=4} MNumber of models VNumber of operators ITotal number of tasks in the combined precedence diagram, (i=1, 2, 3, …,I) MS Number of station levels DThe vector presenting the total demand for each model, D={D 1 ,D 2 ,…,Dm} q0The overall proportion of the number of units of model j P ib Showing the precedence relationship between task b and i. Equal 1 if task bis the precedence for task i su ib Showing the succeeding relationship between task b and i. Equal 1 if task bis a successor for task i J Ind Eng Int (2016) 12:485–497 489 123 o ib A zero–one variable which determines whether or not constraints 2 or constraint 3 is satisfied CCycle time PTotal time in the planning horizon id r Idle time of station r D js Demand of model jin period s t ij Processing time of task iof model j w o Number of workers of skill category o w s o Number of workers of skill category oworking in period s pt s Regular time rate for workers during period s ot s Overtime rate for workers during period s h0Total working hours in a period h0Minimum overtime work for operators h o,s Cost of hiring of a worker with skill level oin period s s o,s Salary of each o-level worker in period s f o,s Firing cost of each o-level worker fired in period s Co,o00,s Training cost of each o-level worker trained for skill level o0in period s a o Productivity of each newly o-level worker hired in period s0\a o \1 b o,o 0Training productivity of o-level worker trained for skill level o00\b o,o 0 \1 a ro Equals 1 if workers of skill category ocan work at processing stage rand zero Decision variables x ir Equals 1 if task iis assigned to station rand equal 0 otherwise y0 r Equals 1 if workstation ris used for assembly and 0 otherwise x0 rs Total number of overtime hours done by workers at station rin period s x rs o Equals 1 if worker from skills category ois allocated to station rin period s U o,s,k Number of o-level workers who are hired and assigned to station level kin period s E o,s,k Number of existing o-level workers who are assigned to station level kin period s UX o’,o,s,k, Number of o0–level workers who were assigned to task level kin period s-1 and now are trained to skill level oand assigned to task level k0in period s UG o 0 ,o Equals 1 if training from skill level o0to skill level ois possible and 0 otherwise meanj¼PI i¼1tij I z1¼^ c1 PM j1q0 jPI i¼1tij vc  0 B B @1 C C AþUWIT þc X R r¼1X M j¼1X I i¼1 tij xij meanj  ! þ2V: Subject to: X R r¼1 xir ¼1ð12Þ X r1 xir1xbr Moib 8i;b;r;r11;pib ¼1 ð13Þ X r1 xir1xib M1oib ðÞ8i;b;rr11su ib ¼1 ð14Þ tjr ¼X I i¼1 xir Max tij 8j;r;sð15Þ X I i¼1X M j¼1 xir Max tij C8rð16Þ C¼p PM j¼1Dj ð17Þ q0 j¼Dj PM j¼1Dj ð18Þ v¼X R r¼1 y0 rð19Þ WIT ¼X M j¼1 q0 jX R r¼1 CX I i¼1 xir tir !ð20Þ idrj ¼CX I i¼1 xir tir 8r;jð21Þ idr¼X M j¼1 q0 jidrj 8rð22Þ X I i¼1 xir y0 r8r:ð23Þ 490 J Ind Eng Int (2016) 12:485–497 123 Minimizing Z2¼X R r¼1X 4 o¼1X T s¼1 xo rs pts  þx0 rs ots  þX T s¼1X 4 o¼1X MS k¼1 ho;suo;s;k þX T s¼1X 4 o¼1X o0X MS k¼1X MS k0¼1 co;o0 ;sUXo0 ;o;s;k0 ;k þX T s¼1X 4 o¼1X MS k¼1 so;sEo;s;k þX T s¼1X 4 o¼1X MS k¼1 fo;sWo;s;k: Subject to: X 4 o¼1 xo rs aro ¼18r;sð24Þ X R r¼1 xo rs aro ¼wo s8o;sð25Þ wo swo8o;sð26Þ tj¼X I i¼1 tij 8jð27Þ X T s¼1 xo rs  h0þx0 rs  X T s¼1X M j¼1 Djs tj 8rð28Þ v¼X R r¼1X 4 o¼1 xo rs ð29Þ Eo;s;k¼Eo;s1;kþUo;s;kWo;s;k þX 4 s¼1X MS k0¼1 UXo0 ;o;s;k0 ;kUXo;o0 ;s;k0 ;k  8o;s;k ð30Þ AX 4 o¼1"Eo;s1;kþaoUo;s;kWo;s;k þX 4 e¼1X MS k0¼1 bo0 ;oUXo0 ;o;s;k0 ;kUXo;o0 ;s;k0 ;k  #X R r¼1X 4 o¼1 xo rs ð31Þ Wo;s;kRaro 8o;s;kð32Þ Uo;s;kRaro 8o;s;kð33Þ UXo;o0;s;k0 ;kRaro 8o0 ;o;s;k;k0ð34Þ UXo0;o;s;k0 ;kRaro 8o0 ;o;s;k;k0ð35Þ UXo0;o;s;k0 ;kRUGo0;oð36Þ X O¼4 o0X MS k0 UXo;o0;s;k0 ;kþWo;s;kEo;s1;k8o;s;kð37Þ X O¼4 o0¼1X MS k¼1 UXo0 ;o;s;k0 ;kRyo;s;k08o;s;k0yo;s;k0¼0;1½ ð38Þ Wo;s;k0R1yo;s;k0  8o;s;k:ð39Þ Methodology Proposed model in this paper is multi-objective, so the methods for solving the problem are NSGA-II and MOPSO. Rabbani et al. (2016a,b) applied these two algorithms for solving a mixed-model assembly line problem, and the results obtained by these two algorithms were compared to each other. NSGA-II is a popular nondomination-based genetic algorithm for multi-objective optimization. It is a very effective algorithm but has been generally criticized for its computational complexity, lack of elitism, and for choosing the optimal parameter value for sharing parameter (Rabbani et al. 2016a,b). Kusiak and Wei (2012) introduced MOPSO for optimizing continuous non-linear functions, Particle Swarm Optimization (PSO) defined a new era in Swarm Intelligence (SI). PSO is a population-based method for optimization. The population of the potential solution is called as swarm and each individual in the swarm is defined as particle. PSO is motivated by social behavior of birds flocking or fish schooling Solutions are represented by particles in the search space. The particles fly in the swarm to search their best solution based on experience of their own and the other particles of the same swarm. PSO started to hold the grip amongst many researchers and became the most popular SI technique soon after getting introduced, but due to its limitation of optimization only of single objective, a new concept Multi-Objective PSO (MOPSO) was introduced, by which optimization can be performed for more than one conflicting objectives, simultaneously. Coello et al. (2002) described the advantages of using MOPSO in solving multi-objective optimization problem rather than the single objective version of the algorithm. Representation of solutions The chromosome is a string of length Iwhich shows the task numbers, where each element represents a task and the value of each element represents the workstations to which the corresponding task is assigned. The maximum number J Ind Eng Int (2016) 12:485–497 491 123 of stations is equal to total number of tasks. For example, for 16 tasks, 9 workstations will be created. In this research, individuals in the initial population are all randomly generated. While a heuristic procedure can provide good initial solutions, it can cause the solutions to be biased. Illustrative example In this section, 5 small-size and 5 large-scale problems are implemented to compare the performance of algorithms with each other in various size problems. Parameters of problems were generated based on Table 3. In this paper, the workers assignment is based on their skill level. Workers with skill level 4 can work on task levels 1, 2, 3, and 4. Workers with skill level 3 can work on task levels 1, 2, 3, and so on. The problem with five tasks is as follows: The precedence diagram of five task problems is shown in Tables 4,5,6,7,8,9. The results from NSGA-II algorithm are shown below: This table shows that task number 1 is assigned to workstation number 4, task 2 and task 3 are assigned to workstation number 3, task number 4 is assigned to workstation number 2, and task number 5 is assigned to workstation number 1. Training should happen according to Table 10: Parameters tuning The efficiency of the meta-heuristic algorithms in finding better solutions in less run time is considerably dependent on their parameters. To setting the MOPSO and NSGA-II parameters, design of experiment (DOE) using Taguchi approach is used in the paper. The performance of NSGA-II is influenced by four parameters, including population size (Np), maximum number of generations (Max Iteration), mutation rate (Pm), and crossover rate (Pc). MOPSO parameters consist of population size (Np), maximum number of iterations (Max Iteration), inertia weight (w), repository size (Nr), personal learning coefficient (c1), and global learning coefficient (c2). After specifying levels for each parameter (factor), design of experiment is performed using the Minitab software to set these two groups of parameters (Figs. 1,2,3). Parameters tuning for both algorithms are done according to the results of large-sized problem (Table 11). The consequences of Taguchi method in tuning of parameters are shown in Figs. 4and 5.In addition, the results are summarized in Table 12. Comparative results Comparison metrics: It is common to compare the performance of the multi-objective algorithms’ performance by means of some specific comparison metrics; to compare proposed algorithms with each other, three comparison metrics are employed (Rabbani et al. 2016a,b). 1. Number of Pareto solutions (NPS): The quantity of non-dominated solutions that every algorithm can discover. 2. Spacing metrics (SM): This kind of metric provides us details about the uniformity of the distribution of the solutions obtained by the way of each algorithm. This metrics are computed as follows: SM ¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 N1X n i¼1 did  2 sð40Þ Table 3 Test problem generation Parameters Value Parameters Value Demand U(5, 10) Hiring cost U(1500, 2000) Processing time (2, 5) Firing cost U(1500, 2000) Training cost U(50, 150) Salary U(100, 500) Table 4 Initial number of workers with skill level 1 in the beginning of the planning horizon Skill level 1 2 3 4 Initial number of workers 5110 Table 5 Level of tasks Task 1 Task 2 Task 3 Task 4 Task 5 Skill level 1 Skill level 2 Skill level 3 * * Skill level 4 * * * Signed cells means that worker with skill level ocan work the task number jconsidering tasks level, and in addition, the training cost from skill level Oto skill level O0in period s is shown in Table 6 Table 6 Cost of training from skill level Oto skill level O0in periods From skill Period/ to skill Skill level 2 Skill level 3 Skill level 4 Skill level 1 1 118 64 127 2 125 69 112 3 118 83 117 Skill level 2 1 97 130 29565 3 106 66 Skill level 3 1 124 286 3 146 492 J Ind Eng Int (2016) 12:485–497 123