Flexible Flow Shop Scheduling Problem with Reliable Transporters and Intermediate Limited Buffers via considering Learning Effects and Budget Constraint

In this study, a new mathematical model is presented to solve the flexible flow shop problem where transportation is reliable and there are constraints on intermediate buffers, budgets, and human resource learning effects. Firstly, the model is validated to confirm the accuracy of its performance. )en, since it is an NP-hard one, two metaheuristic algorithms, namely, MOSA and MOEA/D, are rendered to solve midand large-scale problems. To confirm their accuracy of performance, two small-scale problems are solved using GAMS exact solution software, and the obtained results have been compared with the output of the algorithms. Since the problem in this study is multiobjective, five comparative indices are used to compare the performance of algorithms. )e results show that the answers achieved using the metaheuristic algorithms are very close to the ones achieved via the GAMS exact program.)erefore, the proposed algorithms are validated, and it is proved that they are accurately designed and useable in solving the real-world problems (which have midand large-scale) in logical calculation time. By comparing the obtained results, it can be seen that theMOEA/D algorithm performs better in terms of computational time (CPU time) andMean ideal distance (MID).)eMOSA algorithm also performs better according to the index Spread of nondominated solutions (SNS), diversity metric (DM), and number of Pareto solutions (NPS). Considering the confirmation of precision and accuracy of performance of the proposed algorithms, it can be concluded that MOSA and MOEA/D are useful in solving the midand largescale modes of the problem in the study, which is very applicable in the real world.


Introduction
Scheduling is the gathering of principles, models, and method for making decisions and determining a time schedule. To do so, it is necessary to define the problem, its dimensions, and objectives accurately. Resources and activities are defined in a system in various ways. Resources can be machinery in a production workshop, air lines in an airport, worker in a construction project, and so on. If the process is done only in one stage and there is only one processor in it, the problem is a single machine mode one. Sometimes, more than one processor is used in a parallel way for optimizing the efficiency and speed of processing. In such a situation, the sequence of operation problem transforms into a problem with parallel machines [1,2]. e flexible flow shop problem includes a layout of numerous machines in a linear way so that production begins in one in the initial work stage and finishes in the last one in the last stage. All productions must be put inside one machine in all work stages and must be passed through production and process stages one by one until the last stage [3]. e sequence of operation problem can be categorized based on the way of processing and the placement and characteristics of processors. If a process is done only in one stage, and there is only one processor in it, the problem is a single machine mode one. Sometimes, more than one processor is used in a parallel way for optimizing the efficiency and speed of processing. In such a situation, the sequence of operation problem transforms into a problem with parallel machines. A more complicated mode occurs when data must be processed in numerous consecutive stages. In this mode, machines are placed in a series way, and each activity must be performed on all machines, respectively. e sequence is the same for all activities. Such a problem is called a flexible flow shop one [4,5]. Sometimes, in flow shop problems, check points occur, decreasing their efficiency. It is because processing takes a lot of time in some stages. As a result, numerous processors are used in a parallel way instead of just one to fasten and facilitate processing. Such problems with such characteristics (numerous processing stages in a series way, many processors in a parallel way in each stage, and the identical processing for all activities) are called flexible flow shop ones. e flexible flow shop problem includes a layout of numerous machines in a linear way so that production begins in one in the initial work stage and finishes in the last one in the last stage. All productions must be put inside one machine in all work stages and passed through production and process stages one by one until the last stage.
e objective of such problems is the quantification of the time of processing all jobs on machinery, so that all jobs are processed in the shortest time. As mentioned earlier, in this study, the main problem is studied under the assumption that parameters and data were predetermined. Figure 1 shows a comprehensive categorization of the sequence of operations and the current flexible flow shop problem position.
In this study, a scenario-based multiobjective problem is focused on in a production system based on the flexible flow shop. e probabilities of each scenario are determined. In each of them, the costs of each machine work and the rate of transportation failures are different in every stage. In the model, the flexible flow shop is presupposed. In all stages of the flow shop, a lot of transportation is done for transporting products and processing on machines. Each time the work is moved by the transporter between each station, it has its own transportation time. e preparation times of each job at each station depend on the previous work done at that station. Performing preventive maintenance and repair operations is another operation that can be carried out on work transport devices between workstations. In this situation, each transportation device becomes defective based on its relative failure distribution function and need for fixing. At this time, operations or productions remain on transportation devices until transformation operations finish and the device becomes available. Afterwards, transportation continues. Repairing the transportation devices is done randomly, and each problematic device will be repaired based on a probability distribution function. Regarding the constraints in the number of transportation devices among stages, intermediate buffers with limited capacities are specified for each work stage. Also, each human resource has a learning rate regarding each machine. In addition, the time of processing and preparing decreases based on human resource learning rates. e objective is the quantification of the mean of completion time value at the end of stages and the penalty of waiting time in buffer and also the cost of processing based on jobs' specification. e number of works performed in intermediate buffers is fewer than the maximum capacity of the relative accumulator. Another constraint affecting the proposed flexible flow shop problem is the budget. e budget limit in this issue is applied in such a way that the total budget available for the maintenance costs and staying of the works in the in intermediate buffers before the workstations is a limited and specific budget. Each work waiting for processing on work shop stages' previous machines in intermediate buffers has a certain cost per time units in the problem, and for all works in the intermediate buffers of each work stage, a maintenance cost is determined for the model. e total cost must be less than the total budget of the project.
In the second section, the review of literature is presented. In the third section, the mathematical model is introduced, and the results will be presented in the fourth section. At the end, conclusions and suggestions are mentioned.

