Research Article  Open Access
Multiobjective Fuzzy Mixed Assembly Line Sequencing Optimization Model
Abstract
It can be deduced from previous studies that there exists a research gap in assembly line sequencing optimization model for mixedmodel production lines. In particular, there is a lack of studies which focus on the integration between job shop and assembly lines using fuzzy techniques. Hence, this paper is aimed at addressing the multiobjective mixedmodel assembly line sequencing problem by integrating job shop and assembly production lines for factories with modular layouts. The primary goal is to minimize the makespan, setup time, and cost simultaneously in mixedmodel assembly lines. Such conflicting goals arise when switching between different products. A genetic algorithm (GA) approach is used to solve this problem, in which trapezoidal fuzzy numbers are implemented for variables such as operation and travelling time in order to generate results with higher accuracy and representative of realcase data.
1. Introduction
It is known that an efficient and effective production line sequencing and scheduling are required for manufacturing industries due to the fluctuations in market demand and increasing competition on the global scale. Hence, manufacturers are gearing towards optimizing the design of manufacturing systems in the shortest time possible in order to enhance competitiveness. The production lines of a multistage process can be generally classified as jobshop (JS) and flowshop (FS) line processing, in which each type complies with a set of sequential operations. Production sequential operation is established in the factories in which the materials undergo a refining process in order to produce a product that is suitable for onward consumption and the components are assembled to make finished articles. Raw materials require a sequence of treatments to render them useful. Development of optimization methods for solving JS and FS problems effectively is of utmost importance to the apparel industry and other manufacturing industries which require similar assembly operations at large. JS and FS involve the assignment of a set of tasks to workstations (machines) in a predefined sequence while optimizing one or more objectives without violating the restrictions imposed by the production line [1]. Most of the problems involving line balancing and sequencing problem in JS and FS fall within the NPhard class of combinatorial optimization problems, which results in a critical need to develop efficient algorithms to attain optimum solutions [1–6]. It can be extremely timeconsuming for optimumseeking methods to obtain optimum solutions within a vast search space [2].
Efficient and effective JS and FS sequencing and scheduling is required for recent manufacturing processes in order to achieve higher customer satisfaction and shorter production lead times and enhance competitiveness in the escalating volatile market demand. This emphasizes the significance of justintime (JIT) and agile manufacturing in the modern era. To cope with such circumstances, JS and FS assembly lines are shifting towards mixedmodel assembly lines, in which multiple models are produced on a constant basis. Therefore, sequencing with multiobjective goals is a critical operation scheduling that directly affects the efficiency of mixedmodel production assembly lines for factories with modular layouts. It will be noted that modular layouts are layouts which integrate JS and FS. Owing to the importance of this issue, this study is aimed at developing a multiobjective mixedproduction line sequencing model by implementing genetic algorithms (GAs) and fuzzy logic.
2. MixedModel Assembly Line Sequencing
Mixedmodel assembly lines (MMAL) have been widely used by manufacturers and they play a key role in producing a variety of products. Products with similar characteristics are assembled with different processing times on the same assembly line at very low cost [7–9]. This, in turn, requires the implementation of costefficient and flexible production systems. MMAL reduce setup operations to an extent that various models from a common base product can be manufactured in intermixed sequences. However, the observed diversity of MMAL makes thorough sequence planning essential in order to exploit the benefits of assembly line production. These benefits include diversified small lot production and prompt response to sudden changes in demand for models without large inventories [10]. Mixedmodel sequencing (MMS) aims at avoiding or minimizing sequencedependent work overload based on detailed scheduling which explicitly accounts for operation time, worker movement, station borders, and other operational characteristics of the line [10].
Solving the MMS problem involves determining a production sequence for multiple products along a modular layout. In order to reach this goal, manufacturers aim to achieve optimum scheduling times with minimum cost. Furthermore, it is impractical for manufacturers to constantly alter their machine layouts when producing a new product. For companies with fixed factory machine layouts, optimization of production lines is solely based on altering the sequence of production lines. However, this approach poses the following challenges that need to be addressed:(1)production of various products using existing facilities such as the number and location of machines, travelling time between the machines, and limited number of workers for different products;(2)assignment of parts from different products to the machines, in which the optimum sequence needs to be considered in order to minimize scheduling time. Problems also arise when changing the setup number of the machines;(3)assignment of parts to other parts during product assembly, in which the parts that are ready for assembly need to be considered in order to minimize the scheduling time required to assemble the product;(4)identification of the best minimum setup number and setup cost sequencing, in which the above mentioned issues need to be considered simultaneously.
It can be seen that identifying the best scheduling and sequencing production line is a challenging task, particularly for the production of a mixedmodel from a large number of providers with multiple objectives. MMS is a NPhard problem which requires an effective model. A unique and stable model is needed to facilitate MMS in the production of mixed products in a manufacturing environment with multiple parts, machines, products, and assemblies in order to minimize scheduling time, idle time of the machines, setup number, and setup cost as well as to maximize the number of products and assembly of various products. An effective model provides companies with a means to reschedule their production lines based on the changing demands of a flexible market. The effective utilization of MMAL requires solving two problems in a sequential manner, as follows.(1)Levelling workloads for stations on the line: this is known as line design and balancing.(2)Maintaining a constant rate of usage for each part used on the production line: this involves determining the production sequence for various models.
The first goal involves sequencing the mixedmodels in order to achieve a balanced workload at each assembly station, whereas the second goal involves sequencing the mixedmodels in order to minimize variations in the usage of various parts and components over time [8, 9, 11]. A number of MMS studies attempted at solving mixedmodel sequencing problems using sequencing procedures that will optimize various system measures such as throughput, scheduling time, number of stations, idle time, flow time, line length, workinprocess, and raw material demand deviations [12–16].
Zhu et al. [17] studied the complexity of product variety of MMAL defined by the operators such as selection of parts, tools, and fixtures as well as assembly procedures in sequence planning. The aim was to reduce manufacturing complexity in a multiproduct, multistage, manual assembly line environment. The methodology was developed to determine optimum assembly sequences which will minimize system complexity. Rekiek et al. [18] and De Lit et al. [19] developed an approach which integrates product family design (including assembly sequences) with the assembly system so that multiple products can be assembled on the same line. A new alternative for mixedmodel assembly for lowvolume manufacturing environments was investigated by Heike et al. [7]. The efficient utilization of labour and/or machinery presents significant challenges. The models developed focus on evaluating and understanding different alternatives for mixedmodel assembly in the aerospace industry. These models include one linear and two nonlinear programs which are formulated to evaluate constant and variable scheduling time policies. The models were successfully implemented on an airplane assembly process characterized by lowvolume manufacturing. Bukchin et al. [20] proposed a mathematical formulation for MMAL designs based on the maketoorder production policy. Hence, the arrival sequence is randomly distributed according to demand proportions, rendering the balancing procedure an important aspect of the design process. An integrated machine tool selection and sequencing model was proposed by Moona et al. [21]. The model determines machine visiting sequences for all part types, such that the total production time for the production order is minimized and the workloads between machine tools are balanced. Numerical experiments were carried out using problems of various sizes in order to demonstrate the efficiency of the proposed GA on the integrated machine tool selection and sequencing problem.
Khan and Day [22] introduced a knowledgebased design methodology (KBDM) for automated and manual assembly lines, which can be applied equally well to single, multi, and mixedproduct assembly lines with either deterministic or stochastic operation times. Mendes et al. [23] developed simulation models for varying the levels of demand and line configurations and utilized a mixedmodel PC camera assembly line. The simulation models were used to compare flow time and resource utilization, which provide operational support and help finetune line configurations. Joly and Frein [24] investigated a set of vehicles within an industrial environment in order to minimize manufacturing costs by considering assembly and paint shop objectives. The approach was implemented to solve a problem with a monoobjective function. One heuristic (a progressive, constructionsequence algorithm) and three metaheuristics (simulated annealing, variable neighbourhood search, and an evolutionary algorithm) were described and compared.
The use of GA for mixedmodel sequencing line problems is discussed as follows. Ghosh and Gagnon [25] introduced a mathematical programming model and an iterative GAbased procedure for MALBP with parallel work stations, in which the goal was to maximize the production rate of the line for a predetermined number of operators. An intelligencebased GA was developed by Norozi et al. [26] in order to tackle the complexity of sequencing in parallel MMAL. A spreading and cutting sequencing (SCS) model using GA was implemented by Wong et al. [27] to solve the sequencing problem of a computerized fabric cutting system used in the garment industry. The job sequence of spreading and cutting operation can be optimized by reducing the completion time for daily operation of fabric spreading and cutting as well as improving the utilization of machines. Ponnambalam et al. [28] studied the performance of GAs for a multilevel MMAL sequencing problem. Both single objective and multiple objectives were considered. The main objective of the multiobjective GA was to minimize the total utility work levelling, parts usage, and total setup cost. Development of optimization methods for solving JS and FS that are involving assembly line balancing/sequencing in mixedmodel problems falls under the NPhard class problem. Optimization involves exploiting the capabilities of mathematical programming in order to minimize the overall makespan and setup objectives for sets of single lines [26]. Evolutionary computing is a research area within computer science for solving combinatorial optimization and complex problems based on the principles of generic populationbased heuristic techniques. In the advent of metaheuristic algorithms in recent years, numerous complex scheduling problems have been studied and solved using metaheuristic search techniques such as ant colony optimization (ACO), tabu search (TS), genetic algorithm (GA), and simulated annealing (SA). Metaheuristic algorithms are used to overcome the complexity of sequencing in assembly lines [26].
GA was introduced by David Goldberg based on the procedure of natural mechanism and natural genetics. The population is composed of a collection of chromosomes in which each string is encoded and the problem is solved as a finite length of genes [26]. Good gene characteristics produce new generations. After several evolutions, the genes become highly adaptive to a given environment [29]. GAs have been proven to be highly effective for achieving optimum or nearoptimum solutions in complex realworld optimization problems. However, GAs are limited by the fact that their performance is very sensitive to parameter settings [29]. GA design consists of two key steps, namely, genetic operations and parameter settings [29]. The genetic operations involve choosing a suitable selection method. The typical methods for selection are roulette, tournament, stochastic uniform, remainder, and uniform selection. Selection is the only operation of evolution and thus choosing a suitable selection method is important. Parameter settings involve setting the required parameters and variables for controlling the algorithms such as population size, number of generations, the number of selected candidates, crossover rate, and mutation rate [30].
2.1. Multiobjective MixedProduction Line Optimization
Multiobjective formulations are realistic models of numerous complex engineering optimization problems. In most reallife problems, the objectives often contradict one another and optimizing a particular solution with respect to a single objective may give undesirable results for other objectives. Being a populationbased approach, GA is well suited to solve multiobjective optimization problems [31]. A generic singleobjective GA can be modified to search for a set of multiple nondominated solutions in a single run. The ability of GAs to simultaneously search different regions of a solution space makes it possible to search a diverse set of solutions for difficult problems with nonconvex, discontinuous, and multimodal solution spaces. Multiobjective optimization problems can be found in various fields, such as product and process design, finance, aircraft design, oil and gas industry, automobile design, or wherever optimum decisions are needed, by considering the tradeoffs between multiple conflicting objectives. If more than one criterion is to be treated simultaneously, then it is a multiobjective optimization problem. Maximizing the profit and minimizing the cost of a product, maximizing the performance and minimizing fuel consumption of a vehicle, and minimizing the weight while maximizing the strength of a particular component are all examples of multiobjective optimization problems. If the multiobjective problem is well formed, there will not be a single solution that simultaneously minimizes each objective to its fullest. In each case, the aim is to search for a solution in which each objective has been optimized to an extent that if the objectives are optimized further, other objectives will suffer as a result. The search for an optimum solution and quantifying the extent to which this solution is superior compared to other such solutions (note that there will generally be many solutions) are the goal when setting up and solving a multiobjective optimization problem. Most of the realworld decision problems involve multiple conflicting objectives that need to be tackled while adhering to the various constraints. In multiobjective problems, there may not exist a solution which is the best with respect to all objectives. It will be highlighted that there will be a set of solutions which are superior to other solutions in the search space when all objectives are considered but are inferior to the solutions for one or more objectives. These solutions are called nondominated solutions [32].
The jobshop scheduling problem with sequencedependent setup times (SDSTJSP) is a generalization of the classical JSS, in which the setup operation of a machine is required when the machine switches between two jobs. The common goal of SDSTJSP is to minimize the completion time of the last job (makespan). SDSTJSP has been investigated by a number of researchers. A new multiobjective GA (MOGA) was introduced by Zacharia and Nearchou [33] for solving fuzzy assembly line balancing of SALBP2 comprising multiple objectives, that is, minimizing the fuzzy scheduling time, balancing the fuzzy delay time, and smoothing the fuzzy index of the line. McMullen [34] proposed a technique to generate production sequences, whereby the number of setups and scheduling flexibility are considered. These objectives are inversely correlated, which poses a challenge for optimizing both objectives simultaneously. An efficient frontier approach was exploited for problems where simultaneous maximization of flexibility and minimization of setups are desired. McMullen [35] presented a relatively new approach using ACO to address a twoobjective productionsequencing problem by stimulating the artificial intelligence agents of virtual ants to obtain desirable solutions for manufacturing logistics. The first objective involves minimizing the number of setups, while the second objective involves optimizing the stability of material usage rates. McMullen [36] developed a technique which addresses a JIT productionscheduling problem consisting of two objectives, namely, minimization of setups between product variants and optimization of schedule flexibility. These objectives are inversely correlated, which poses a challenge in attaining desirable results for both objectives. An efficient frontier approach was employed to address the problem, which yields desirable sequences for both objectives. The artificial neural network approach using Kohonen selforganizing map (SOM) was implemented to determine sequences which are desirable in terms of the number of setups and flexibility. McMullen and Tarasewich [37] introduced an efficient frontier approach to support sequencing decisions. A beamsearch heuristic approach was used to effectively generate efficient frontiers. RahimiVahed et al. [8] addressed a multiobjective sequencing problem using a mathematical formulation which considers three objectives simultaneously. These objectives are minimizing the total utility work, minimizing variations in the total production rate, and minimizing the total setup cost. The scatter search methodology was used to generate a set of solutions that approximates the nondominated frontiers where simultaneous minimization of the above mentioned objectives is desired. Javadi et al. [9] proposed a fuzzy goal programming approach to solve a multiobjective mixedmodel assembly line sequencing problem for a JIT production system. Three objectives are solved simultaneously as mentioned in [8]. These objectives comprise minimization of the total utility work and variations in the total production rate and total setup cost. The fuzzy goal programming based approach was proposed to solve the problem because of the conflicting objectives. Bard et al. [38] investigated the sequencing of MMAL, which was perceived as a multiobjective problem. They developed a method to minimize the overall line length and maintain a constant rate for part usage. The problem was solved using the weighted sum and tabu search methods. McMullen [39] focused on minimizing the number of setups and maintaining a constant rate for part usage, in which the multiobjective problem was solved using the TS method. McMullen and Frazier [40] developed an SA method for the model used by McMullen [39] as well as TS method. Ponnambalam et al. [28] compared the performance of selection mechanisms (i.e., Pareto stratumniche cubicle and selection based on scalar fitness function values). The objectives of their study were to minimize variations in part usage, total utility work, and setup cost. The results indicated that the GA which incorporates Pareto stratumniche cubicle exhibits superior performance compared to GAs with other selection mechanisms. Giard and Jeunet [41] presented an integer programming formulation for sequencing in mixedmodel assembly lines, whereby the number of temporarily hired utility workers and the numbers of sequencedependent setups are optimized simultaneously via a cost function. The objective function is a cost function involving two elements, namely, the cost associated with additional utility workers and setup cost.
2.2. Application of Fuzzy Numbers in Assembly Line Models
Fuzzy set theory accounts for the ambiguities involved when assigning processing and scheduling times as well as the uncertainty contained within such time variables [42]. The data of realworld problems are imprecise, vague, or uncertain, in which the input data can only be estimated within this uncertainty. The uncertainty may be represented by a fuzzy number. The membership function of fuzzy data represents the grade of satisfaction of the decisionmaker. The concept of fuzzy numbers was introduced to treat imprecise data such as the processing and travelling time of each job [43, 44]. The time variables were represented by triangular fuzzy numbers, which allow practitioners to account for ambiguities when assigning processing and scheduling times, while maintaining the variability of the stochastic environment [42].
Gen et al. [43] and Tsujimura et al. [44] proposed a new method for treating the data of realworld problems by introducing a fuzzy assembly line balancing based on GA. Fuzzy variables were implemented by Fonseca et al. [42] in order to solve assembly line balancing. Fuzzy set theory is a viable alternative method for modelling and solving stochastic assembly line balancing problems. The variability and uncertainty in assembly line balancing problems have been traditionally modelled by statistical distribution. A novel fuzzy extension of the simple assembly line balancing problem (SALBP) was proposed by Zacharia and Nearchou [33]. The fuzzy job processing times reflect the uncertainty, variability, and imprecision with which realworld production systems are afflicted. The job processing times are formulated by triangular fuzzy membership functions. A new MOGA was introduced for solving the fuzzy SALBP in order to minimize the fuzzy scheduling time, fuzzy balance delay time, and fuzzy smoothness index of the line.
2.3. Summary of Previous Studies and Findings
It can be observed from the literature review that the majority of works are focused on an assembly line balancing and only a few papers address mixedassembly line sequencing. Hence, mixedassembly line sequencing needs to be investigated in detail. Most papers focused on the number and order of mixedmodels, neglecting the assignment of a product’s subparts to their respective machines (job sequence). In addition, most papers overlook the designation and adaptation of parts with other parts as well as their final assembly. It is found that the production of mixedmodels is influenced by a number of criteria, which limits the applicability of the research results in real production line conditions. It will be highlighted that only one criterion was considered in previous works such as operation line, and therefore other factors such as the travelling time of the conveyor were neglected. Deterministic time has also been of interest in most studies. It will be highlighted that only assembly line balancing was carried out in previous works in order to minimize the prediction time for input data, which include the travelling and processing times for each available job. Although the results from mixedassembly line sequencing studies are applicable in real manufacturing environments, little is known about fuzzy mixedmodel assembly line sequencing. Hence, it is evident that there is a lack of studies which implements the concept of fuzzy time in order to minimize the prediction time of input data.
In general, the research papers can be classified into two groups, whereby the first group focuses solely on objective criteria, while the second group focuses on multiobjective investigations. A critical evaluation of previous works clarifies that addressing these objectives involves the development of various methods, in which MMAL sequencing is the ideal method for a single objective. Comparison of various GA methods reveals that multiobjective studies have not been investigated extensively, unlike singleobjective studies. Much effort has been made to intensify and accelerate the running of GA methods to achieve optimum results. Although mixedmodel assembly line sequencing is of prime importance, there is a lack of studies which focus on this topic. In this paper, a hybrid method is proposed, in which mixedmodel assembly line sequencing is integrated with the operating and travelling time in the form of fuzzy numbers of a multiobjective optimization problem. The hybrid method will provide manufacturers with an ideal and feasible alternative since mixedmodel assembly line sequencing is integrated into a software application. This has been ignored in previous works and much effort is needed to address the above issue.
3. Problem Definition and Assumptions
In a highly competitive market, companies try to produce various kinds of products to meet market demand. One manufacturing problem that is often associated with scheduling practices is mixedmodel production line sequencing as shown in Figure 1. Solving this problem involves determining the optimal sequence for jobs () which includes the assigned part () to machine (), considering each of the jobs with the availability of the machines () in jobshop. Moreover, with the need to solve the problem of determining the jobs () which includes assembly of the parts () in the assembly line for multiple products along a modular layout. The main solid arrows in Figure 1(a) represent the sequence to produce the parts, followed by a discrete part manufacturing assembly, leading to the final product. Likewise, the elapsed traveling time between machines for manufacturing a part is given by (). To satisfy the nowait restriction, the completion time of the operation () must be equal to the earliest time to start the operation. There must not be any waiting time in the processing of any consecutive operations of each of the jobs. Figure 1(b) represented by gray lines provides insight into problems involving the determining of the optimal number of setup tool changes sequencing () of jobs () assigned to “” machines, throughout the jobshop and assembly production line.
(a) The precedence problem in parts of the raw material process for producing a product
(b) The precedence problem of parts assigned to each machine based on changing the machine setup
The following assumptions are taken in the multiobjective evaluation: (i) all jobs are available at zero time; (ii) the conveyor (operator movement) moves at a constant speed. If job overlapping occurs, the remaining work will be accomplished by temporary operators; (iii) the job operation begins when the part enters the machine. Once the job is completed, the operator will move the part to the next machine; (iv) at a given time, only one job can be processed on one machine; (v) the operator is assigned to each selected part that is assigned to each machine; (vi) the position of the machines on the assembly line varies from one another depending on the user input, which is based on the travelling time; (vii) once operation has started on a machine, it cannot be interrupted before completion, either on or between the machines; (viii) the assembly line can process jobs for a product family, which is described by a joint priority matrix; (ix) the processing time varies for different jobs and these jobs are allocated to the same machine, no matter what the optimum processing time for each job might be; (x) the demand for all products as well as the sequence of the products entering the assembly line is predetermined; (xi) the completion time for all jobs is represented by a fuzzy number.
The goal is to minimize the makespan or total completion time and minimize the number of change tools setups for each machine in mixedmodel production line sequencing of all the jobs without waiting between successive operations of a job by using the fuzzy operation and travelling time which makes the solution more realistic.
4. Model Derivation
The objectives of the function value, for multiobjective mixedproduction assembly line sequencing, are categorized into two main objectives: minimizing the total makespan and minimizing the total setup number as discussed in the following.
4.1. Minimizing the Total MakeSpan
The step involves evaluating each chromosome based on the makespan by considering the product number and assembly sequence for a gene code (). The , , , and () are computed in which () are calculated as follows. is determined by checking the start time of the parts () entering the machines (), based on the sequence assigned to the machines. This is expressed by is determined by checking the start time machine’s availability () assigned to produce the parts (), as given by is then determined by checking the start time of the parts () entering the machines () based on the sequence of the part’s assembly () for each product () as given by (5) and (6). For each , where , note that the product number (()) and assembly sequence () are important. Equations (3) through (6) are given below.
(i) Consider
(ii) If , then
(iii) If when searching gene code , where , , go backwards one by one until is met as follows: The start time is obtained using (7), in which are determined from the previous steps. Consider Note. In the first run, if , then .
The makespan fitness function is then calculated for all chromosomes using
4.2. Minimizing the Setup Number
In many industries, sequencedependent setups are considered as an important item in assembly operations. The model considering sequencedependent setups is considered in this paper. The machine number and the setup number are defined for the element in the chromosome, where ; the total setup number at the beginning is . All the genes , where , in the chromosomes one by one will be checked and total setup number () based on the and has been determined for all as shown in
(i)
(ii)
5. Fuzzification of the MixedModel Line Sequencing
In order to develop a fuzzy version of the mixedmodel line sequencing, operation and travelling time had to represent a fuzzy number, which is applied for the processing and scheduling time of each robot and machine. The fuzzy sets that are most appropriate to represent these values are triangular fuzzy numbers (TFN) due to the nature of the processing and scheduling times. Rather than assigning a specific value for each time variable as in the case of deterministic time methods, TFN establish extreme points to represent the most and least likely values of the individual variables. The use of TFN intime modelling enables users to account for variability and ambiguity, similar to statistical distributions. However, TFN differ from statistical distributions such that they do not require prior knowledge or historical data when establishing their values. This is a major advantage of using TFN compared to statistics. TFN are achieved by defining the existing deterministic time as the most likely values for the new TFN. This implies that the deterministic value represents the fuzzy element in the set, with a membership value of 1. A value of 1 is added or subtracted from the most likely values in order to get the least likely values for the new TFN, with a membership value of 0. This same procedure is applied to adapt the deterministic scheduling times in the fuzzy domain [42]. No major changes are required when adapting these heuristics to handle fuzzy processing and scheduling times.
6. Solution Algorithm
In this study, a multiobjective optimization model is developed for a mixedmodel process, comprising the manufacture of parts, assignment of parts to their respective machines, and assembly of parts to form the final products. The multi objective mixed model involves minimizing the makespan and setup cost. The operation time and travelling time are represented by fuzzy numbers. It will be highlighted that the steps for GA modelling include model, constraints, fitness function, and satisfying multiple objectives within the domain of model restrictions, fitness function, and multiobjective cases. These steps are carried out for various types of applications and include input data, fuzzy variables, and initialization of parameters, evaluation, crossover rate, mutation rate, and termination. The input data are coded by considering the sequence of part production stages, sequence of product assembly, and sequence of part allocation to each machine during the initialization of parameters. The proposed model consists of input data and fuzzy variables. Care is taken to ensure that the setup and operation time for each machine as well as the travelling time of the parts between various machines are accounted for. The evaluation stage involves developing the multiobjective fitness functions. Nine rules are used to control the crossover rate during the crossover stage, whereas eight rules were used to control the mutation rate during the mutation stage. Two rules are designed for the termination stage in order to control termination and gain optimum results. Sorting, selection of tournament candidates, elitism, and new population size were performed based on the existing conditions.
The model is developed based on the works of [1, 5, 8, 28, 34, 45, 46]. A contribution of this paper is the development of a multiobjective fuzzy mixed production line sequencing optimization model using the genetic algorithm approach by integrating jobshop and assembly production lines for factories with modular layouts. The significance of this model is for those factories who want to produce the various kinds of products with fixed machine just by changing the sequencing of the products. The model helps the manager to sequence and schedule the production line easily and accurately by taking the market demand into consideration. The line sequencing is optimized based on the makespan and also the setup number of the production line. The flow chart of the proposed multiobjective fuzzy mixedproduction assembly line sequencing model is shown in Figure 2.
7. Development of Multiobjective Genetic Algorithm
The procedure used in the development of the model is described in detail as follows.
Step 1 (input data). The first step involves identifying the machine number () such as CNC, NC, and robot, as well as assigning the parts () to their respective machines and robots based on the production and assembly line sequence. Also, it needs to identify the fuzzy completion time of the operation () and fuzzy traveling time between machines in order to manufacture a part by ().
Step 2 (initialization of parameters). Initialization involves setting the parameters of the GA, creating the scores for the simulation, and creating the first generation of chromosomes based on the Notations section. The results of the model will depend on the difficulties encountered during production planning and sequencing. The general schematic for reading data in mixedproduction assembly line problems is presented in Figure 3. Figure 3 shows the number of finished products, the number of tasks involved in the total production line (), the order of parts in the assembly (), and the setup of each machine (). Figure 3 also shows the number of parts () and the order for manufacturing each part. For example, the user proceeds to the next item (part) once the order of the manufacturing process is completed for each part. represents the number of machines used for manufacturing parts and various products. The dashed arrows illustrate the coherence between these machines () during various phases of the production line. The fuzzy operation time () and fuzzy traveling time () are presented in Figure 3. The first column of Figure 3 (i.e., gene code ()) shows the coding of each row, which will be discussed later.
From the data shown in Figure 3, the number of existing gene codes in the first column represents empty chromosomes (), which are formed based on the job numbers. Once the chromosomes have been formed, the chromosomes will be filled with numbers via stochastic repeating. The numbers vary between 0 and the maximum number of genes in the chromosomes. The gene codes are filled up by random numbers, neglecting the encoding sequence of the produced chromosomes. Once the chromosomes have been filled with random numbers, the classified gene code numbers are ranked as a possible sequence which can be used by the GA. The steps are described briefly as follows.
Identification of the First of Each Part Type’s Gene Code. The first of each part type’s gene code is identified using (11), whereby each part has a sequence of gene codes with elements, where , , have the same part number. Consider
Classification Based on the Product’s Part Sequence. Each gene code number in the chromosome line is compared to the previous one. Three conditions are used when comparing the gene code numbers and are described below.(1)The first or existing gene code number fills up the new chromosome without any changes in the gene code number.(2)If the gene code number is the same as other gene code numbers in the chromosome, the chromosome is filled up with the addition of one number out of the existing gene code numbers.(3)If the gene code number is different from other gene code numbers within the chromosome, the existing gene code number fills up the chromosome without any changes in the gene code number.
Classification Based on the Assembly Sequence. Once the product’s part sequence has been classified, the chromosome is filled up using the assembly sequence number relating to a gene code () of the specific chromosome. Finally, the final gene code numbers are ranked by two sequence filters, which accounts for the assembly sequence code and product’s part sequence.
Step 3 (multiobjective evaluation). After initializing the population size, each chromosome must be evaluated during each generation of the selection process. The fitness function value involves evaluating each chromosome by minimizing the makespan and considering the product number and assembly sequence for a gene code () which are defined as (1) through (8). Moreover, the evaluation of each chromosome for the second objective function to minimize the setup number in mixedmodel line sequencing is defined as (9) through (10). Finally the total fitness values of the efficient frontiers are calculated based on these two objectives [35]. This is accomplished by looking up the score of each gene in the chromosome, by means of adding and averaging the scores for the chromosome. This step is repeated for each possible chromosome () in the population size. The elite chromosome of the generation is determined as part of the evaluation process.
Step 4 (sorting). The total fitness value is calculated for the population size based on Step 3. Calculations are carried out in ascending order from “Min” to “Max” for the total fitness values of all chromosomes.
Step 5 (tournament candidate selection). The chromosomes are selected for the next generation using the tournament selection scheme, which is a proportionate random selection method. Two parent chromosomes are selected from a population according to their fitness function values.
The main characteristics of tournament selection are summarized as follows.(i)Tournament selection is quite useful in certain situations, such as multiobjective optimization.(ii)Tournament selection uses only local information.(iii)Tournament selection is easily implemented with low time complexity.(iv)Tournament selection can be easily implemented in a parallel environment.However, tournament selection also suffers from selection bias, which means that the best one will not be selected if it is very unlucky.
Step 6 (crossover selection). Crossover probability crosses over parents to form new offspring (children). In the crossover phase, all chromosomes (except for the elite chromosome) are paired up and crossed over with a probability crossover rate. Crossover is accomplished by choosing a site randomly along the length of the chromosome and exchanging the genes of two chromosomes (parents) for each gene past this crossover site. The steps for crossover are as follows.(1)Identify the number of different requirements for manufactured products ().(2)Create the one random number between product type numbers () as follows: (3)Identify the gene code number () for the selected product obtained from the previous step.(4)Search and identify the gene code from parent A based on random product selection and transfer these gene codes to child B, which is exactly at the same gene location.(5)Search and identify the gene code from parent B based on random product selection and transfer the gene code to child A, which is exactly at the same gene location.(6)Transfer the remaining gene code from parent A to the gene blank of child A.(7)Transfer the remaining gene code from parent B to the gene blank of child B.(8)Calculate the number of crossovers based on the crossover rate () from Step 9 and population size (PS) using (9)Once new offspring have been created, the new offspring will have previous chromosomes in the current generation.
Step 7 (elitism). In crossover operation, the worst or weakest chromosomes will fade away, whereas the characteristics of the chromosomes will change continuously during mutation operation. The elite chromosome will not be subjected to mutation in the next generation. Consequently, GA does not lead to annihilation since several chromosomes (one, two, or three) from each generation are transferred directly to the next generation. Mutation is not applied on chromosomes which are immune. It is possible to maintain a fixed fitness value in some generations, but they will never deteriorate. The first three best chromosomes are selected in the elitism step.
Step 8 (mutation). Following the crossover operation, the genes will mutate to any of the codes with a mutation rate for each gene in the chromosomes, with the exception of the elite chromosome. When the crossover and mutation operations are complete, the chromosomes will be evaluated for another round of selection and reproduction. Considering elitism and after identifying the parts in which mutation will be applied, the number of mutations in each generation is calculated using (14) based on the mutation rate (), population size (PS), and maximum gene code () as follows:
After identifying the number of gene mutations, a set of rules needs to be devised for the mutation of genes from point A to point B and vice versa, while focusing on the stability of the chromosome sequence. It is important that the sequence of the chromosomes is not displaced. There are eight sets of rules for this step which are classified into two groups as follows.(1)Four rules are used to check the mutation based on the part sequence (Rules 1 and 2) and product assembly (Rules 3 and 4) of genes from A to B.(2)Four rules are used to check the mutation based on the part sequence (Rules 5 and 6) and product assembly (Rules 7 and 8) of genes from B to A.
The first group of gene mutation from A to B is illustrated in Figure 4. Likewise, the second group of gene mutation from B to A is illustrated in Figure 5.
Note.(i)Care should be taken to ensure that the position of A is before B.(ii)The above eight mutation rules should be checked thoroughly to ensure their ideal applications. In the adverse case, the above eight rules ought to be done over.
Step 9 (new population size). Having performed crossover, elitism, and mutation operations, the most ideal chromosomes of the current generation are compared and evaluated to identify its total value, after checking its termination in the following step.
Step 10 (termination). The loop of chromosome generations is terminated when certain conditions are met. The elite chromosome is returned as the best solution once the termination criteria are met. The termination criteria are listed below. (1)If the number of generations reaches its maximum, the loop of chromosome generations is terminated.(2)If there are no changes in the elite solution (i.e., no changes in fitness function value), the loop of chromosome generations is terminated using
8. Computational Results
The results of the multiobjective fuzzy mixedmodel assembly line model are presented in this section. Each metaheuristic model is verified, validated, and tested using a testbed. It seems natural that one of the crucial elements for comparing heuristics is the testbed onto which the heuristics are tested. Hence, testbeds are essential when comparing two metaheuristics. Silberholz and Golden [47] studied how to create and classify instances in a new testbed and how to ensure that other researchers have access to the problems for future metaheuristic comparisons. Silberholz and Golden [47] proposed that two types of testbeds may be used for comparing metaheuristics, namely, existing testbeds and new testbeds. Existing testbeds are used when comparing a new metaheuristic to an existing one. It is advantageous to test problem instances that have been tested in previous studies. The results will be compared on a byinstance basis, allowing relative gap calculations between the two heuristics. New testbeds are developed when an existing testbed is insufficient. For instance, there are no testbeds available for cases where metaheuristics are written specifically for a new problem and thus a new testbed needs to be developed for this purpose. A new testbed also needs to be developed if the existing testbeds are insufficient, which are often due to the fact that the testbeds are too small to effectively test a heuristic. There are two points which need to be addressed when developing new testbeds, that is, the purpose of developing the testbeds and the accessibility of new test instances [47]. The purpose of a problem suite is to emulate realworld problem instances with a variety of test cases and difficulty levels. When creating a new testbed, the focus is to provide others with accessibility to problem instances. This enables other researchers to make comparisons easily, while ensuring that the problem instances are widely used. One way to ensure this is to create a simple generating function for the problem instances. Capturing the real aspects of a problem is particularly significant when developing a new testbed. In this section, the results are presented based on the development of new testbeds due to the fact that the existing testbeds are inadequate. This is followed by a discussion on the identification of the optimum scheduling time for a fuzzy mixedmodel assembly line sequencing problem.
8.1. Input Data
A hypothetical numerical example is designed to test the fuzzy mixedmodel assembly line sequencing problem. The input data of the hypothetical numerical example is given in Table 1, consisting of 50 jobs and 20 parts in order to produce four products. There are five machine tools (one lathe, two CNC, and two robots) assigned to assemble four products.

