An Effective 4–Phased Framework for Scheduling Job-Shop Manufacturing Systems Using Weighted NSGA-II

: Improving the performance of manufacturing systems is a vital issue in today’s rival market. For this purpose, during the last decade, scientists have considered more than one objective function while scheduling a production line. This paper develops a 4-phased fuzzy framework to identify effective factors, determine their weights on multi-objective functions, and, accordingly, schedule manufacturing systems in a fuzzy environment. The aim is to optimize product completion time and operational and product defect costs in a job-shop-based multi-objective fuzzy scheduling problem. In the ﬁrst and second phases of the proposed framework, it was shown that the existing uncertainty of the internal factors for the studied cases causes the weights of factors to change up to 44.5%. Then, a fuzzy-weighted NSGA-II is proposed (FW-NSGA-II) to address the developed Non-linear Fuzzy Multi-objective Dual resource-constrained scheduling problem. Comparing the outcomes of the proposed method with other solving algorithms, such as the Sine Cosine Algorithm, Simulated Annealing, Tabu Search, and TLBO heuristic, using seven series of comprehensive computational experiments, indicates the superiority of the proposed framework in scheduling manufacturing systems. The outcomes indicated that using the proposed method for the studied cases saved up to 5% in the objective function for small-scale, 11.2% for medium-scale, and 3.8% for large-scale manufacturing systems. The outcomes of this study can help production planning managers to provide more realistic schedules by considering fuzzy factors in their manufacturing systems. Further investigating the proposed method for dynamic product conditions is another direction for future research.


Introduction
The job-shop facility location can be considered among the industry's most frequently used facility location mode, applied in many industry sectors.
With the Industrial 4.0 revolution, the need to use heuristics and meta-heuristics to schedule manufacturing systems became increasingly obvious [1]. As a result, in recent years, many references used meta-heuristics as a core scheduling mechanism or part of it [2][3][4].
In a job-shop scheduling problem (JSSP), each in-process product passes through different shops, according to its operational requirements, where, in each shop, several parallel machines can serve an activity [5]. Traditionally, the main objective of job-shop scheduling problems was to find the best part routes of raw material throughout different shops until a product is completed [6]. However, today, various objective functions are taken into account in more complicated manufacturing environments [7]. Product completion time, quality, job tardiness, human resource scheduling, and operational cost can be considered the most frequently used objective functions during the last two decades.

•
The first is to identify the influential internal factor influencing the scheduling of dual-resource-constrained job-shop-based manufacturing systems under uncertainty.

•
To determine the effects of each uncertain factor on scheduling the multi-objective function of this research (time, defect, and HR cost).

•
To propose a new method to schedule the dual-resource constrained manufacturing system using the found fuzzy weights.

1.
This research deals with fuzzy scheduling of job-shop scheduling which has previously been less developed. 2.
In this study, a multi-objective non-linear mixed integer programming method is developed to minimize the makespan and completion cost and maximize the quality of products simultaneously.

3.
A hybrid NSGA-II-SA, weighting, and Lp-metric method are proposed to solve the model. Note that the Lp-metric and hybrid NSGA-II-SA have not previously been used for job-shop scheduling problems.

Managerial Implications
The outcomes of this research will help managers of small and medium-scale manufacturing systems to schedule their job-shop-based systems more effectively by considering the fuzzy internal factors that can influence their systems. Moreover, the proposed method in phases 3 and 4 of this framework will help the managers find near-optimum solutions that can simultaneously improve completion time, product defect, and human resource and operational costs.

