Optimal Allocation Method of Discrete Manufacturing Resources for Demand Coordination between Suppliers and Customers in a Fuzzy Environment

Discrete manufacturing products are often assembled from multiple parts through a series of discrete processes. How to effectively configure resources in a discrete manufacturing environment is an important research topic worthy of attention. Based on an in-depth analysis of the discrete manufacturing operation model and the manufacturing resource allocation process, this paper fully considers the uncertainty factors of the manufacturing resource customers and the interests of the manufacturing resource suppliers and proposes a bilevel planning model under a fuzzy environment that comprehensively considers the customers’ expectation bias and the suppliers’ profit maximization. The method firstly uses a language phrase to collect the language evaluation of the customers and suppliers for manufacturing tasks and uses a trapezoidal fuzzy number to convert the language evaluation phrase into a value that can be calculated. Then, we use the prospect theory to optimize the constraint indicators based on the language evaluation of customers and suppliers. Next, the bilevel planning model for optimal configuration of manufacturing resources in discrete manufacturing environment is established under the consideration of the respective interests of both the customers and the suppliers, and the fast nondominated sorting genetic algorithm (NSGA-II) is used to solve the model. Finally, an example is given to verify the validity and feasibility of the model.


Introduction
With the rapid economic growth, customer demand for products has become more diverse.How to effectively grasp customer demand, shorten lead time, lower production cost, and increase product quality are key factors for companies to achieve sustainable development [1].Due to unpredictable market changes, this requires the manufacturing systems to be able to rapidly reconstruct in response to rapid market changes.Discrete manufacturing has gradually become the mainstream model of manufacturing industry because of its advantages such as noncontinuity and reconfigurability [2].With lots of manufacturing units that provide the same functionality but have different parameters in a discrete system, effective resource configuration which can reflect customer needs is often considered as a key technology [3].
Manufacturing resource configuration (MRC) plays a very important role in discrete systems, especially when manufacturing systems have to cope with shorter product life cycles [4].In order to fulfill the dynamic customer needs, it always needs discrete systems to invoke several manufacturing units in sequence and combine them together fast.Due to complexity and diversity of manufacturing resources, resource optimal configuration has become a key issue in discrete systems and has been widely studied in both industrial community and academia.Despite of significant progress achieved by the researchers in manufacturing resource configuration, grey relational analysis [5], manufacturing grids [6], the idea of Pareto [7], graph theoretic methods [8], artificial intelligence-based methods [9], and other methods have been proposed successively.Most of current MRC algorithms where the data of MRC are in the form of real number are not suitable for discrete systems environment because the MRC of discrete systems is often fuzzy and uncertain.For example, when MRD describe the quality of products, they can better express their customers' perception by using language evaluation words of "good," "bad," and so on.Therefore, this paper proposes a resource optimization configuration that considers the customers' expectation bias and considers the suppliers' profit maximization under fuzzy environment.
The remainder of this paper is organized as follows.After reviewing the related literature in Literature Reviews, some basic theories such as prospect theory and bilevel programming model are introduced in Preliminary Knowledge.Problem Description and Symbol Introduction presents the issues of the thesis research and some basic symbols.In The Proposed Method, it proposes a novel method for discrete manufacturing resource optimization configuration under fuzzy environment.An illustrative example is given to demonstrate the use of the proposed method in Illustrative Example.Conclusions and discussions are drawn with brief comments in Conclusions and Discussion.