8.2. Initialization of Parameters and Fuzzy Variables
The initialization of parameters for the mixedmodel assembly line sequencing example is shown in Table 3. It can be seen that 50 jobs (; ) are required to produce 20 parts (; ) and these parts are assembled to produce four types of product (; ). The number of tool changes is 5 (; , , , , and ) and these tools are assigned to five machines (; ). The job sequence is dependent upon the part and product assembly and is described as follows. First, the job number is assigned to produce the first product, ranging from 1 to 10 (, ; , ). It will be highlighted that there are 10 jobs in this case and they are sequenced to produce three parts according to the following order. Job numbers 1, 2, and 3 are assigned to produce Part (0) ((),), while job numbers 4 and 5 are assigned to produce Part ((), ). Job numbers 6 and 7 are assigned to produce Part ((), ), whereas job numbers 8, 9, and 10 are assigned to produce Part , which is a subpart of the product assembly ((), ). Production of the second, third, and fourth products is based on the sequence described for Product , as shown in Table 2. The fuzzy processing time of each job (, ) is defined as a triplet (). The total operating time is based on the fuzzy triangular time and is required to complete the jobs sequentially when producing each part, in which each part is assigned to a machine (, ). The total operating time is defined as (; ; ; = time). The total travelling time based on the above information is defined as (; ; ; = time). The operation and travelling time are fuzzy numbers, which are indicated by “,” “,” and “.” The parameters “,” “,” and “” represent the optimistic time, normal time, and pessimistic time, respectively.