Importance and Necessity
In a real industrial environment, internal factors can influence the efficiency of the scheduling process, resulting in unrealistic schedules. However, such factors are usually ignored while scheduling a dual-constrained manufacturing system. In addition, the internal factors are uncertain, and such uncertainties for small elements may consequently cause increased uncertainty in a production schedule of a DRCSP [9]. As shown by Section 2 Mathematics 2022, 10, 4607 3 of 28 (literature review), there is no evidence of internal factors' role in scheduling efficiency. Therefore, it is crucial to propose a new multi-step framework to identify influential factors, classify them, measure the levels of uncertainties, and schedule a DRCSP accordingly.
The rest of the paper is organized as follows: (i) In the first section, an in-depth review of the literature on scheduling manufacturing problems will be carried out. (ii) In the next section, a new 4-phased framework will be proposed; where, in the first phase, the influential factors in scheduling manufacturing systems will be identified, followed by determining the weights of each factor using data analytics methods (phase 2). In continuation, a new multi-objective mathematical model will be developed (phase 3). In the last phase, a new fuzzy-weighted NSGA-II will be proposed for scheduling manufacturing systems in the fuzzy environment (Phase 4). (iii) The proposed framework will then be applied to several case studies, (iv) followed by statistical measuring to determine the proposed framework's performance.

Literature Review
During the last two decades, many scientists have focused on job-shop scheduling issues, where different objective functions were taken into account under various conditions. In the following, several opted references will be reviewed to determine the influential factors for the first phase of the proposed framework in the next section. Completion time is the main objective of job-shop scheduling, where the aim is to find the fastest production cycle by determining the best series of machines in different shops to complete parts.
Recently, a new concept in job-shop facility design emerged which focused on performing operations served by more than one machine. A vital issue for manufacturing systems is to respond to customers' needs as quickly as possible. Therefore, scheduling processes are required to create more real-time schedules. In this regard, what makes job-shop scheduling important is the way of completing the required operations according to the operational process chart (OPC), where, for each operation, an appropriate machine should be selected among a set of parallel machines in a shop and then materials are transferred to visit a set of determined machines, accordingly. Delgoshaei, Ariffin, et al. focused on material movement problems in manufacturing systems by giving special attention to voids and idle machines [10]. Additionally, important factors play a vital role in scheduling job-shop-based manufacturing systems. In the following, the important factors will be reviewed; such factors will be used in the next section as the important factors for the proposed 4-phase framework.

Makespan
Product completion time is the first objective in classic job-shop scheduling problems. Liu et al. proposed a multi-objective-based Genetic algorithm framework (MPMOGA) for addressing multi-objective job-shop scheduling where product completion and total tardiness were considered [11]. In continuation, Kianpour et al. focused on due dates in jobshop scheduling while uncertain processing times were considered [12]. Later, Ying et al. proposed meta-heuristics for minimizing makespan in rapid job-shop-based manufacturing systems [13]. Then, Wang et al. proposed fuzzy job-shop scheduling where multi-objectives were considered [14]. Afsar et al. investigated the multi-objective scheduling of green job-shop scheduling by proposing hybrid Memetic and NSGA-II algorithms [15]. With multi-processing speed, Zheng et al. applied an NSGA-II to a bi-objective DRCJSP issue to reduce makespan and overall electricity consumption [16]. Their research discussed that a stochastic machine's working process leads to non-deterministic job production time. This means that each activity's processing time is determined by the machine deployed and the worker allocated to it.

Setup Times
Loading and unloading parts to a device other time-consuming activities that are usually ignored during the scheduling. This challenge was addressed by Akbar et al. using Mathematics 2022, 10, 4607 4 of 28 a GA in dual-resource constrained scheduling problems where the aim was minimizing the completion time. [17]. In addition, workers vary in their loading and unloading times [18]. Such differences may happen as a result of skill, age, and experience. In many manufacturing systems, the setup of machines is based on product sequences (sequencedependent) which influence the total completion time of a product accordingly [19].

Machine Amortization
Machine amortizations are the other important issue that may cause serious problems for effective scheduling [20]. Gupta et al. proposed an integrated method for scheduling job-shop-based manufacturing systems while machine failures played a vital role in scheduling [21]. Amelian et al. addressed the failure-prone job-shop scheduling using an NSGA-II algorithm [22]. Souza et al. focused on preventive and corrective maintenance in scheduling job-shop manufacturing systems where uncertain constraints existed [23].

