Abstract
To improve the satisfaction of both service demanders (SDs) and service providers (SPs) in the matching of cloud manufacturing (CMfg) tasks and services, a two-sided stable matching model of CMfg tasks and service considering the nonlinear relationship between satisfaction and expectations is proposed. As the expectations of SDs and SPs are difficult to be quantified directly, an evaluation method based on interval-valued hesitant fuzzy linguistic sets (IVHFLSs) is first presented. Next, a nonlinear model of satisfaction and expectations is built to quantify the satisfaction, which achieves accurate quantification of satisfaction. Then, a two-sided stable matching model of CMfg tasks and service is built, which takes the satisfaction of SDs and SPs as the optimization goals and considers the individual rationalities and blocking pairs. Finally, an adaptive genetic algorithm (AGA) is designed to solve the proposed two-sided matching model. A practical application and comparison analysis is used to verify the effectiveness and superiority of the research.
1. Introduction
Cloud manufacturing (CMfg) is a new manufacturing mode that organizes the online service released by service providers (SPs) and allocates them to service demanders (SDs) on demand [1, 2]. SDs expect to cooperate with SPs of low service cost and high service quality. SPs also expect to serve SDs with fast payment speed and high credibility [3, 4]. Hence, both SDs and SPs have expectations on the matching scheme of CMfg tasks and service. The comparison between the expectations of SDs for service and the actual situations of the service forms the satisfaction of SDs [5]. Similarly, the comparison between the expectations of SPs for tasks and the actual situations of the tasks forms the satisfaction of SPs. The higher the satisfaction of SDs and SPs, the stronger the competitiveness of CMfg platforms [6, 7].
Many studies concerning the matching of CMfg tasks and service have been performed. The research can be divided into two categories: single-sided matching and two-sided matching. Single-sided matching research has provided an important reference for two-sided matching research. For example, Lartigau et al. optimized the manufacturing cost, time, reliability, maintainability, and availability in research on CMfg service composition [8]. Zhang et al. developed an optimization model for service configuration in CMfg that took manufacturing cost, time, and quality as the optimization objectives [9]. Hsieh and Lin proposed a dynamic scheme of scheduling complex collaborative tasks for minimizing time and cost [10]. Xiong et al. maximized manufacturing efficiency and balanced service load during the service scheduling process [11]. Joglekar and Ford proposed a service allocation matrix to shorten the project duration [12]. Wang et al. optimized the cost, time, quality, and risk with a mathematical model [13]. Because single-sided matching research has mainly considered the expectations of SDs but ignored the expectations of SPs, and its effectiveness and feasibility need to be improved. To compensate, Zhao and Wang proposed a two-sided matching model of CMfg tasks and service [14]. First, SPs and SDs assessed each other based on linguistic information, and the assessment results were transformed into numerical values. Then, the satisfaction was calculated by using the variable fuzzy recognition method. Finally, a multiobjective optimization model was established, and the multiobjective functions were converted into a single objective function. Li et al. built a two-sided matching model with hesitant fuzzy preference information for configuring CMfg tasks and service [15]. First, SPs and SDs gave their expectations using the hesitant fuzzy element. Then, a two-sided matching model was constructed, and the optimal configuration results were obtained by solving the model. Ren and Ren put forward a one-to-many two-sided matching model, and the expectation utility theory was applied to calculate the satisfaction of SDs and SPs [16]. Zhao and Ding researched the two-sided resource matching mechanism and stability of a CMfg platform [17].
In summary, the two-sided matching research of CMfg service and tasks is in the exploratory stage at present. Relevant research has supposed that the relationship between satisfaction and expectations is linear. The more the matching scheme exceeds the expectations of SDs and SPs, the higher the satisfaction. They took the values of satisfaction indexes or the difference between the values and the expectations as the optimization objectives. However, satisfaction is a psychological feeling that forms after expectation is compared with the actual situation. The functions of expectations and satisfaction are nonlinear and segmented [5]. In some segments, the satisfaction is 0 or a fixed number, whereas in other segments, satisfaction increases or decreases between 0 and the fixed number.
When assessing satisfaction, some information such as expectations, the lower threshold, and the upper threshold needs to be determined in advance. The complexity of satisfaction indexes and the limitation of human cognition, hesitation, and fuzziness always exist in the expression process of information. Under such a situation, qualitative linguistic terms instead of precise quantitative numbers are more suitable for SDs and SPs to express information of satisfaction indexes [18]. Consequently, a number of linguistic decision-making methods emerged, such as the probabilistic linguistic term set (PLT) [18], hesitant fuzzy linguistic sets (HFLSs) [19, 20], interval-valued hesitant fuzzy linguistic sets (IVHFLSs) [21, 22], fuzzy preference relation with self-confidence (FPR-SC) [23], etc. These methods have been applied widely to various fields. In this article, the interval-valued hesitant fuzzy linguistic sets (IVHFLSs) are introduced, which are the extension of the linguistic term sets and interval-valued hesitant fuzzy sets (IVHFSs) [22]. Compared to PLT and HFLSs, IVHFLSs make the membership degrees no longer just relative to some crisp number, thus depicting the hesitancy of SDs and SPs, comprehensively. Compared to IVHFSs, IVHFLSs express the preference of SDs and SPs with linguistic terms, thus qualitatively representing the uncertain information. In a word, IVHFLSs are more convenient for SDs and SPs to express the information related to satisfaction and can effectively describe uncertainty, hesitancy, and inconsistency inherent in the decision-making process.
To bridge the gap, a two-sided matching model of CMfg tasks and service based on the quantification of satisfaction is proposed. Compared with the traditional matching model (as shown in Figure 1(a)), the main contributions of the article are as follows:(1)A nonlinear model of satisfaction and expectations is built (as shown in Figure 1(b)), which achieves accurate quantification of satisfaction and lays the foundation for the improvement of satisfaction in the matching of CMfg tasks and service.(2)A two-sided stable matching model of CMfg tasks and service is put forward (as shown in Figure 1(b)), which directly takes the satisfaction of SDs and SPs as the optimization goals and considers the individual rationalities and blocking pairs, thus making the optimization of satisfaction more effective.