Literature Reviews
MRC has emerged because of the need for manufacturing organizations to cope with shorter product life cycles, timeto-market, and a shift to respond to demands for MRD [10].In the past years, many approaches, models, and methodologies have been proposed for solving manufacturing resource configuration problems.Among them, the main research algorithms are as follows: (i) Grey relational analysis (ii) Manufacturing grid (iii) The idea of Pareto (iv) Graph theoretic methods (v) Artificial intelligence-based methods Zhang used grey relational analysis to further study manufacturing machine and manufacturing cell (MC) of multigranularity resource configuration process.During resource modeling, advanced information and sensor technologies are adopted to construct the information models of resources, which make the traditional production process more transparent, traceable, and on-line controllable [5].Based on the quantum evolution theory, Zhang and Hu proposed a hybrid chaotic quantum evolutionary algorithm (CQEA) for resource combinatorial optimization (RCO) problems.Using an example to prove the proposed CQEA is effective, efficient, and scalable for the RCO problem in manufacturing grid system [6].Li et al. gave a resource configuration method based on binary decision diagram (BDD) which is a directed acyclic graph (DAG) based on Shannon's decomposition.This method extends the scale of the reliability system.Through the results of three case studies, it is found that the decision graph expansion method is more computationally efficient than the traditional BDD [7].Xiang et al. introduce a new multiobjective optimization algorithm based on the combination of the idea of Pareto solution and group leader algorithm (GLA), which study of quality of service (QoS) and energy consumption assessment (EnCon) [8].Tao et al. gave a parallel intelligent algorithm of resource configuration, which can minimize implementation time and cost and maximize the reliability of MGrid resource service composition paths [9].
Although there are many research results on the optimal configuration of resources, there are two obvious shortcomings in the current research results.First of all, the existing research results can promote and facilitate quicker and smarter decisions for service composition, but it cannot play any role in the fuzzy problem in the discrete system.That is, most of the existing research methods use real number, but the evaluation of manufacturing resources in discrete systems is often fuzzy and uncertain.Secondly, most of the existing research results are single objective, which can only consider the interests of one of the MRD or the MRP, and there are a few multiobjective research results, which consider the interests of both the MRD and the MRP, but when modeling, the multiobjective is converted into a single objective by using the weighted operator, and it is not truly multiobjective.
This paper proposes a bilevel programming model considering customer expectation under a fuzzy environment.Based on the prospect theory, we calculate the expected deviations of cost, quality, time, and green indicators of manufacturing resource customers and establish the objective function, then considering the efficiency, coordination, agility indicators of the manufacturing resource suppliers, we establish the objective function.Based on bilevel programming theory, the objective function and constraint conditions of MRD and MRP are established, respectively, and the fast nondominated sorting genetic algorithm (NSGA-II) is used to solve the model.