Literature Review
e first research regarding the flexible flow shop problem was presented in 1970s by [6], in which the modeling of the production system in a textile supplying factory is studied. After the introduction of this branch of sequence of operation, researchers became gradually interested in it and published numerous articles [7].
Johnson [8] proved that the problem of flexible flow shop in which the function of the quantification of the make span is under the set of NP-hard ones. As a result, researchers were encouraged to render various heuristic and metaheuristic methods to approximately solve problems. Tavakkoli-Moghaddam et al. [9] proposed a two-way model of the work flow problem of parallel machines considering independent initial set up times and the precedence constraint and used a genetic algorithm to solve it.
Naderi et al. [10] presented two algorithms based on electromagnetic one and simulated annealing one to solve the two-objective problem of job shop flow. e objectives were quantifying the completion total time of jobs and the whole delays.
e results show that the electromagnetic algorithm had better answers compared with other methods. Wang and Choi [11] presented a multiobjective model of parallel machine problem considering the maintenance operation and preventive multi fixings. e objective function included the quantification of the completion time of all jobs; the first case is the unreliability of machines and the second case is the unreliability of preventive maintenance and repairs. To solve the problem, they used a nondominance sorting genetic algorithm (NSGA). Shahidizadeh et al. [12] presented a two-objective model of a job shop flow with parallel machine, in which the process was a group production one. e objective function included the quantification of the completion time of all jobs, the reduction of delay, and the cost of purchasing machinery. Since their problem was NP-hard and had large-scale, they used a multiobjective harmony search algorithm to solve it. e results showed that their proposed method rendered more qualified answers compared with other algorithms. Liao et al. [13] presented a multiobjective model of the group scheduling problem in the job shop flow system considering the maintenance operation and preventive repairing. e objective function included the quantification of the completion time of all jobs, maintenance costs, and preventive repairing. Zandieh et al. [14] presented a multi objective model of hybrid job shop flow considering the maintenance operation and preventive repairing. eir solution included a hybrid dominance sorting genetic algorithm with two approaches. e results showed that their proposed method was very efficient in hybrid job shop flow problems. Mollaei et al. [15] presented a two-objective model for multistage problems in production systems. To solve it, they used a hybrid of Taguchi and Monte Carlo methods. eir study was about an automotive company in France. e results proved the quality of their solution in the mentioned study. Ozsoydan and Sagir [16] proposed an integrated Greedy algorithm enhanced by metaheuristic-based learning for flexible flow shop problem, considering independent initial set up times. Cheng et al. [17] presented a group scheduling problem in a job shop flow one. ey presented two algorithms of multistage simulated annealing and local searchbased variance ones. Also, by solving the problem in the previous literature and studies using their method, they indicated that most local search-based variant algorithms render better results, both considering the objective function and computational time. Shen et al. [18] presented a twoobjective model for the flexible flow shop problem in the macaroni production industry in Belgium.
e objective function of his study included quantifying the time of job completion and the energy and human resources costs. ey proposed a customized two-phase genetic algorithm. e results of solving their proposed problem indicated that their algorithm renders better results compared with the Pareto power-based evolutionary one. Pagnozzi and Stuzle [19] rendered a hybrid probability local optimization algorithm to solve their proposed job shop flow problem. e objective function included quantifying the time of job completion and the total of delays.
Lotfi et al. [20] have researched renewable energy. e most important innovations in their study included the use of robust two-level programming techniques and game theory (Stackelberg Competition) for locating renewable energy sites. e results show that the combination of uncertainties can increase energy production and supplier profits. In addition, the objective functions of the proposed model are compared with those under uncertain conditions. e sensitivity analysis of the main parameters is performed to validate the proposed model. As uncertainty increases, the energy produced decreases, and the supplier's profit increases. Supplier profits gradually decrease as the discount rate increases. In addition, as the scale of the problems increases, the energy produced and the profit of the supplier increase.
According to the above-mentioned studies, a comprehensive categorization of studies on the flexible flow shop is summarized in Figure 2. e problem in this study is a scenario-based flexible flow shop one in which reliable transportation devices and intermediate limited buffers are considered in addition to budget constraints, dependent setup time, the learning effect, and also the objective functions of "the mean of completed time weight value," "the mean of the penalty of waiting time in buffer," and the total specified processing costs. Reviewing the literature showed that the problem in this study has not yet been focused on in any research and is completely new.