Raw Material Quality
In their novel review, Mohan et al. considered raw materials the other important factor that can cause a job-shop manufacturing system's failure [20]. For this purpose, many scientists considered the study of raw material requirements and specifications as the objective of their scientific activities. Zhang et al. focused on energy efficiency in the optimal control problem by developing a multi-objective model [24]. Ramya et al. integrated material requirements planning in job-shop scheduling using the Shuffled Frog Leaping Algorithm. For this purpose, they developed a new application with VB.NET and MYSQL [25]. Abderrahim et al. proposed a heuristic for scheduling job-shop systems where processing resources and transportation were considered the main objectives [26]. Vaez developed a new method of scheduling industrial systems that considers job etiquette and fines for lateness [27]. H. Li et al. argued that most scheduling references focused on the time factor [28]. However, essential objectives should not be ignored while studying job-shop-based systems. For this purpose, they proposed a meta-heuristic considering energy-consciousness during job-shop system scheduling.

Material Transferring
Material transferring (handling) is another factor influencing completion time and system costs in job-shop systems. Saidi-Mehrabad et al. developed a mathematical model for scheduling process products in job-shop systems by considering movement times between and within shops [29]. Karthikeyan et al. proposed a meta-heuristic to solve the FJSP problem where machines can operate along many pathways [30]. Delgoshaei, Ali, et al. argued that improper scheduling in job-shop systems might cause machine imbalance in a shop, consequently decreasing system productivity while increasing machine emergency service costs. Moreover, they demonstrated that dynamic cost conditions influence material routing in the process and may cause machine load variance [31]. While re-entrant portions are allowed, Rabbani et al. employed NSGA-II for scheduling difficulties [32]. Later, Huang focused on the role of automated guided vehicles in scheduling job-shop systems [33]. Burdett et al. developed an operating coal export terminal as a flexible job-shop system where the aim was a material recall from terminals [34]. Sun et al. focused on deadlock and robot movements in scheduling automated job-shop-based manufacturing systems [35].

Workforce
In this research, human resources are considered the main workforce category, consisting of two parts: human resource skills and human resource costs. Although Industry 4.0 is built on quick allocation and sequencing decisions, it is asserted that the workers should still be considered a crucial component for coordination and management, necessitating exact schedules for each worker [36]. Worker flexibility is the main concept in the human resource scheduling of job-shop-based manufacturing systems [37]. Delgoshaei et al. addressed a meta-heuristic for solving worker assignment problems in manufacturing systems while worker skills were considered [38]. Dhiflaoui et al. addressed worker flexibility in dual-constraint job-shop scheduling problems [39]. In the same year, Gong et al. proposed a GA and NSGA-II to solve the double flexible job-shop scheduling problem where workers and machines were considered [40]. Shahgholi Zadeh et al. used the Bee Colony optimization method to address double flexible job-shop scheduling problems while various activity processing times were possible [41]. On the other hand, Kress et al. argued that workers' skills should not be ignored during a system's job-shop scheduling [42].

Literature Review Gap and Problem Statement of the Research
An in-depth review of the opted references indicated that the issue of scheduling jobshop problems while optimizing completion time, product defect, and human resource and operational costs using a fuzzy environment to address the uncertainty of factors has not been developed before. Therefore, this paper will provide a base to generate an alternative schedule for optimizing the multi-objectives above using a fuzzy inference system. This research continues with Saidi-Mehrabad et al. who determined that scheduling is vital when product completion time, product defect, and human resource and operational costs in a fuzzy environment are considered in a fuzzy scheduling problem [29]. This research aligns with the core of Industry 5.0, where human resources are assigned to work better and faster using an advanced scheduling method [43].