Trapezoidal Fuzzy Number and Language Evaluation Phrases
3.1.1.Trapezoidal Fuzzy Numbers.Let A be a fuzzy set, and a fuzzy subset a of A is defined with a membership function v a x that maps each element x in a to a real number in the interval [0, 1].The function value of v a x signifies the grade of membership of x in a [11].A trapezoidal fuzzy number a represented with four points as follows: a = a 1 , a 2 , a 3 , a 4 , a 1 ≤ a 2 ≤ a 3 ≤ a 4 (see Figure 1).Its membership function v a x is defined as Complexity where a 1 , a 2 , a 3 , and a 4 are real numbers and these constants reflect the fuzziness of the evaluation data [12].As shown in Figure 1, the trapezoidal fuzzy numbers can be denoted by a 1 , a 2 , a 3 , a 4 .The x in interval a 2 , a 3 gives the maximal grade of v a x , that is, v a x = 1, and it is the most probable value of the evaluation data.a 1 and a 4 are the lower and upper limits of the available area for the evaluation data, and they reflect the fuzziness of the evaluation data.If a 2 = a 3 , then a = a 1 , a 2 , a 3 , a 4 is reduced to a triangular fuzzy number a = a 1 , a M , a 4 , where When a > 0 and b > 0, the basic arithmetic operations are as follows [13]: Definition 1 (see [14]).There are two trapezoidal fuzzy numbers a = a 1 , a 2 , a 3 , a 4  Definition 2 (see [14]).There are two trapezoidal fuzzy numbers a = a 1 , a 2 , a 3 , a 3.1.2.Language Evaluation Phrases.Let G = g 0 , g 1 , … , g l be the preestablished finite and totally ordered linguistic term set with odd cardinalities, where g i denotes the ith linguistic term of set G, and l + 1 is the cardinality of G.The middle term in linguistic term set G is thought to be represented an assessment of "approximately 0.5" and the remaining terms of G are thought to be placed around it symmetrically [15].
For example, a linguistic term set with seven terms can be expressed as follows.
Therefore, g i can be approximately expressed as a triangular fuzzy number a = a 1 , a 2 , a 3 , a 4 using the following formula [16]: For example, by (3), each linguistic term in a linguistic term set with seven terms, G = g 0 , g 1 , g 2 , g 3 , g 4 , g 5 , g 6 , can be expressed as the corresponding trapezoidal fuzzy numbers listed in Table 1, where l = 6 and i = 0, 1, 2, … , 6.
3.2.Prospect Theory.Tversky and Kahneman believe that the actual decision-making behavior of an individual under uncertainty is deviated from the basic principle of expected utility.The uncertainty foreground f is a function from the natural state set S to the result set X, that is, f S → X.The foreground f is a sequence of x i , S i sequences, where S i is a division of S which is called S i event.When event S i occurs, it produces a result of x i .If i > j, then x i > x j .The result x i of each foreground is arranged in ascending order, that is, If you choose x 0 as the reference point and its value is "0," the profit is x i > 0 and the loss is x i < 0.
In 1992, Tversky and Kahneman proposed the cumulative prospect theory [17] based on the prospect theory [18].Cumulative prospect theory has introduced capacity theory, which can better solve the problem of dominant advantage and deal with problems with multiple results.The capacity value can be expressed as w + and w − , then the value of the foreground f = x i , S i , −m ≤ i ≤ n is expressed as where V f + and V f − can be obtained by ( 5) and ( 6), respectively.
Figure 1: A triangular fuzzy number a. 3 Complexity Here, If the result of event S i is probabilistic, the decision problem can be seen as a probabilistic prospect, that is, f = x i , S i → f = x i , p i .Here, p S i = p i represents the probability of occurrence of event S i .Under these circumstances, The core content of the cumulative prospect theory is the value function v and the weight function p, which can be expressed as The Weight Function.The weight function converts the probability into the decision weight, so the calculation formulas for the probability weight of the profit and loss are where p is the probability; γ and δ are that parameters that indicate the degree of curvature of the probability weight function.

The Value Function.
A great breakthrough in the expectation theory is to replace the traditional utility function with a value function, so that the carriers concerned can be implemented in the value change rather than the final amount [19].The value function is to convert surface value into decision value.The specific form of the value function of Tversky and Kahneman is [17] Here, when x ≥ 0, v x indicates profit; when x < 0, v x indicates loss.α and β indicate the degree of roughness of the value function in the region of profit and loss, that is, the rate of decline in the sensitivity of decision-makers, 0 < α < 1 and 0 < β < 1 [20].λ indicates that the loss area of the value function is steeper than the income area, that is, it reflects the degree of loss avoidance of decision-makers.If λ > 1, it denotes that the decision-maker is more sensitive to the loss.People tend to risk gambling when faced with conditions of considerable loss but tend to accept certainty profits when faced with fairly favorable earnings.The happiness caused by profit is not equal to that caused by the same amount of loss.The latter is greater than the former [21].
When people evaluate a thing or make a choice, they always intentionally or unintentionally compare it with a certain reference, which is called the reference point from the definition of mathematics.The reference point is a very important feature of the value function, because the profit or loss is always compared with a certain reference point.The value of the profit or loss from the reference point is located on the right side of the profit or loss to indicate positive evaluation, while the left side of the profit or loss indicates negative evaluation [22].The reference point is used as an evaluation criterion.It is subjectively determined by an individual and will change due to different evaluation topics, environmental time, and the like.
This paper selects the expectations of MRD as a reference point for the value function.At the same time, we divide the decision indicators into benefit type and cost type.The profitability index means that the larger the index value, the better; and the cost index means that the smaller the index value, the better.According to (7) and Definition 2, the value function of the profitability index and cost index is, respectively, where θ r ik is the reference point (i.e., expectation) given by the MRD B k for the r index of the manufacturing task TS i ; h r jk indicates that the r index of the selected manufacturing unit MS ij gives the manufacturing resource customers' perception after completing the manufacturing task.When ijk are referred to as the value of the profit generated by the MAD relative to the reference point θ r ik ; on the contrary, when ijk are referred to as the value of the loss generated by the MAD relative to the reference point θ r ik ; α k and β k are the degree of roughness of the profit area and the loss area of the value function, and λ k is the loss aversion coefficient [20,21].In actual decision analysis, these parameters are usually obtained by nonlinear regression of experimental data [23].