(a)

(b)
The remaining sections of this article are organized as follows: Section 2 establishes a nonlinear model of satisfaction and expectations. Section 3 details a two-sided stable matching model of CMfg tasks and service on the basis of Section 2. Section 4 presents an example application and a comparison analysis of this article. Section 5 provides conclusions.
2. Quantification of Satisfaction Based on the Nonlinear Relationship
2.1. Expression of Expectations Based on IVHFLSs
To make the problem specific, suppose that SDs published I manufacturing tasks on a CMfg platform, and SPs issued J matching service. and indicate the CMfg task and the CMfg service, respectively.
Let be the m-th satisfaction index of SDs that is difficult to quantify directly, be a linguistic set, be the linguistic evaluation result of CMfg task Di for index , and be the membership interval of , then the expectation of the service demander published task Di on index based on IVHFLSs can be expressed as . Moreover, , where is the number of intervals. Similarly, if is a satisfaction index of SPs that is difficult to be quantified directly, the expectation of the service provider published Pj on index based on IVHFLSs can be expressed as .
2.2. Quantification of Expectations Based on Score Function
According to literature [21], the expectation function of can be denoted as follows:
Thus, the score function of can be represented as follows:
Similarly, the score function of can be represented as follows:
2.3. A Nonlinear Model of Satisfaction and Expectations
By fitting data, Liu et al. proposed that the piecewise exponential function is the subject function to the expression of customer satisfaction [5]. Especially, as the values of satisfaction indexes increase or decrease, the change in customer satisfaction can be divided into four areas: no-response area, defect area, overflow area, and saturation area. In the no-response area, the customer satisfaction is 0; in the defect area, the customer satisfaction increases or decreases rapidly and exponentially; in the overflow area, the customer satisfaction increases or decreases slowly and exponentially; and in the saturation area, satisfaction reaches its maximum and will not change. Again by fitting data, Zhou et al. found that the customer satisfaction has exponent relation to a satisfaction index [24]. Based on the above studies, a nonlinear model of satisfaction and expectations is built.
If is a positive correlation index, the satisfaction degree of the service demander published task Di with service Pj for index can be expressed as in the following equation (the corresponding satisfaction curve is shown in Figure 2(a)):

(a)