A 4-Phased Framework
In this stage, a 4-phased framework will be developed. Figure 1 shows the flowchart of the proposed framework. That the main concept of the research is inspired by Saidi-Mehrabad et al., w dressed a new method for optimizing the transportation times between machines shop scheduling problem [29].  This method starts with identifying the effective factors influencing the multi-objective fuzzy scheduling of job-shop-based manufacturing systems. For this purpose, a series of opted references were reviewed. Moreover, an interview with academic gurus and industrial owners will be performed. In continuation, statistical data analysis will be carried out to determine the weights of the factors. Then, in the next step (phase 2), a fuzzy inference system will be set up using Matlab ® 2016 to determine the fuzzy weights of each influential factor. A new mathematical modeling program will be developed using fuzzy weights where multi-objective decisions are considered. The complexity of the developed model will be measured to propose the appropriate solving method. Then, in phase 4, a weighted NSGA-II will be presented to solve the model where small, medium, and large-scale case studies are considered. The performance of the proposed method will then be measured using some indexes.

Assumption
That the main concept of the research is inspired by Saidi-Mehrabad et al., who addressed a new method for optimizing the transportation times between machines in a job-shop scheduling problem [29].

Assumption
Several assumptions are taken into account as follows to determine the optimum multiobjective function for the job-shop-based manufacturing systems under fuzzy conditions:

•
The manufacturing system's facility design is a job-shop, which means parallel machines exist in each shop that can perform a similar service.

•
The number of machines in a shop is fixed and cannot be changed during manufacturing. • Each product must receive services from parallel machines in the related shops.

•
The internal factors are uncertain and might be varied due to various production circumstances.

•
Each part can receive one service at a time.

•
Each service can be completed by only one machine at a time.

•
Part re-entrant is not considered in this model but can be regarded as an opportunity to expand the model for future studies. • Each machine has a defect rate considered during services (service rate).

•
Raw materials that have a defect rate may cause product failure.

•
The period's beginning inventory should be zero, and the period's last inventory should also be zero.

Identify the Effective Factors (Phase 1)
This section will identify the influential internal factors that can influence the multiobjective function. For this purpose, several factors that are mostly taken into consideration by scientists are considered. Moreover, interviews with academic scientists and industry owners in Malaysia were conducted. Figure 2 shows a list of the internal factors that are identified by the literature review and interviews with the scientists and industry experts: This section will identify the influential internal factors that can influence the multiobjective function. For this purpose, several factors that are mostly taken into consideration by scientists are considered. Moreover, interviews with academic scientists and industry owners in Malaysia were conducted. Figure 2 shows a list of the internal factors that are identified by the literature review and interviews with the scientists and industry experts:

Questionnaire Designing and Data Gathering
After identifying the factors, a statistical survey is performed to determine the crisp weight of each factor. For this purpose, a questionnaire is designed. The Cronbach's alpha values showed the reliability of the questions before gathering data from the real-world, which consisted of various automotive, house appliance, and metal forming industries that use job-shop facility locations. Table 1 shows Cronbach's Alpha:

Questionnaire Designing and Data Gathering
After identifying the factors, a statistical survey is performed to determine the crisp weight of each factor. For this purpose, a questionnaire is designed. The Cronbach's alpha values showed the reliability of the questions before gathering data from the real-world, which consisted of various automotive, house appliance, and metal forming industries that use job-shop facility locations. Table 1 shows Cronbach's Alpha:

Pearson Correlation Test
The correlations between factors is then calculated using the Pearson Correlation Test. Table 2 indicates the Pearson Correlation Test results. The outcomes indicate that positive correlations exist between most of the factors. In the next step, those factors with strong correlations (p ≥ 0.6) and whether the correlations are meaningful in the real world will be considered in the regression formula developed in the next section.