Methodology
In this study, a multiobjective scenario-based problem that exists in a production system is presented. e system is based on a flexible flow shop. In this model, the flexible flow shop problem is presupposed, and inside each job stage, there are some devices for transportation, each having its own specific times. Preventive repairing and maintenance is one of the other operations performed on transportations used in transporting productions. In this situation, each transportation device breaks down based on its relative distribution failure function and need for fixing. At this time, operations remain on transportation devices until transformation operations finish and the instrument becomes available. Afterwards, transportation continues. Repairing the transportation devices is done randomly, and each problematic device will be repaired based on a probability  In the following section, the mathematical model and all its parameters and variables will be explained in detail. Considering Figure 3 of the model, suppose that the job j is being processed in the stage S. In each stage, N s number stages are placed in a parallel way, and all jobs must be processed on one machine in each of job stages. Preventive maintenance and repairing operation can be performed on all machineries in job stages. Also, there is an intermediate limited bu er where jobs have zero start time or end time. According to Figure 3, which is a proposed Gantt Chart sample of the problem with three job stages and 5 productive machines, jobs 1, 2, and 3 are in the work stages 1, 2, and 3, respectively, with machines. As observed, on machine 1 in work stage 1, the preventive maintenance operation is done after the processing of job 3 was completed. Another preventive maintenance operation was done on job stages 2 and 3. Job 2 was done in the job stage 1 between times 2 and 3. It started from time 6 in the job stage 2. erefore, this job remains in the intermediate bu er of job stage 1 until the process in job stage 2 begins. If there is no capacity to accept job 2 in the intermediate bu er, this job will remain in the machine, and as a result, it will block. In some cases, it is possible that more than one job is stored in the intermediate bu ers. For instance, in the intermediate bu er of job stage 2, jobs 4 and 5 are stored between times 13 and 14.
e mathematical model used in this study is derived from the one in Zabihzadeh and Rezaeian [21]. is model is developed here. is model includes two objective functions and 17 constraints. e variables in this study are of positive and binary types. ere are 6 positives and 2 negatives. All indices, parameters, and variables (whether dependent, independent, or control) are de ned in the basic model.
Assumptions and constraints in this study are as follows: (i) All jobs j that are indicated by J {j|j 1, 2, . . ., J} are available at the start of jobs (time-zero) and no job stops before nal production. (iii) e machinery group shown by j is placed in a parallel way in each job stage to be processed and used for producing materials. (iv) All productions must be put on all job stages and speci ed on each of the stages for processing and producing to be done on them. In other words, when processing begins on a job in one job stage, it must continue onstage by stage until the last one. (v) Each job has a preparation time depending on the previous one.

Binary Variables
X jomswt : if the j th job in the o sequence is specified to the m th machine and the w human resources in the s stage in the t th time period, it equals 1. Otherwise, it is zero. U jrtk : if the r th transportation device is specified to the j th job in the t th time period under the k scenario, it is 1. Otherwise, it is zero.