(b)
If is a negative correlation index, the satisfaction degree of the service demander published task Di with service Pj for index can be expressed as in the following equation (the corresponding satisfaction curve is shown in Figure 2(b)):
In equations (4) and (5), , , , and represent the expected values, the tolerable values, the upper threshold, and the lower threshold of task for index , respectively. means the actual values of service for index .
Similarly, if the satisfaction index of SPs is a positive correlation index, the satisfaction degree of the service provider published service Pj with task Di for index can be expressed as follows:
If is a negative correlation index, the satisfaction degree of the service provider published service Pj with task Di can be expressed as follows:
In equations (6) and (7), , , , and denote the expected values, the tolerant values, the upper threshold, and the lower threshold of service for index , respectively. means the actual values of task for index .
In equations (4)–(7), K is the overflowing coefficient of satisfaction. According to literature [5], the value of K can be set by experts, and its range is usually between 1 to 1.2. Based on equations (4)–(7), the satisfaction degree of Di with Pj and the satisfaction degree of Pj with Di can be expressed as in the following first and second equations, respectively:
In equations (8) and (9), and are the weights of satisfaction indexes and , respectively.
3. A Two-Sided Stable Matching Model of CMfg Tasks and Service
3.1. Model Construction
CMfg tasks can be divided into different types, such as machining, heat treatment, welding, and metal forming. Let q denote the q-th type of task and the 0-1 variables , , and denote the type of CMfg tasks, the type of cloud service, and the matching variable of tasks and service, respectively. If task Di belongs to the q-th task type, then ; otherwise, . If service Pj can operate the q-th task type, ; otherwise, . If task Di is matched with service Pj, then ; otherwise, .
Definition 1. The two-sided matching scheme between CMfg tasks and service is defined as mapping : . If , satisfies the following conditions: (1) ; (2) ; (3) , then is a matching pair of CMfg tasks and service.
If a service provider is incapable of the q-th type of task, the service provider would rather not provide the service to the q-th type of task. This phenomenon is called the individual rationality of SPs. Its specific definition is given below.
Definition 2. If service satisfies one of the following conditions: (1) ; (2) , , then it is called the individual rationality of SPs.
Definition 3. If all the SPs in the matching scheme are rational, then is an individual rational matching scheme of SPs.
If the SDs would rather not match service than match the unacceptable service, then this phenomenon is called the individual rationality of SDs. Its specific definition is given below.
Definition 4. If task satisfies one of the following conditions: (1) ; (2) , , and ; (3) (M is a large number), , then it is the individual rationality of SDs.
Definition 5. If all the SDs in matching scheme are rational, then is an individual rational matching scheme of SDs.
Definition 6. If is a matching scheme of both SPs’ individual rationality and SDs’ individual rationality, then is an individual rational matching scheme and satisfies the following equations:If there are blocking pairs in a matching scheme, the matching two sides may match privately, thus affecting the stability and effectiveness of the matching scheme [25]. The blocking pairs in a matching scheme of CMfg tasks and service are defined as follows.
Definition 7. For the matching schemes of CMfg tasks and service , if satisfy one of the following conditions: (1) , , and ; (2) , , , and , then is a blocking pair in the matching scheme .
Definition 8. If the matching scheme of CMfg tasks and service is individually rational and contains no blocking pairs, then is a stable matching scheme and satisfies the following equation: Based on the above analysis, the two-sided matching model of CMfg tasks and service can be established as in equation (13). The first and second optimization objectives are to maximize the satisfaction of SDs and satisfaction of SPs, respectively. The third optimization objective is to maximize the matching number of CMfg tasks and service. The meaning of constraints refers to Definitions 6 and 8,
3.2. Model Solution
To solve the proposed two-sided model, an adaptive genetic algorithm (AGA) was designed.
3.2.1. Coding
Real coding is adopted in this algorithm. Di = j (j = 0, 1, 2, …, J) denotes the task in which Di and service Pj form a matching pair.
3.2.2. Fitness Function
The objectives and in equation (13) are mutual restrictions. The optimization solution satisfying the two objectives is difficult to obtain, but the ideal optimal solution and the ideal worst solution of each objective are easy to get. To this end, the technique for order preference by similarity to an ideal solution (TOPSIS) is adopted to convert the optimization of and into the following equations:
In equation (14), represents the t-th objective function value of the r-th individual in population. and indicate the ideal optimal value and the ideal worst value of the t-th objective function, respectively.
3.2.3. Selection Operation
Given that tournament selection has better convergence than other selection operations [26, 27], the tournament selection strategy is adopted. First, select Ns individuals from population. Next, compare the fitness values between these selected individuals. Finally, put the individuals with the highest fitness values into a crossover pool. Cycle the above process until the crossover pool is full.
3.2.4. Adaptive Crossover
A two-point crossover is adopted here [28]. First, select individuals to be crossed as parents according to the crossover possibility pc. Second, randomly generate two integral points, rnd1 and rnd2, within the length of individuals as two crossover points. Finally, exchange the genes between the two integral points, rnd1 and rnd2, of two parents. Hereby, two children individuals are obtained. To avoid prematurity and local optimality, the adaptive crossover probability instead of a fixed one is adopted. The calculation method of the adaptive crossover probability is expressed as follows:
In equation (16), pc denotes an adaptive crossover probability, pcmin indicates the minimum crossover probability, pcmax denotes the maximum crossover probability, fmin indicates the minimum fitness values in contemporary populations, favg indicates the average fitness values in contemporary populations, and f ′ denotes the larger fitness values of the two individuals involved in the crossover.
3.2.5. Adaptive Mutation
Likewise, an adaptive mutation strategy is adopted here. The calculation method of the adaptive crossover probability is expressed as follows:
In equation (17), pm denotes an adaptive mutation probability, pmmin indicates the minimum mutation probability, pmmax means the maximum mutation probability, fmin denotes the minimum fitness values in contemporary populations, favg indicates the average fitness values in contemporary populations, and f ′ denotes the larger fitness values of the two mutated individuals.
4. Example Application and Comparison Analysis
4.1. Example Application
There are 6 metal forming tasks on a CMfg platform. After preliminary screening, a total of 9 services can operate these welding tasks. The sets of welding tasks and service are expressed as and , respectively. The satisfaction indexes of SDs include cost (z1), delivery time (z2), and credibility (z3); the satisfaction indexes of SPs consist of payment time () and credibility (). The delivery time represents the days required for processing and transportation; payment time is how many days SPs would get 95% of the full amount after a welding task is accepted (the other 5% is usually a warranty deposit).
4.1.1. Expression of Expectations Based on IVHFLSs
First, the SDs and SPs evaluate their expectations and self-conditions. For cost, delivery time, and payment time, which are easy to quantify directly, SDs and SPs can evaluate them with numbers. For credibility, which is difficult to quantify directly, SDs and SPs evaluate it with seven-level linguistic variables: L = {poorest, poorer, poor, general, good, better, best}. Then, the SDs and SPs use intervals to express the membership of linguistic levels. The evaluation results are shown in Tables 1 and 2, respectively. After the discussion of experts, the lower thresholds of cost, delivery time, and payment time are set to 200, 10, and 3, respectively. The overflowing coefficient of satisfaction K = 1.2.
4.1.2. Quantification of Expectations Based on Score Function
According to the linguistic scale function and equations (2)–(7), the satisfaction of SDs D1–D6 with service provider P1 can be obtained, as shown in Table 3.
Similarly, the satisfaction of SDs D2–D6 with the SPs P1–P9 and the satisfaction of SPs P1–P9 with SDs D1–D6 can be obtained. If the weights of cost, delivery time, and credibility are set to , respectively. The weights of payment time and credibility are set to , respectively. Then, the satisfaction matrices and can be obtained, as shown in Tables 4 and 5.
4.1.3. Obtaining the Soundest Alternatives
Based on the stable two-sided matching model shown in equation (13) and the proposed AGA, the soundest alternative was obtained.
The proposed AGA was programmed in Matlab 2014a. The initial population was set to 100, the maximum generation to 200, the minimum crossover probability to 0.7, the maximum crossover probability to 0.9, the minimum mutation probability to 0.08, and the maximum mutation probability to 0.1. After solution, the ideal optimal values of and were 6.517 and 6.670, respectively. The ideal worst values of and were 2.761 and 1.990, respectively. After solution, the best fitness values of every generation are shown in Figure 3.