Determining the Crisp Weights of the Variables
After determining the correlations between factors, the crisp weights for the variables (factors) are determined. For this purpose, a regression equation is developed by SPSS. The coefficients of variables are shown in Table 3: As shown by Equation (1), various variables impact efficient scheduling differently. In this regard, the part route, as question 6 found, is at the top of the list of influential factors that can influence scheduling efficiency. In addition, the statistical society believed that Workers' Skills and Production Technology are the other most important factors that can significantly impact the scheduling of job-shop-based manufacturing systems. Meanwhile, they also believed that Machine Amortization and Raw Material Quality are vital factors that should not be ignored during the scheduling process. By contrast, it is found that Product Formula and Procedures have fewer impacts on the dependent variable.

Fuzzy Inference System (Phase 2)
The weight factors calculated in the last section were crisp, meaning that the uncertainty has not been considered for them. However, in reality, the factors in a dynamic condition of the manufacturing system are not crisp and may be affected by uncertainty. For this purpose, in this section, a fuzzy inference system will be developed to involve the uncertainty level of the crisp factors, thus making them fuzzy weights. The fuzzy weights will be directly used in the mathematical model presented in the next section.

Fuzzy Inference System Model
In this research, the responders' confidence level will involve uncertainty. For this purpose, the level of certainty was asked of responders after each question. This approach is inspired by Alazemi et al. [44]. For this purpose, the responders selected the level of confidence from a range of uncertainties for each factor. Then, the FIS model, where a multi-input single-output (MISO) design is used for each factor, is developed as follows (     Table 4 indicates the designed FIS model according to the structure explained in Figure 3.

Fuzzy Domain Function
Trapezoidal Distribution Functions Fuzzy Rule 9 Rules

Fuzzy Function
Mamdani rule:

If (FP is i) and (FP is j) then (FO is k)
where FP can be considered as the first variable's input and FP as the second variable's input. FO will be a fuzzy output accordingly.  The trapezoidal distribution function is a commonly used function for expressing fuzzy ranges for a data series by considering breaking points to identify the range of data. It is selected in this article based on the scoring range in Section 3.3. Other distribution functions may be considered, such as the triangular distribution function.
The model design and fuzzy membership functions are shown in Figure 5. The model is then expanded and represented in Figure 6. The trapezoidal distribution function is a commonly used function for expressing fuzzy ranges for a data series by considering breaking points to identify the range of data. It is selected in this article based on the scoring range in Section 3.3. Other distribution functions may be considered, such as the triangular distribution function.
The model design and fuzzy membership functions are shown in Figure 5. The model is then expanded and represented in Figure 6. The trapezoidal distribution function is a commonly used function for expressing fuzzy ranges for a data series by considering breaking points to identify the range of data. It is selected in this article based on the scoring range in Section 3.3. Other distribution functions may be considered, such as the triangular distribution function.
The model design and fuzzy membership functions are shown in Figure 5. The model is then expanded and represented in Figure 6.         The fuzzy weights can be calculated using Matlab's fuzzy weights output window. Table 5 indicates the results of the fuzzy weights. The weights are then normalized to be applied to the mathematical model in the next section.  The fuzzy weights can be calculated using Matlab's fuzzy weights output window. Table 5 indicates the results of the fuzzy weights. The weights are then normalized to be applied to the mathematical model in the next section.  Figure 9 indicates the 3D graphical view of the FIS model, which shows the areas where uncertainty can affect the risk factors. The fuzzy weights can be calculated using Matlab's fuzzy weights output window. Table 5 indicates the results of the fuzzy weights. The weights are then normalized to be applied to the mathematical model in the next section.  Figure 9 indicates the 3D graphical view of the FIS model, which shows the areas where uncertainty can affect the risk factors.
The above risks can then be used as the fuzzy weights of the related objective functions of the mathematical model in the next section.

Mathematical Modeling (Phase 3)
This section will construct a new mathematical model with fuzzy weights in which each risk factor is considered a variable. The mathematical model's parameters will then be described.

Indexes
The indexes of the model will be explained as follow: i: number of products (3) P: number of machines (4) q: transportation counter (5) k: number of shops, each shops has series of parallel machines (6) but with difference capacities (7) t: time slot (8) r: material type (9)