e Objective Function and the Constraints of the Problem.
In this section, the suggested mathematical model is presented as follows: e first objective function calculates the mean of the time weight value at the end of a stage as follows: e second objective function shows the penalty of waiting time on buffer and the cost of the j th job (production) process based on its specification.
Subject to o∈O m∈M s w∈W s j∈J m∈M s w∈W s o∈O m∈M s o∈O j∈J m∈M s∈S 8 Complexity e constraint (3) specifies each job (production) to one sequence and one machine of each stage. Constraint (4) guarantees that each sequence must be specified to one job in each stage. Constraint (5) guarantees that each human resource can only be specified to one machine. e constraint (6) guarantees that the next sequence can only be specified only when the previous one is done. Constraint (7) calculates the availability of the transportation device of each time period under any scenario. Constraint (8) calculates the mean of transportation time based on the availability of transportation devices. Constraint (9) determines the type of the transportation device in order to transport the j th job in the t th time period. Constraint (10) determines the maximum of investment for each time period. Constraint (11) calculates the time of completion of each job in the first stage and for the first sequence. Constraint (12) calculates the time of completion of each job in the stages, except the first one, and for the first sequence. Constraint (13) calculates the time of starting of the j th job as soon as it is completed in the previous stage. Constraint (14) calculates the time of j′ th job start based on the completion of the previous one in the previous stage and the previous stage in the same machine. Constraint (15) guarantees that if X jomswt becomes zero, the time of completion will become zero too. Constraint (16) calculates the time of completion of the j th job in the previous stage. Constraints (17) and (18)

The Numerical Results
Since the mathematical model in this study is an NP-hard one, to solve problems of the real world, two metaheuristic algorithms, namely MOSA and MOEA/D, are used to in this section, and the results will be analyzed and compared. To evaluate the quality and data sparsity of the multiobjective metaheuristic algorithms, there are various metrics, which are different from the ones of the one objective metaheuristic ones. e reason is that there are answers with no priority over each other. e metrics used in this study are as follows: the number of Pareto solutions (NPS), diversification metric (DM), mean ideal distance (MID), spread of nondominance solutions (SNS), and computational time (time). To make these comparisons, 10 experimental problems are made, which are small, medium, and largescale, and the results of their solutions via algorithms are mentioned in tables. Based on Table 1, the parameters of the problem are randomly made and used as primary data for problems. e primary data of the problem are as follows: process time, transportation time, probability of each scenario, preparation time, the job costs, the amount of breakdown and repairing of the transportation devices, the learning coefficient of human resources, the penalty (cost) of jobs, the coefficient of the importance of jobs, the maximum investment in each time period, and the time for preventive  1, 15)) CAPDE t � U (10,20) α wm � U(0.1, 0.3) p k � U(0.1, 0.8) Tr rs � round (U(1, |T| − 1)) RR jj′t � round(U(500, 700)) ST jj′t � round(U(1, 10)) F � round (U(8000, 9500)) Complexity 9 repairing.
e distribution function of each of them are mentioned in Table 1.

Adjusting the Parameters for the Solution Algorithms.
Adjusting parameters is vital for the e ciency of metaheuristic algorithms. erefore, determining apposite parameters for metaheuristic algorithms is an important step in making them. erefore, a set of calibration experiments are usually performed for nding the optimized mode of various amounts of control parameters of algorithms. Since the increase in test levels and the number of factors leads to an exponential increase in time and cost, the Taguchi method is used to perform tests. For the two metaheuristic algorithms used in this study, the control parameters are as follows: in the case of the MOSA algorithm, the parameters are the maximum of iteration (MaxIt1), maximum of internal iteration (MaxIt2), primary temperature (T0), the coe cient of heating (Alpha), the maximum capacity of answer archive (nRep), number of grids in each dimension (nGrid), the leading choice coe cient (Beta), and the choice coe cient for deleting archive (gamma). In the case of MOEA/D, parameters include the maximum of iteration (MaxIt), number of population (nPop), maximum number of archive (nArchive), and the coe cient of creating harmony in the intersection operator (gamma). Each parameter considered for the proposed algorithms has a salient e ect on the quality of answers and time. For example, if a large number of population is considered for MOEA/D, the algorithm will have a long range of answers, and consequently, it can take a long time to achieve them. If a small number of populations is considered, the optimized answer could be locally optimized. Regarding the experiments and results, the e ective range of each control parameter of algorithms has been determined.