When the algorithm converges, the optimal value of fitness function was 0.837, the optimal value of fitness function was 6, the corresponding optimal solution was , and the corresponding optimization objectives were .
4.2. Comparisons and Discussion
To further verify the superiority, the proposed two-sided matching model was compared with the traditional matching model proposed by Li et al. [15], which takes the values of satisfaction indexes as optimization objectives. Let represent the satisfaction expression of the service demander published task Di on service Pj for index based on IVHFLSs, indicate the satisfaction expression of the service provider published service Pj on task Di for index based on IVHFLSs, and mean the score of and , respectively. Hence, the satisfaction optimization goal of SDs based on the thought proposed by Li et al. [15] can be expressed as . In a similar way, the satisfaction optimization goal of SPs can be expressed as . The function of optimal matching numbers can be expressed as . Combining the constraints in Li et al. [15], the compared matching model is described as follows:
Then, the proposed AGA was programmed to solve the matching model proposed by Li et al. [15]. The parameters here are set same as the parameters used when solving the matching model proposed by this article. After solving, the optimal solution was . The corresponding optimization results were . Specifically, the satisfaction of SDs on SPs in the matching scheme is , , , , , and the satisfaction of SPs on SDs in the matching scheme is , , , , , .
As shown in Figure 4, the satisfaction of SDs optimized by this article is equal to the satisfaction of SDs optimized by Li et al. [15]. However, the satisfaction of SPs optimized by this article is increased 0.099 than the satisfaction of SPs optimized by Li et al. [15]. This indicates that the matching model proposed by this article is superior to the traditional model proposed by Li et al. [15].