Parameters
The parameters of the model will be explained as follow: TH: Time horizon in a manufacturing period (10) τ pk : The time required to perform an operation by the machine p in the shop k (11) λ qk : The time required to transporting one product using transportation type q (12) D it : Demand of product type i in period t (13) OPC ik : Operation process chart of product (if product I needs service type k) C pk : Capacity of machine p in shop k (15) DP PK : Defect rate of machine p in shop kth (16) Dr r : Defect rate of raw material type r (17) α i,r : The usage coefficient of raw material type r for producing one unit of product i

Decision Variables
The decision variables of the model will be explained as follow:

Mathematical Model
In this section, the mathematical model will be developed as follow: Min Z 1 : ω 1 .
The first objective function is to minimize completion time (Equation (26)). This objective includes two main parts: in the first part, the operational time is minimized, and in the second part, transportation time is minimized. The second objective targets operational and transportation costs within and between shops to complete all essential services (Equation (27)). The third objective tries to minimize the defect rate of the products by considering the raw material defect rate and machine defect rate (Equation (28)). The first constraint guarantees that in each production period, the market demands must be fulfilled (Equation (29)). For this purpose, the number of total operations are divided into the total number of required operations for each product. For example, if the number of total operations for product 1 is determined as 300 (∑ K k ∑ P p X i,p,k,t ) in a manufacturing period, and for completing each unit of product type one, five operations are suggested in the OPC matrix ∑ K k ∑ P p OPC ik , then the total number of product type 1 in the period will be 60 (∑ K k ∑ P p X i,p,k,t /OPC ik ). The second constraint is ensuring that the number of operations for all products in a production period will not exceed the machine's capacity (Equation (30)). The third constraint supports the time horizon in each manufacturing period by ensuring that the total times for production and transferring of all product types will not exceed the available time horizon in that period (Equation (31)). The fourth constraint shows that the total amount of each raw material type used for manufacturing the products will not exceed the total available raw material in that period (Equation (32)). The fuzzy weights (ω i ) extracted in the previous part are used here as the weights of the objective functions (Equation (33)). Table 6 compares the contribution and novelty of this study to those of other similar studies in the literature. Focuses on the different meta-heuristics in DRCSPs.

Weighted NSGA-II (Phase 4)
The developed model in phase 3 is multi-objective. Usually, multi-objective models are complex, and there is no exact optimum point that can optimize all objectives simultaneously. For this purpose, in most cases, scientists used heuristics and metaheuristics to solve multi-objective models [8,[56][57][58]. However, the classic NSGA-II does not consider the uncertainty variables (factors) in its calculation. As a result, to address this issue, a new fuzzy-weighted NSGA-II will be proposed in this research to solve the developed model. NSGA-II was firstly proposed by Srinivas et al. [59] and has been widely used to solve multi-objective models many times [60][61][62]. Verma et al. provided a novel review on using NSGA-II for combinational optimization problems [63]. Recently, fuzzy NSGA-II has attracted the attention of scientists due to its high capability to consider uncertainty in multi-objective optimization problems [64,65]. This research proposes a new fuzzy weight NSGA-II to solve the multi-objective job-hop-based scheduling problem in a fuzzy environment. The flowchart of the proposed fuzzy-weighted NSGA-II is shown in Figure 10. Then, the pseudo-code for the FW-NSGA-II is shown in Algorithm 1.
The algorithm starts by entering algorithm parameters (including population size, generations, and mutation rate) along with the case study parameters (including the number of products, machines, raw material type, shops, and demand). According to OPC, the algorithm then recognizes the possible list of solutions for the machines that can serve the required operations consecutively. At this stage, using the crossover operator, the part routes for the possible solutions will be determined (Equation (41)). Then, a list of possible solutions, called a part route list, will be further processed by calculating fuzzy weighted objective functions. In continuation, those solutions that achieved greater multi-function objective function will be stored in the tournament list for the next generations. In the next generations, new solutions will be generated by changing the possible part routes using parallel machines in the same shop. The process will be repeated until the maximum number of generations is achieved (stop criterion).

Cross-Over Operator
The cross-over operator for the proposed method will be explained as follow: For i ∈ I For p ∈ OPC i List the required machines for the operation i Find potential machines in the related shop If C pk < C pk : Remove p from the machine list: Select a machine from the machine list Generate the new part route list * C pk is the remained capacity of a machine 3.9.2. Setting NSGA-II Parameters Before using the proposed NSGA-II for the solving algorithms, the algorithm parameters such as the number of generations, population size, cross-over rate, mutation rate, and neighborhood ratio should be set. The aim is to tune up the NSGA-II to achieve maximum efficiency. For this purpose, a Taguchi model is designed using Minitab ® 18.0. The outcomes of the best values that were estimated for each of the solving algorithm parameters are shown in Table 7: Algorithm 1 Pseudo code of the proposed FW-NSGA-II algorithm)