e Design of Experiments Using the Taguchi Method.
After determining the e ective range of the control parameters, a set of design experiments (DOE) is designed using the Taguchi approach to nd the e ects of these parameters on the performance of the proposed algorithm and achieving their optimized mode.
In this study, the number of Pareto solutions (NPS), diversi cation metric (DM), mean ideal distance (MID), spread of nondominance solutions (SNS), and time (Time) are used to compare the performance of the suggested algorithms.
e necessary experiments for analyzing the various modes of parameters and relative answers are presented in Table 2 of each algorithm. e data were analyzed using MINITAB20, and the results are shown in Figures 5 and 6 and Table 3.
Subsequently, the same procedure is repeated for MOEA/D. Table 4 indicates the test level of each parameter, and Figures 7 and 8 and Table 5 shows the results.

Evaluating the Performance of Multiobjective Algorithms.
To measure the quality and diversity of multiobjective metaheuristic algorithms, there are various comparison metrics. In multiobjective optimization, evaluating the algorithm performances could be di cult since there are di erences presented in the nal answers and paradoxes in objectives. When the visual analysis of results is very difcult, it is vital to have various metrics to evaluate the performances to recognize the best set of nondominance solutions.
To compare the capability of various algorithms in solving multiobjective problems, various methods and instruments are used. Some of them, which were used more in the multiobjective literature, will be discussed in the following session; Number of Pareto solutions (NPS): this metric calculates the number of nondominance solutions achieved in every run of algorithm. According to it, the more the number of nondominance solutions, the better the algorithm performance. Diversi cation metric (DM): this metric is about the range of Pareto solutions of algorithms. It can be calculated by the equation (23). e more the DM, the better the algorithm's performance. Signal-to-noise: Smaller is better   Mean ideal distance (MID): the closeness of the achieved Pareto set is calculated by the algorithm to the optimized Pareto edge. Since achieving the Pareto edge is impossible for most problems, the distance between Pareto points and the ideal ones, (0,0), is calculated via this metric. e equation is as follows:

Complexity
In the above-mentioned equation, n is the number of achieved nondominance solutions, and i � 1,2, . . ., n Regarding MID, the more it is, the farther the algorithm is from the ideal Pareto edge. Spread of nondominance solution (SNS): the spread among the solutions in the nondominance solutions achieved by the algorithm is calculated by this. e equation is follows: e more SNS in the algorithm, the more the algorithm spreads, which is more favorable. Computational time (time): the time for the performance of an algorithm. e lesser time the algorithm takes, the better.
In the following sessions, examples are categorized as small-, medium-, and large-scale ones regarding the model nature. For this purpose, two small scale examples, four mid-dimensional, and 4 high-dimensional ones were randomly made to evaluate the model when its design is real and parameters are random. ese problems are listed in Table 6. e examples one and two are solved using metaheuristic algorithms, and their results will be compared with the ones achieved using the GAMS exact program to validate the answers. e purpose is to confirm whether to continue with metaheuristic algorithms. Also, the Pareto front chart of each problem is presented, and the five metrics (MID, DM, SNS, NPS, and CPU time) are considered in it.
e results are shown in Tables 7 and 8 and Figures 9-14. F 1 is the weight objective function of the completion time, and F 2 is the objective one of the total specification cost and waiting time for jobs.

e Comparison of Algorithm Results for Small-Scale
Problems. Each algorithm was run for five times, and the best results were considered. e results of comparison between the proposed algorithms and one of the GAMS exact programs used for low-dimensional problems are mentioned in Table 9.
As mentioned earlier, the more DM, SNS, and NPs and the fewer MID and time, the more efficient the algorithm with respect to the metric. According to the results mentioned in Table 9 for problems 1 and 2, for time and SNS, the results of GAMS are better. In case of MID and DM, GAMS and MOEA/D have proved more fruitful. Regarding the NPS metric, the results achieved from GAMS and MOSA are relatively better. Nevertheless, since there is not much difference between results achieved from the GAMS and the ones of the proposed algorithms regarding any metrics, it can be concluded that they are well-designed and efficient, and thus, they are useable in solving problems with high dimensions.  1  2  3  3  2  2  3  2  3  2  2  4  3  3  3  4  3  3  3  4  5  3  3  3  5  4  4

e Comparison between Algorithms for Mid-and High-Dimensional Problems.
e results of the proposed algorithms used on mid-and high-dimensional problems after ve times of running are shown in Table 10 and Figures 15-19. ese results are based on the best answers achieved from algorithms.        According to Table 10 and Figures 15-19, MOEAD has a better performance regarding MID and time. On the other and, MOSA performs more fruitfully regarding SNS, DM, and NPS.