The Bilevel Programming.
Manufacturing unit configuration is the main form of discrete manufacturing tasks.The process of optimization selection belongs to the typical multiobjective optimization problems (MOP).The traditional MOP solution is to convert MOP into single-object problem solving.Common methods include the main target method, linear weighting method, and hierarchical optimization method.But the discrete manufacturing unit configuration 4 Complexity optimization process involves the interests of the customers and suppliers, because the customers and suppliers each have a part of the optimization variables, each represents their own interests, and each variable exists between influences and constraints, and it cannot be solved using traditional methods.The bilevel programming [24] model is a hierarchical model with a master-slave hierarchical structure.In the bilevel programming model, the upper and lower decisionmakers have their own objective functions and constraints.
The upper decision-maker makes decisions firstly.According to their objective functions and constraints, the lower decision-maker obtains the optimal solutions within the possible range and sends back their own optimal solutions to the upper decision-maker.Then, the upper decision-maker obtains the global optimal solution of the problem within the possible scope based on the optimal solutions of the lower decision-maker.The mathematical description of the bilevel programming model is as follows: where (U) is the upper plan and (L) is the lower plan; Z is the objective function of the upper plan, x is the determinate variable of the upper plan, and G is the constraint condition of the decision variable x; z is the objective function of the lower plan, y is the decision variable of the lower plan, g is the constraint condition of the decision variable y, and the lower decision variable y is the function of the upper level decision variable x, that is, y = y(x).
The upper and lower optimization problems are relatively independent, and their optimization processes are dependent on each other, so the bilevel programming problems cannot usually be independently solved layer by layer.Not only the interests of MRD but also the interests of MRP must be taken into account in a discrete manufacturing environment.The idea of bilevel programming is used to solve the problem of optimal configuration of manufacturing unit in discrete manufacturing environments, as shown in Figure 2.