Insert Dataset
Input data sets (TH; S; τ pk ; D it ; OPC ik ; C pk ; Q PK ; CO ipk ; M k ) Input parameters (population size, mutation rate, generations) Input Fuzzy Weights (ω1, ω2, ω3) For Generation 1:1:N: for Population 1:1: pop.size: If t is 1: Find appropriate machines that can serve the first operations Else Choose a member from the candidate list List the required consecutive operations to complete a product (OPC iK ) Find tournament list i Activate Cross-over Operator Choose the best part route Call Mutation Operator Rand a number (Y : Rand (1) Min Z 1 : ω 1 .    The Lp-metric method is a well-known and robust technique for solving multiobjective algorithms. The Lp-metric algorithm starts by solving models and their routes (p = 2, 3, etc.). This strategy enables the algorithm to find solutions with higher quality. The algorithm stops whenever the solutions of two consecutive routes are the same (for example, the results of solving the second and third routes of the model became equal).
In continuation, similar to the previous algorithm, the algorithm is applied for a smallsize case study. The results showed that the solver (Branch and Bound Algorithm) stopped in iteration 22 by falling into a local optimum point. Table 8 shows that the Lp-metric algorithm found the same (and best) solutions in the third route (p = 3). A weighting method coded with Lingo 12.0 ® is used to solve the multi-objective job-shop-based scheduling model that enables us to compare the results with the fuzzyweighted NSGA-II. The algorithm is then applied to a small case study. Lingo 12.0 can easily recognize the model as mixed integer programming (MIP) and find the global optimum solution in 31 iterations. The lingo used the branch and bound method to solve the model using the multi-objective weighting technique.
The algorithm is then applied to seven data sets generated considering different data sets. The outcomes in Table 9 show that the method could find global optimum solutions for small-scale problems. However, it could only find the local optimum for large-scale problems. In addition, the elapsed time for solving the case studies increased dramatically by increasing the size of the variables.

Results and Discussion
In order to check the performance of the proposed method, a number of case studies on a small, medium, and large scale will be solved. The outcomes' performance will then be examined using some metrics (Table 10). Note that all case studies were solved by a personal laptop supported by Intel(R) Core(TM) i7, CPU @ 2.40 GHz, and 8 GB RAM.
The findings showed that the proposed FW-NSGA-II could solve all the case studies in a reasonable time. Then, the model is solved by other methods, namely the Sine Cosine Algorithm, Simulated Annealing, Tabu Search, and TLBO heuristic. The results showed that, in many cases, the proposed method could provide better solutions (Table 11).

Verification of the Model by Application to a Real Manufacturing Firm
In the next step, the proposed FW-NSGA-II was applied to a medium-scale metalforming company located in Malaysia. The manufacturing plant has four workshops to produce various products. The firm contains cutting, bending, welding, and painting shops. Figure 11 shows the workshop view: In the next step, the proposed FW-NSGA-II was applied to a medium-scale metalforming company located in Malaysia. The manufacturing plant has four workshops to produce various products. The firm contains cutting, bending, welding, and painting shops. Figure 11 shows the workshop view: Figure 11. Graphical view of the case study. * The numbers in the parenthesis show the quality ratio of outcomes by using a machine.
The following table shows the OPC products (Table 12). In Table 12, the operation sequence of the products is displayed. Intersections in the table with a zero number indicate that the product is not needed.
The demand for products for the three planning periods is as follows (Table 13): Figure 11. Graphical view of the case study. * The numbers in the parenthesis show the quality ratio of outcomes by using a machine.
The following table shows the OPC products (Table 12). In Table 12, the operation sequence of the products is displayed. Intersections in the table with a zero number indicate that the product is not needed.
The demand for products for the three planning periods is as follows (Table 13): In addition, the amount of time required (minutes) for each operation is listed in Table 14. Then, the best combination of production and sequence of operations for such a combination is found. The quality coefficient of each machine in the given layout is specified in blue and is given in Table 15. The capacity of each manufactured machine is shown below:

Managerial and Practical Implications
It was found that strong correlations exist between the total completion time, product defect costs, human resources, and operational costs that should not be ignored during the scheduling of a dual-resource-constrained manufacturing system. In addition, it was shown that the uncertain environment has a strong effect on the crisp weights of the studied factors, changing the completion time, human resource and operations, and material risk factors from 0.413, 0.366, and 0.221 to 0.582, 0.203, and 0.215, respectively.
Moreover, minimizing the total objective function consisting of completion time, cost, and product defect to 11.2% using the proposed method can be considered another implication of using the proposed method.

Technological Opportunities
The proposed framework can be used as the core of scheduling software to improve scheduling efficiency by considering factors' real values in a DRC manufacturing system. For this purpose, a cloud-based system can be designed to obtain real-time data from the manufacturing system and forward it to a scheduling software where the proposed solving method can be used to schedule the system.

Conclusions
This paper presents a 4-phased fuzzy framework for the multi-objective scheduling of job-shop-based manufacturing companies where product completion time, operational cost, and product defect cost were considered as the objective functions. For this purpose, an NL-IP mathematical model is developed, which is then solved by a proposed fuzzy-weighted NSGA-II method for seven case studies.
In the first phase of the framework, the influential factors in scheduling manufacturing systems in the fuzzy environment were identified. Then, using a fuzzy inference system, the fuzzy weight of the factors was determined. A multi-objective mathematical model was developed in the next phase to optimize the product completion time, operational cost, and product defect costs using the determined fuzzy weights. In the last phase, a fuzzy-weighted NSGA-II was proposed to solve the fuzzy multi-objective scheduling of job-shop-based manufacturing companies.
The outcomes of the research are threefold. First, it was demonstrated that the unpredictability of the environment had a significant impact on the crisp weights of the factors under study, causing the completion time, operational, and material risk factors weights to change from 0.413, 0.366, and 0.221 to 0.582, 0.203, and 0.215, respectively. Then, it was found that while using the proposed framework, the total objective function is boosted up to 5% for small-scale, 11.2% for medium-scale, and 3.8% for large-scale studied cases. Second, it was found that there are interactions between the studied objective functions, wherein minimizing the total completion time and the product defect costs have a strong adverse correlation (up to 7.5% for the studied cases). Additionally, there was a negative interaction between product completion time and human resource and operational costs (up to 14.3% for the studied cases). Finally, the correlations between product defect and human resource and operational costs had super strong correlations of 23.2%.
Third, although the proposed weighted NSGA-II method was slower in terms of solving time, at the same time, it showed superiority when applied to the fuzzy multiobjective scheduling of job-shop-based manufacturing companies.
Further expansion of the proposed framework can be carried out by considering the external factors influencing a manufacturing system. Moreover, the integration of the proposed model with real-time data application is suggested for future research.