Conclusion
According to the studies on scenario-based flexible flow shop systems, a multiobjective problem has never been considered for the following characteristics: reliable transportation devices and constraint in capacity of intermediate buffers, intermediate limited buffers, dependent set-up time, and learning effect and budget constraints. Considering the reliability constraints of transportation devices and one of the finances has resulted the problem in becoming closer to the ones in the real world. e budget issues, like production costs, which are important for most companies and production factories, are also considered in this study. ey can be considered in many of such problems in the real world. In the proposed model, scenarios are supposed with definite probabilities in each of which the cost of jobs and the amount of failure of transportation differ for each machine in different stages. In this study, it is supposed that in job stages, there are many transportation devices on machines for transporting and processing jobs and productions. e transportation of jobs and production in different stages has a specified time. e time of transportation of productions, which can be done by conveyors, lift-tracks, or automatic guided vehicles (AGV) is different. e reliability of transportation is also considered in this problem. It means that not all transportation devices are available all the time and preventive maintenance operations can be done on them. Also, because of the limited number of means of transportation between workstations, an intermediate limited buffer has been considered for each workstation. e number of jobs in each buffer must be less than the maximum capacity of them. Each time the jobs cannot be processed on machinery since machines being busy or under preventive maintenance operations, they remain in buffers. For each time unit that jobs wait to be processed in intermediate buffers, an amount of maintenance cost is considered for the project. e total budget for maintaining jobs in intermediate buffers is limited, and their total cost must be less than the total available budget. Each job is done depending on the previous one in each stage, and there is a preparation time for them. ere is a learning rate for human resources regarding each machine, based on which processing and preparation time decreases. According to the above-mentioned issues, it can be concluded that parameters affecting the objective functions in the problem are as follows: processing time on productive machines, time of transportation using transportation devices, dependent preparation times, learning effect, waiting time of jobs in intermediate buffers, and reliability of transportation devices. As a result, the importance of discovering the best sequence of jobs and optimizing the objective functions becomes clear.
In this study, a new mathematical model with reliable transportation and constraints in intermediate buffers regarding budget and the effects of human resources learning was presented for solving the flexible flow shop problem. Firstly, the model was validated to confirm its accuracy of performance. en, since it is an NP-hard one and the exact methods are not useful in logical calculation time, two metaheuristic algorithms, namely MOSA and MOEA/D, were presented to solve mid-and high-dimensional problems. To ensure the proposed algorithms accuracy of performance, two small-scale problems were solved by GAMS exact solution software, and the results obtained were compared with the output of the algorithms. Since the problem in this study is multiobjective, five comparative indices, nsmely the number of Pareto solutions (NPS), diversification metric (DM), mean ideal distance (MID), spread of nondominance solutions (SNS), and computational time (CPU Time), were used to compare the performance of algorithms. e results show that the answers achieved using the metaheuristic algorithms are very close to the ones achieved via the GAMS exact program. erefore, the proposed algorithms were validated, and it was proven that they were accurately designed and useable in solving the real-world problems (which have mid-and high-dimensions) in logical computational time. Comparing the answers with the ones in the real-world, it was observed that regarding the time and MID indices, the MOEA/D algorithm, SNS, DM, and MPS indices, MOSA had relatively better performance. Considering the confirmation of the accuracy of the proposed algorithms previously done by comparing the answers with the ones of GAMS, it can be concluded that both MOSA and MOEA/D are useful in solving the mid-and high-dimensional modes of the problem in the study, which is very applicable in the real world.

Recommendations.
e flexible flow shop problem can be developed considering other assumptions and parameters in machinery, in which, among them, the following are mentioned. e mathematical model of the present problem can be developed by considering assumptions, such as independent time for initial set up, failure, maintenance operation, process of operations considering the constraints in availability, preventive repairs on productive machinery regarding maintenance and repairing costs, the time needed for machine installation, process of operations with availability constraints to machinery, availability time the legal change in sequence, intervals for restarting, and so on. Also, another objective function can be considered for the quantification of delays or time consumption in intermediate buffers or total time of operations. In addition, in solving the problem, other multiobjective metaheuristic algorithms can be used, such as neural network, ant colony, tabu search, water flow-like, etc. Metaheuristic and heuristic algorithms, such as branch and bound, Lagrange or a hybrid of two or more algorithms, and a hybrid coding, can also be used to solve such problems.

Data Availability
Data will be available upon request to the corresponding author. 18 Complexity