Problem Description and Symbol Introduction
The optimal configuration of manufacturing resources in discrete environments refers to the process of configuring reasonable manufacturing resources or manufacturing resource combinations according to different manufacturing tasks [5].After a discrete manufacturing task is issued, the MRP decomposes the total manufacturing task into a set of manufacturing subtasks according to the decomposition preferences.According to the matching principle of manufacturing resources, we find all manufacturing units that can meet the manufacturing subtasks and form a set of manufacturing units.Considering the profit of MRD and MRP, we select the optimal set of manufacturing units to form the optimal manufacturing resource configuration, as shown in Figure 3.With the emergence and development of networked manufacturing technologies, manufacturing units are not only confined to the internal structure of a single company but are also composed of multiple companies in different geographic locations, so the manufacturing resource optimization configurations are distribution, heterogeneity, and dynamics.Compared with the existing method for optimal configuration of manufacturing resources, this paper proposes an optimization algorithm for manufacturing resources in a fuzzy environment that considers the interests of both MRD and MRP.Firstly, the MRP decomposes the MRD's product demand for the MRP (i.e., the total manufacturing task) according to design preferences and obtains a set of manufacturing subtasks.Secondly, we use the existing infrastructure to match the manufacturing subtasks with the existing manufacturing units.A set of manufacturing units for each subtask is obtained.Next, a questionnaire is used to obtain the MRP and MRD's linguistic evaluation of each manufacturing unit, and u a trapezoidal fuzzy number is used to translate the language phrase into numerical values.Then, we obtain the customer's expected deviation value for each manufacturing unit based on the prospect theory.On this basis, a bilevel planning model that considers the (ii) TS = TS 1 , TS 2 , … , TS n is a set of manufacturing subtasks, where TS i denotes the ith manufacturing subtasks, i = 1, … , n.
(iii) MS i = MS i1 , MS i2 , … , MS im i is a candidate manufacturing unit set of the ith manufacturing subtasks, which can provide a similar function.
Where MS ij denotes the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and (vi (vii) G i = G i1 , G i2 , … , G im i is a green set of the ith manufacturing subtasks, which is the set of environmental protection, where G ij denotes green of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .
(viii) E i = E i1 , E i2 , … , E im i is an efficiency set of the ith manufacturing subtasks, where E ij denotes efficiency of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .
(ix) F i = F i1 , F i2 , … , F im i is an agility set of the ith manufacturing subtasks, where F ij denotes agility of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .
is a coordination set of the ith manufacturing subtasks, where R ij denotes coordination of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .6 Complexity greater than the joy after harvest.Therefore, the prospect theory can be used to indicate people's satisfaction with products or equipment.The prospect theory focuses on the difference in customer psychology, not just the customer's final value.In this paper, MRD expectation is used as a reference point to optimize the MRD indexes and MRP indexes.

Index Constraint Optimization of MRD
(1) Cost Constraint.In a discrete manufacturing system, the cost indicators of each manufacturing unit include not only the fixed costs of manufacturing units but also the logistics costs between manufacturing units.The optimization of the cost constraint index is based on the prospect value of the MRD for the cost index.It can be calculated as Here, V C denotes the total cost prospect value of manufacturing units selected by MRD; H ij can be expressed as The traditional cost constraint is that the cost of completing the task for the selected manufacturing units' combination cannot be greater than the maximum cost required by the MRD.On the basis of the prospect theory, the traditional cost constraint conditions are transformed into that the combined cost prospect value of the selected manufacturing units cannot be greater than the maximum cost prospect value given by the MRD, which can be expressed as (2) Time Constraint.In discrete manufacturing systems, the time of each manufacturing unit mainly refers to the total time from which the raw material to workshop to finished product leaving workshop.Each manufacturing unit time involves the time of the unit running and the product twisting time between the manufacturing units.The optimization of the time constraint index is based on the prospect value of the MRD for the time index.It can be calculated as Here, V T denotes the total time prospect value of manufacturing units selected by MRD; H ij can be expressed as The traditional time constraint is that the time required for the selected manufacturing unit combination to complete the task cannot be greater than the longest delivery time required by the MRD.On the basis of the prospect theory, the traditional time constraint conditions are transformed into that the combined time prospect value of the selected manufacturing units cannot be greater than the maximum time prospect value given by the MRD.

V T ≤ V T max 19
(3) Quality Constraint.In the discrete manufacturing system, the quality index of each manufacturing unit refers to the quality of each manufacturing unit to complete the relevant manufacturing tasks, that is, the quality qualification rate for each manufacturing unit to complete the relevant manufacturing tasks.The optimization of the quality constraint index is based on the prospect value of the MRD for the quality index.It can be calculated as 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 15 H ij = 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 18 7 Complexity where V T denotes the total quality prospect value of manufacturing units selected by MRD; H ij can be expressed as The traditional quality constraint is that the quality qualification rate of any manufacturing unit in the selected manufacturing units' combination should not be less than the minimum quality qualification rate required by the MDR.On the basis of the prospect theory, the traditional quality constraint conditions are transformed into that the quality prospect value of any manufacturing unit in the selected manufacturing unit combination should not be less than the minimum quality prospect value required by the MRD, which can be expressed as Green is based on the requirements of the current low carbon environmental protection.The evaluation of the environmental indicators includes carbon emissions and processing material losses in the manufacturing process of manufacturing units.In discrete manufacturing systems, the green index of each manufacturing unit is the degree of pollution to the environment by which the selected manufacturing units are combined to complete the related manufacturing tasks.The optimization of the green constraint index is based on the prospect value of the MRD for the green index.It can be calculated as where V G denotes the total green prospect value of manufacturing units selected by MRD; H ij can be expressed as The traditional environmental protection constraint is that the pollution rate of any manufacturing unit in the selected manufacturing unit combination cannot be greater than the highest manufacturing pollution rate required by the MRD.On the basis of the prospect theory, the traditional environmental constraints are translated into that the prospect value of pollution of any manufacturing unit in the selected combination of manufacturing units cannot be greater than the maximum prospect value of pollution required by the MRD, which can be expressed as 5.1.2.Index Constraint Optimization of MRP.Manufacturing tasks in discrete manufacturing systems are susceptible to changes in manufacturing tasks and changes in manufacturing resources, making it impossible to complete manufacturing tasks with high quality and efficiency.Therefore, it is very important for manufacturing tasks to consider the factors in the process of optimizing the configuration of manufacturing resources.
(1) The Efficiency of the Manufacturing Units' Configuration.
The efficiency of the configuration of manufacturing units in a discrete manufacturing system refers to the number of qualified products that the manufacturing unit produces within a unit time after the manufacturing task is reached.The efficiency index of manufacturing units' configuration mainly includes the functional efficiency of each manufacturing unit E n ij and the decomposition capability of the manufacturing unit E f ij , which can be expressed as 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 21 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 24 8 Complexity Here, The functional efficiency E n ij of any manufacturing units cannot be less than the minimum functional efficiency E n min required by MRP; and the decomposition capability of the manufacturing unit E f ij of any MRP in the manufacturing unit's configuration cannot be less than the minimum decomposition capability of the manufacturing unit E f min required by MRP.They can be expressed as

28
(2) The Agility of the Manufacturing Units' Configuration.
The agility of the manufacturing unit configuration in a discrete manufacturing environment refers to the manufacturing units' configuration ability to react quickly and successfully complete the manufacturing task when the content of the manufacturing task changes or a manufacturing resource withdraws for some reason.So that enterprises can cope with the rapidly changing and unpredictable market demand, thus obtaining long-term economic benefits enterprises.The agility index of manufacturing units' configuration mainly includes the functional diversity F d ij and the manufacturing resource types F z ij of each manufacturing unit, which can be expressed as where The functional diversity F d ij of any manufacturing unit cannot be less than the minimum functional diversity F n min required by the MRP; the type of manufacturing resources F z ij provided by any MRP in the manufacturing units' configuration should not be less than the minimum type of manufacturing resource F z min required by the MRP.They can be expressed as 3) The Coordination of the Manufacturing Units' Configuration.The coordination of the manufacturing units' configuration in discrete manufacturing environment refers to manufacturing units that can coordinate and efficiently complete manufacturing task when the manufacturing task comes down.The coordination index of manufacturing units' configuration mainly includes the reliability of the manufacturing unit R k ij and the compatibility of the manufacturing unit R x ij , which can be expressed as

32
where H ij = 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 27 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 30 1, the ith manuf acturing subtask is completed by the manuf acturing unit MS ij , 0, the ith manuf acturing subtask is not completed by the manuf acturing unit MS ij 33

Complexity
The reliability of any manufacturing unit R k ij in the manufacturing units' configuration cannot be less than the minimum reliability R k min required by the MRP.The coordination of any manufacturing unit R x ij in the manufacturing unit combination cannot be less than the minimum coordination R x min required by the MRP.They can be expressed as  In discrete manufacturing, it is difficult for MRD and MRP to express their preference for a certain characteristic with an explicit or implicit "utility function," but it is generally easy to determine what level of expectation a certain target achieves and according to the actual situation to adjust the level of the object expected to reach.In these elements that influence multiobjective decision-making in discrete manufacturing environment, they may contain both integer variables and fractional variables and may be continuous or discontinuous [25].
To solve this kind of decision-making problem, this paper proposes a bilevel programming decision method based on the object expectation: selecting index cost (C), time (T), quality (Q), and green (G) as the upper level optimization objectives that affect the MRD of manufacturing resources and then, selecting the index efficiency (E), agility (F), and coordination(R) as the lower level optimization objectives that affect the MRP of manufacturing resources.The bilevel optimization model under the discrete manufacturing environment is established as follows: Here, (35) denotes the upper-level optimization object given by the MRD, and (36) and (39) denote the quality, cost, time, and green indicator constraints of the upper-level optimization objectives; (40) denotes the lower-level optimization object given by the MRP, and (41) and ( 46  11 Complexity efficiency, coordination, and agility indicator constraints of the lower-level optimization object.

Model Solving.
There are more famous algorithms for resource configuration optimization, such as genetic algorithm, ant colony algorithm, particle swarm optimization, and simulated annealing algorithm, but these algorithms have certain constraints when solving the bilevel programming model of manufacturing resource optimization configuration in this paper.The bilevel planning is a NP-hard (nondeterministic polynomial, NP) problem.The process of solving such problems is very complicated.The fast nondominated sorting genetic algorithm (NSGA-II) is an improvement of the NSGA algorithm and is one of the best evolutionary multiobjective optimization algorithms [26,27].
This paper uses the NSGA-II algorithm to solve the bilevel planning model under the discrete manufacturing environment.The general flow is shown in Figure 4.

Illustrative Example
6.2.Model Solving.NSGA-II algorithm was used to solve the model study (48).The initial population of the algorithm is 50; the crossover probability is 0.6; the mutation probability is 0.03; and the largest genetic algebra is 200.In the MATLAB 2016a calculation environment, the average fitness of each population under different fitness functions is calculated [29], as shown in Figure 6.In the 200-generation evolution process, the average fitness value after 50 generations tends to be stable.Therefore, after 50 generations of evolution, the Pareto optimal solution of the lower level optimization object of the manufacturing resource optimization configuration model is obtained.Calculate the Pareto frontier of optimal solution sets consisting of 50 solutions, as shown in Figure 7.Each point in the graph represents a Pareto optimal solution.The entire Pareto optimal solution set is located on the first-level Pareto frontier and distributed uniformly.That is, the ideal Pareto optimal solution set is obtained.The Pareto optimal solution of the lower-level optimization of the manufacturing unit optimization selection model is regarded as the feasible solution of the upper optimization goal, the corresponding target value of the upper optimization goal is calculated, and the global optimal solution of the optimal selection bilevel programming model is obtained, according to the advantages and disadvantages.The degree is ranked and the first five groups of solutions are listed.The group with the smallest target value is the global optimal solution that satisfies the optimal choice between the customer and the manufacturer.
The Pareto optimal solution of the lower level optimization object of the manufacturing unit optimization model is the feasible solution of the upper level optimization object; then, the corresponding the upper level optimization object values are calculated.Further, we get the global optimal solution of the bilevel programming model and sort them according to their pros and cons and take the first five solutions listed, as shown in Table 7.The group with the smallest object value is the global optimal solution of the optimization problem.The upper-level optimization objective of the bilevel programming model, presents the objective of the traditional manufacturing resource configuration model.The exhaustive method is used to calculate the functional target decisionmaking scheme, and all the obtained data are ranked according to the degree of pros and cons.The top 10 groups are listed, as shown in Table 8.
From Tables 7 and 8, the five optimal solutions (groups 1, 3, 4, 7, and 9 in Table 7) of the model study in this paper all contain the optimal top 10 optimal solutions.Therefore, the results of the model calculation in this paper not only satisfy the upper level optimization object but also satisfy the lower level optimization object of the model.

Discussion
In this paper, when establishing a bilevel programming of the discrete manufacturing resource optimization configuration model, the upper level optimization object is to select the smallest comprehensive value of the MRD and the lower level optimization object is to consider the maximum demand of the MRP.The upper level optimization goal is to select the upper level to consider the resource suppliers' comprehensive foreground value to be the smallest and the lower layer to consider the resource supply-side demand maximization, while the upper and lower layer constraints are.There are no interactions between the constraints imposed by the resource suppliers and the resource customers.The upper level and lower level give their own constraints, but the constraints are given by the MRD and the MRP, respectively, so there is no interaction between the constraints.We add a constraint condition which is strongly related to the upper level optimization object to the lower-level.A new bilevel programming model is established as follows: L max z = E, F, R T 50 s t V ij ≤ V max , 51 In the lower constraint, add (51) which represent the foreground value constraint of the MRD for each manufacturing unit.By adding (51) to the lower constraint conditions, when the MDP is optimizing the manufacturing resource configuration, it is a priority to discharge the prospect value of the MRD, thus simplifying the complexity of computing, speeding up the calculation speed, and further improving the interactivity between the upper level and the lower level.

Conclusions
The traditional discrete manufacturing object decisionmaking method only considers the requirements of the MRD and ignores the practical difficulties of the MRP or only from the perspective of the MRP to provide self-perceived the MRD satisfaction and then through the linear weighting method to convert the multiobjective optimization problem into a single-objective optimization problem to solve.This paper analyzes the importance and current deficiencies of manufacturing resource optimization configuration in a discrete manufacturing environment and considers uncertainties such as manufacturing resource changes and manufacturing task changes.On this basis, according to the process and characteristics of manufacturing resource configuration in a discrete manufacturing environment, and from the interests of different participants, a bilevel programming mathematical model for manufacturing resource optimization configuration was built.It not only ensures the interests of different participants but also ensures the smooth progress of the manufacturing service.Finally, it uses NSGA-II to solve the model.The method presented in this paper has the characteristics of clear concept and simple calculation process and has strong operability and practicality.It provides a new way to solve the problem of multiobjective resource optimization configuration.16 Complexity and b = b 1 , b 2 , b 3 , b 4 .When a 1 ≥ b 4 , it can be said that a is greater than b and recorded as a > b; when a 1 = b 1 , a 2 = b 2 , a 3 = b 3 , and a 4 = b 4 , it can be said that a is equal to b and recorded as a = b; when b 1 ≥ a 4 , it can be said that a is less than b and recorded as a < b.

4
and b = b 1 , b 2 , b 3 , b 4 .The distance between the trapezoidal fuzzy numbers a = a 1 , a 2 , a 3 , a 4 and b = b 1 , b 2 , b 3 , b 4 is defined as d a, b . d

Figure 2 :
Figure 2: Schematic diagram of the bilevel programming.
a cost set of the ith manufacturing subtasks, where C ij denotes the cost of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .(v)T i = T i1 , T i2 , … , T im i is a time set of the ith manufacturing subtasks, where T ij denotes the time of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .
a quality set of the ith manufacturing subtasks, where Q ij denotes quality of the jth manufacturing units of the ith manufacturing subtasks, i = 1, … , n and j = 1, … , m i .

6. 1 .
Model Establishment of Discrete Resource Configuration.The MDR submits the manufacturing task to the MDP according to the market demand.On the basis of certain

Table 2 :
Set of candidate manufacturing units.
) denote the

Table 5 :
The prospect value of MRD (one MRD's example).

Table 6 :
Evaluation information of MRP.