Overall, satisfaction of SDs and SPs is an important property for improving the competitiveness of CMfg platforms. Although the two-sided matching problem of CMfg tasks and service has been studied in some articles, they took satisfaction indexes or expectations as the optimization goals of the matching model. Whereas, satisfaction is not linear with expectations. In fact, there is a limit to satisfaction, the bigger the satisfaction index, the higher the satisfaction is not necessarily in a range. In this article, the nonlinear model of satisfaction and expectations is built first. The satisfaction of SDs and SPs is taken as the optimization objectives of the two-sided matching model of CMfg tasks and service, thus making the optimization of satisfaction more effective.
Second, stability is necessary to consider when building the two-sided matching model, or else SDs and SPs may give up the current matching results [15]. In this article, stable matching conditions are considered when establishing the matching model, which not only can maximize satisfaction degrees of SDs and SPs but also can derive the stable matching results.
Third, IVHFLSs are introduced to SDs and SPs to express their expectations, tolerant values, and so on, which effectively describe the uncertainty, hesitancy, and inconsistency inherent in the decision-making process.
Therefore, the proposed matching model enriches the matching studies of CMfg tasks and service. However, the article also has some limitations. First, it cannot be used to deal with the conditions when some complex linguistic expressions are needed, such as multigranular unbalanced linguistic expression [29]. Second, it is assumed that the weight of satisfaction indexes is given by crisp numbers in advance. However, sometimes SDs and SPs may tend to provide criteria weight information with linguistic expressions, like “satisfaction index A is more important than in index B.”
5. Conclusions
The matching of CMfg tasks and service is an essential issue in the field of CMfg. Improving the satisfaction of SDs and SPs is significant to enhance the competitiveness of CMfg platforms. A nonlinear model of satisfaction and expectations and a stable two-sided matching model of CMfg tasks and service are proposed in this article. The proposed nonlinear model achieves the accurate quantification of satisfaction. The proposed two-sided matching model not only takes the satisfaction of SDs and SPs as the optimization objectives, but also considers the rationality of SDs and SPs and blocking pairs in the matching schemes, thereby improving the effectiveness.
In the future research, the following studies can be further conducted:(1)For the matching of CMfg tasks and service, due to different cultural and knowledge backgrounds of SDs and SPs, complex linguistic expressions like multigranular unbalanced linguistic expression may be needed. Therefore, it is necessary to develop new models with complex linguistic expressions.(2)In the proposed matching model, the weight of satisfaction indexes are given by SDs and SPs with crisp numbers. However, sometimes SDs and SPs may tend to provide criteria weight information with linguistic expressions in the real-world applications. Hence, how to deal with the weights of satisfaction indexes with linguistic expressions needs further study.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (No. 71971130).