Next Article in Journal
Orthogonal Optimal Design of Multiple Parameters of a Magnetically Controlled Capsule Robot
Next Article in Special Issue
Research on Task-Service Network Node Matching Method Based on Multi-Objective Optimization Model in Dynamic Hyper-Network Environment
Previous Article in Journal
Diffusivity Measurement by Single-Molecule Recycling in a Capillary Microchannel
Previous Article in Special Issue
A Layer-Arranged Meshless Method for the Simulation of Additive Manufacturing with Irregular Shapes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Cloud 3D Printing Order Task Execution Based on the AHP-TOPSIS Optimal Set Algorithm and the Baldwin Effect

1
School of Electrical and Electronic Engineering, Shandong University of Technology, Zibo 255000, China
2
School of Mechanical & Vehicle Engineering, Linyi University, Linyi 276000, China
3
Shandong Longli Electronic Co., Ltd., Linyi 276000, China
4
School of Mechanical and Electrical Engineering, Wuhan City Polytechnic, Wuhan 430064, China
5
College of Mechanical Engineering, Anyang Institute of Technology, Anyang 455000, China
*
Authors to whom correspondence should be addressed.
Micromachines 2021, 12(7), 801; https://doi.org/10.3390/mi12070801
Submission received: 19 March 2021 / Revised: 7 June 2021 / Accepted: 28 June 2021 / Published: 6 July 2021
(This article belongs to the Special Issue 3D Printed Implants for Biomedical Applications)

Abstract

:
Focusing on service control factors, rapid changes in manufacturing environments, the difficulty of resource allocation evaluation, resource optimization for 3D printing services (3DPSs) in cloud manufacturing environments, and so on, an indicator evaluation framework is proposed for the cloud 3D printing (C3DP) order task execution process based on a Pareto optimal set algorithm that is optimized and evaluated for remotely distributed 3D printing equipment resources. Combined with the multi-objective method of data normalization, an optimization model for C3DP order execution based on the Pareto optimal set algorithm is constructed with these agents’ dynamic autonomy and distributed processing. This model can perform functions such as automatic matching and optimization of candidate services, and it is dynamic and reliable in the C3DP order task execution process based on the Pareto optimal set algorithm. Finally, a case study is designed to test the applicability and effectiveness of the C3DP order task execution process based on the analytic hierarchy process and technique for order of preference by similarity to ideal solution (AHP-TOPSIS) optimal set algorithm and the Baldwin effect.

1. Introduction

Intelligent algorithms are the most commonly used tool to solve NP-complete combination optimization problems. After years of development, many different random search strategies have emerged [1]. They all form their own iterative search mechanisms by simulating the behavior and evolution of natural ecology. They are simple, universal, robust, and easy to use in parallel processing [2]. To solve the problem of computing resource allocation in a cloud 3D printing service (C3DPS) resource pool, it is necessary to consider both the constraints of the task graph and the nodes. In the neighborhood search algorithm, the typical simulated annealing and tabu search algorithms are used due to the strong randomness and the need for only a single iteration based on an individual search. There is a very low probability of finding an optimal solution within the feasible solution set during a short iteration period. Presently, this algorithm is not suitable for solving combination optimization problems such as scheduling when combined with another algorithm application [3]. However, only the heredity and immunity algorithms show good performance in various applications and other improvements in evolutionary learning algorithms. Others, such as artificial neural network and DNA calculation algorithms, are less mature in the computational process of solving problems in mechanics as applied to combination optimization. In addition, a hybrid of ant colony and particle swarm optimization (PSO) algorithms was proposed to solve the multi-objective flexible scheduling problems based on the analysis of objectives and their relationship. This self-learning algorithm describes the application of the particle swarm optimization (PSO) metaheuristic to the continuous version of the p-median problem. To a certain extent, most of the improved algorithms convert the original numerical changes into individuals in a particle swarm, and the optimum crosscurrent operation or sub-optimal exchange has no obvious practical application. The optimum crosscurrent operation is also called genetic PSO [4].
The evaluation of C3DPS quality belongs to the field of multiple-criteria decision-making (MCDM) [5]. A study of relevant data around the world reveals that research on service quality is primarily based on scientificity, reliability, a combination of quantitative and qualitative aspects, and other principles. Therefore, there are three primary methods for comprehensively evaluating C3DPS quality: the classical comprehensive evaluation method of cloud services (CSs), the CS evaluation method based on multi-attribute decision-making, and the quality of service (QoS) [6,7]. These methods are used in the analytic hierarchy process (AHP), fuzzy mathematics, the effectiveness function, principal component analysis (PCA), and so on [8]. Among these methods, analytic hierarchy process (AHP) is a combination of qualitative and quantitative decision analysis methods, and it has a wide range of practicality for decision analysis of various types of problems. For the AHP decision-making problem of multi-expert evaluation, the judgments given by experts may be very different due to the differences in cognition of different experts, the differences in preferences, and the ability of experts. Fuzzy mathematics is an extension or a super-set of the Boolean logic, aimed at maintaining the concept of the partial truth. The utility function method is usually employed to represent the relationship between the utility obtained by the consumer in consumption and the quantity of the product portfolio [9]. The PCA is a statistical method. The main concept behind the PCA is to consider the correlation among features. If the correlation is very high among a subset of the features, PCA will attempt to combine the highly correlated features and represent this data with a smaller number of linearly uncorrelated features. The value of the PDF at any given point in the set of possible values taken by the random variable can be interpreted as providing the relative likelihood that the value of a random variable will equal to the sample. The TOPSIS method sorts potential solutions according to the closeness among a limited number of negative ideal solutions and the preference by similarity to ideal solution (PSIS), which is a comparison of the relative advantages and disadvantages of existing objects [10]. The DEA method is a quantitative analysis method that can input and output several indicators. It measures the linear programming efficiency of the amendment of a decision-making unit, by using data envelopment analysis (DEA) [11].
This study is conducted with the aims to develop a new ranking method for multicriteria decision-making problems with conflicting criteria. These methods include the AHP method, fuzzy mathematics, the effectiveness function, the PCA method, and the PDF method [12]. It performs a literature review of common multicriteria decision making methods, examines the advantages and disadvantages of the identified methods, and explains how their common applications relate to their relative strengths and weaknesses. The TOPSIS method is suitable for small restrictions on the number and types of evaluation indicators. It can be objectively optimized for several decision-making options [13]. The need for optimization of multiple objectives to reach a set of optimal solutions has recently been increased in many applications. Therefore, various types of multi-objective optimization algorithms have been developed for solving multi-objective problems. The Pareto optimal set algorithm has the advantages of simplicity, good convergence, and fast search speed. This algorithm retains the advantages and addresses the disadvantages of the improved local optimal solution or global optimal solution algorithm; its use in solving the resource allocation problem is relatively suitable and has potential.
The research on index systems is mainly divided into the reputation evaluation index, QoS evaluation index, service resource combination evaluation index, and other evaluation indexes [14]. According to the characteristics of the evaluation index, it is a cloud service evaluation method for multi-attribute decision-making. This method puts fewer restrictions on the characteristics of the evaluation index and the type of index and can optimize multiple decision-making schemes. An intelligent optimization algorithm has a simple algorithmic principle. The convergence and search speed are good, the evaluation efficiency is high, it is relatively easy to obtain satisfactory evaluation results, and the related research fields are expanding. Therefore, this paper proposes a multi-attribute decision service evaluation method based on an intelligent optimization algorithm, that is, a multi-attribute decision evaluation method based on a Pareto optimal set algorithm. To solve the problem of computing resource allocation, it is necessary to consider the constraints of both the task graph edges and nodes in the solution process of the C3DPS resource pool [15]. Among neighborhood search algorithms, typical simulated annealing and abut search have strong randomness and use only a single individual based on the iterative search, so the probability of finding a feasible solution is extremely low in a short time (iteration time). At present, they are mostly combined with other algorithms, but these are not suitable for solving optimization problems such as scheduling combinations.
Among evolutionary learning algorithms, only immune algorithms and genetic algorithms have been proven to show good performance in various applications and improvements. Other algorithms are less developed or more immature in the neighborhood of optimization, such as artificial neural networks and DNA computing. In addition, the ant colony algorithm and particle swarm algorithm are widely used in scheduling problems. However, the self-learning mechanism of PSO was originally designed for continuous numerical optimization problems [16]. For discrete combination optimization problems, the speed and position learning mechanisms of these algorithms are not suitable. Most of the improvements convert the original numerical change into a crossover operation or exchange sub-operation between an individual and the current optimal solution. To a certain extent, this method is transformed into a genetic particle swarm algorithm, and the effect is not obvious in practical applications. In summary, the intelligent optimization algorithm for a Pareto optimal set is based on the above-mentioned advantages and disadvantages [17,18,19]. It is a newly improved local optimal solution or global optimal solution algorithm that is relatively suitable and has great potential for solving the problem of computing resource allocation.

2. Framework for Cloud Service Evaluation Based on a Hybrid Multi-Objective BM-MOPSO Evaluation Model

Pareto optimality is a cornerstone concept in the field of optimization. Based on Pareto optimality theory and the AHP-TOPSIS evaluation model, this framework is objectively determined to the weight of each evaluation set and indicator of the AHP method MATLAB software. The TOPSIS method is analyzed to the closeness of the optimal solution and determined by the overall similarity between an optimal and an ideal solution [20]. Therefore, it is possible to construct a new AHP-TOPSIS evaluation ideal value approximation model for decision-making, which can be used in decision-making.
This cloud service evaluation is as follows: The task-order requirements are normalized by task decomposition and detailed to task-order requirements in this C3DPS platform. This task-order requirement uses a parsing function and an analysis task that forms a one-to-one atomic task [21]. According to the atomic task, search matching is used to perform supply-and-demand matching of service resources, form a dual feasible solution for the C3DPS resource candidate service set, and provide feedback to users. Therefore, Pareto optimal optimization algorithm is objectively determined to the weight of each evaluation of set and indicator and construct a comprehensive AHP-TOPSIS evaluation model based on a Pareto optimal set. Here, Pareto optimal optimization algorithm and many local search strategies perform within the multi-objective Evaluation method. This framework of the evaluation method of the C3DPS order task execution process based on the AHP-TOPSIS optimal set algorithm and Baldwin effect is shown in Figure 1.
Here, task-order requirement agents: each agent is registered and released to all the requirements as the condition attribute and forms a meta-model corresponding to the task-order demand; task-order decomposition agent: According to the task decomposition strategy of the C3DPD Platform, this meta-model O R m o d e l is expressed by the discrete attributes values of the decision table; Search-matching agent: this agent is a domain O R S = ( O R S 1 , O R S 2 , , O R S n ) of feasible solutions of the corresponding service resource candidate service subset [22]. Therefore, O R S i is a subset of candidate services, and these are matched by the atomic task O R i ; Service resource agent: according to the weight indicators of the domain O R S = ( O R S 1 , O R S 2 , , O R S n ) of feasible solutions, it corresponds to some subsets of candidate services by the Pareto optimal optimization algorithm, also supported to this C3DPS platform; Evaluation agent: the domain O R S = ( O R S 1 , O R S 2 , , O R S n ) of feasible solutions are solved that it corresponds to some subsets of this candidate services.

3. Intelligent Optimization Algorithm for the Pareto Optimal Set and AHP

3.1. Intelligent Optimization Algorithm for Pareto Optima

The Pareto optimal condition is as follows: In the process of information resource allocation, the marginal rate of technical substitution of two kinds of information resource X and Y in any computer resource is equal to the marginal rate of technical substitution of producing these two kinds of information resources X and Y; that is,
M R T X Y = M R S X Y A
where M R T is an ideal state of a computing resource allocation and M R S X Y A is a marginal rate of technical substitution of the two kinds of information resource X and Y in any computing resource.
The intelligent optimization algorithm for Pareto optima is derived from the above concept [23]. Based on the analysis of the existing intelligent optimization search strategies, a multi-objective optimization algorithm for Pareto optimal and AHP-TOPSIS evaluation models is proposed. It is a multi-objective evolutionary algorithm that combines evolutionary computation and a multi-local search strategy.
Definition 1.
Assuming that the outer set is in the i-th generation, each individual is assigned an intensity value, namely,
S i ( k ) = t i ( k ) n ( k ) + 1 i = 1 , 2 ... , n T ( k ) ,
where t i ( k ) = | { ( u | u i ( k ) < u )     s . t . u i ( k ) T ( k ) } | , T ( k ) = D o m S ( k ) T ( k 1 ) , n T ( k ) is the size of set T ( k ) , and t i ( k ) is the number of individuals who dominate the set.
A strength value is assigned such that each individual can be computed and multi-locally searched for the evolution of the Pareto optimal set intelligent optimization algorithm, and this value is the fitness determined by its advantages and disadvantages, namely,
f i ( k ) = 1 + u j T ( k ) u j u i S i ( k ) ,
Here, the fitness of an individual is equal to the sum of all the external individuals that dominate it. In addition, the size of the set is chosen so that the individual is infinitely close to the fitness value.

3.2. Analytic Hierarchy Process

In the 1970s, Professor T. L. Saaty, an American operational research scientist, proposed a multi-objective decision analysis method called the AHP, which is mainly applied to decision-making problems under multiple evaluation criteria [24].
The modelling steps of the AHP are as follows:
(1) Establish a hierarchical structure model.
The hierarchical structure model is decomposed into the various factors that are contained in the problem, which form several levels from top to bottom according to different attributes [25]. The structure model framework includes the following levels: the highest level (the overall goal of the complex system); the middle level (the planning and decision-making, the measures and adopted policies, and the criteria for achieving the goals), and the lowest level (various strategies and constraints).
Definition 2.
Assuming that the set of evaluation indicators is T = { t 1 , t 2 , , t n } , the evaluation indicators of each factor are compared with a pair of importance degrees, and the determined value α i j is set to represent the importance degree of t i . Then, the hierarchical structure model is as follows:
α i j > 0 , α i j = 1 α j i ,
where α i j is the scale.
Then, the judgment matrix D is as follows:
D = [ X 11 X 12 X 1 n X 21 X 22 X 2 n X m 1 X m 2 X m n ] = [ X 1 X 1 X 1 X 2 X 1 X n X 2 X 1 X 2 X 2 X 2 X n X n X 1 X n X 2 X n X n ] ,
It can be seen from the above that for the positive definite reciprocal judgment matrix D, the maximum characteristic root λ max exists and is unique, and the weight w i is composed of positive elements of vectors and is unique.
(2) Construct a judgment matrix.
The multi-level system is divided into several hierarchical levels according to different goals and functions [26]. Among them, the judgment matrix is constructed with the pairwise comparison method and comparison scale at the lowest level.
The elements in each row of the judgment matrix are multiplied as follows:
M i = Π i = 1 n α i j ,
M i is calculated by taking the n-th power root:
W i = M i n ,
W i in the vector is normalized:
W i = W i / i = 1 n W i ,
Here, the values of vector W i are normalized for the consistency check.
(3) Determine the hierarchical order and perform a consistency check.
Here, the method determines the correlation degree between adjacent level elements in the above judgment matrix. Through the construction of two comparison judgment matrices and the mathematical method of matrix operations, the importance order of the related elements is determined for a certain element in the previous level.
Definition 3.
In checking the consistency of the judgment matrix, knowledge and experience can meet different conditions. The specific formula is as follows:
C R = C i R i = λ max n ( n 1 ) R i ,
where C i is the indicator of the consistency check, n is the order of thejudgment matrix, and R i is the average value of the consistency check.
(4) Determine the hierarchical total ranking and perform a consistency check.
The combination weight vector of each layer element is calculated and ranked by the formula of the combination consistency check. Therefore, it determines the importance degree of each element at the bottom of the hierarchical structure model. In the traditional AHP method, it is difficult to verify the consistency of the judgment matrix, and this matrix is greatly affected by expert knowledge and preferences, which limits its promotion and application [27]. The concept of Multi-indicator Entropy is proposed, and the indicator weight solution method solves the problems of the traditional AHP method. The evaluation entropy refers to the confusion degree of the evaluation value of each candidate C3DPS set in a comprehensive evaluation.
Definition 4.
Assuming that for a certain service request, the set of candidate C3DPSs is C S = ( S 1 , S 2 , , S n ) and Q T 11 , Q T 12 , , Q T 1 n are defined as the evaluation values of the candidate C3DPSs on the indicators, the evaluation entropy of the indicators Q T 1 is
D i f f ( Q T 1 ) = i = 1 n ( Q T 1 i Q T 1 ¯ ) 2 n 1 Q T 1 ,
The solution formula for index evaluation entropy is introduced into other secondary indicators of the QoS index, and D i f f ( Q T 2 ) ,..., D i f f ( Q T i ) are obtained; then, the relative weight corresponding to Q T x     ( x = 1 , 2 , , n ) is
W ( Q T x ) = D i f f ( Q T x ) i = 1 n D i f f ( Q T x ) ,
where W ( Q T x ) is the weight value, D i f f ( Q T x ) is the indicator evaluation entropy of D i f f ( Q T x ) , and i = 1 n D i f f ( Q T x ) is the sum of the entropy values of Q T x     ( x = 1 , 2 , , n ) .
Based on Pareto optimal theory, the improved AHP-TOPSIS evaluation model is used to objectively determine the weight of each evaluation set and index [28]. The closeness of the evaluation model is determined by analyzing and making decisions to obtain the optimal solution that brings about an overall similarity between the best alternative and the ideal scheme combined with the TOPSIS optimization method.

4. Mathematical Model of C3DPS Order Task Execution Evaluation Based on the AHP-TOPSIS Evaluation Model

In the process of multi-attribute decision-making, the weight of each attribute reflects the relative importance of the attributes, which directly affects the result of decision-making. Therefore, the weight of each attribute is one of the key issues in multi-attribute decision-making that determines the results of decision-making. Aiming to solve a multi-attribute decision-making problem in which the attribute weight information is determined and the attribute value is an intuitionistic fuzzy number, a decision analysis method for the AHP-TOPSIS evaluation model is proposed. This analysis method is widely used in multi-attribute decision-making problems. The AHP-TOPSIS comprehensive evaluation model is a comprehensive evaluation and optimization method that combines the AHP and TOPSIS. The weight vector of each evaluation indicator is objectively determined and calculated to achieve comprehensive superiority by this evaluation model [29]. The basic principle of TOPSIS is to sort the evaluation objects by the distance between the fuzzy positive ideal solution and its similarity to fuzzy comprehensive attributes in the multi-objective decision-making problem.

4.1. Establish an Initial Evaluation Matrix A

Definition 5.
Suppose that a set A = { A 1 , A 2 , , A m } of schemes is composed of m schemes A 1 , A 2 , , A m and that each scheme also corresponds to a number of evaluation indicators X 1 , X 2 , , X n . The set of evaluation indicators is X = { X 1 , X 2 , , X n } . Then, the initial evaluation indicator matrix can be expressed as follows:
A = ( X i j ) m × n = [ X 11 X 12 X 1 n X 21 X 22 X 2 n X m 1 X m 2 X m n ] ,
where X i j is the j - th evaluation indicator in the j - th scheme.

4.2. Establish a Weighted Standardized Decision Matrix

The evaluation indicators can be divided into two categories: consumption indicators and profit indicators. The higher the profit indicator is, the lower the consumption indicator is. These are two kinds of indicators that can also be divided into measurement indicators and nonmeasurement indicators with different dimensions and dimensional units [30]. On this basis, it is necessary to strictly define the meaning of indicators and provide reference standards, which will eliminate the incontestability of the resulting indicators and carry out dimension normalization of the evaluation indicators in the evaluation of nonmetric indicators.
The standardized decision calculation method for the initial evaluation matrix is described below.
(1) The consumption indicator is
b i j = X i j min j ( X i j ) max j ( X i j ) min j ( X i j ) ,
where X i j min j ( X i j ) is the difference between the evaluation indicator value and the minimum evaluation indicator value in the initial evaluation matrix; max j ( X i j ) min j ( X i j ) is the difference between the maximum value and the minimum value in the initial evaluation matrix.
(2) The profit indicator is
b i j = max j ( X i j ) X i j max j ( X i j ) min j ( X i j ) ,
where max j ( X i j ) X i j is the difference between the maximum evaluation indicator value and the evaluation indicator value in the initial evaluation matrix, and max j ( X i j ) min j ( X i j ) is the difference between the maximum value and the minimum value in the initial evaluation matrix.
According to the life cycle evaluation indicator system of C3DPSs, a standardized decision matrix for multi-attribute decision-making is constructed. According to the accumulation process of the performance indicator M O pt ( T ,   Q ,   Mat ,   R ,   Rl ,   Flex ,   C ,   Ft ,   SF ,   Sa ) , the objective function of multi-objective optimization is determined, and the performance of the C3DPSs can be evaluated through four sub-evaluation indicator systems [31]. Assuming that the service resources of C3DPSs are evaluated, O R S = { r 1 , r 2 , , r n } will be selected as a candidate set of service resources in the C3DP order execution process. A standardized decision matrix for multi-attribute decision-making is shown in Table 1.
Here, ( μ i T , ν i T ) , ( μ i Q , ν i Q ) , ( μ i M a t , ν i M a t ) , ( μ i R , ν i R ) , ( μ i R l , ν i R l ) , ( μ i F l e x , ν i F l e x ) , ( μ i C , ν i C ) , ( μ i F t , ν i F t ) , ( μ i S F , ν i S F ) , and ( μ i S a , ν i S a ) are, respectively, intuitionistic fuzzy numbers representing the duration time T, cost C, matching degree Mat, service response R, quality Q, reliability Rl, service fault tolerance Ft, flexibility Fl, safety SF, and customer satisfaction Sa. For the performance evaluation of the C3DPSs life cycle evaluation indicator system, multi-attribute decision-making is standardized as a decision-making matrix.
(3) Establish a weighted standardized decision matrix R .
The weighted standardized decision matrix R is multiplied by the column vector of the matrix X with the weight values determined by the AHP method, and the weighted standardized decision matrix R can be obtained as
R = ( r i j ) m × n = [ w 1 b 11 w 2 b 12 w n b 1 n w 1 b 21 w 2 b 22 w n b 2 n w 1 b m 1 w 2 b m 2 w n b m n ] ,
Here, w is a coefficient of the matrix.
(4) Calculate the closeness B i + of the evaluation objects.
The ideal solution of the profit indicator set J 1 is the maximum value of the row vector, and the negative ideal solution is the minimum value of the row vector [32]. The ideal solution of the consumption indicator set J 2 is the opposite.
{ R + = { ( max ( w n b m n | m J 1 ) ) , ( min ( w n b m n | m J 2 ) ) } R = { ( min ( w n b m n | m J 1 ) ) , ( max ( w n b m n | m J 2 ) ) } .
The matrix of the evaluation object and ideal solution is as follows:
{ D + = j = 1 n ( r i j r j + ) 2 D = j = 1 n ( r i j r j ) 2 ,
where D + and D are the distances between the evaluation object and the positive and negative ideal solutions, respectively; r j + and r j are the elements corresponding to R + and R .
The formula for the closeness of the evaluation object is as follows:
B i + = D i D + i + D i ,     0 D + i 1 ,
(5) Construct an AHP-TOPSIS comprehensive evaluation model.
The evaluation matrix is constructed from the proximity analysis of the TOPSIS method, and the result vector F of the AHP-TOPSIS comprehensive evaluation is as follows:
F = W × B ,
In the formula, B is an evaluation matrix formed from the closeness value of each evaluation object, and W is the weight calculated by the analytic hierarchy process.

5. The Solution of the C3DPS Quality Evaluation Model

5.1. Hybrid Multi-Objective Particle Swarm Optimization (PSO) Algorithm Based on the Baldwin Effect (BM-MOPSO)

BM-MOPSO is an intelligent algorithm based on multi-objective PSO that combines Baldwin’s learning strategy idea, a population global target value scalar parameter, a scalar parameter for the population global objective value, intuitionistic fuzzy membership, and a ranking method. In the application of the hybrid multi-objective PSO algorithm, it solves the above-mentioned C3DPS quality evaluation problem [33]. This is the key to the problem of adopting the learning of the Baldwin effect within a certain period so that the global search and learning strategies can ensure the interactive operation of the algorithm in a fixed period.
(1) Baldwin effect learning strategy.
As a learning method, the Baldwin effect learning strategy can effectively reduce the selection pressure. This not only affects the characteristics of the search space but also increases the polymorphism of the genetic process and transforms the shape of the dominant search space. A local search can also be carried out based on the Baldwin effect that improves the nondominated solution of PSO [34]. The mathematical formula of this learning function is as follows:
Y i j = { p z j + k × ( p i j p z j ) ,     p z j < p i j p z j + k × ( p z j p i j ) ,     p i j < p z j p i j ,         O t h e r w i s e ,
where i , z = { 1 , 2 , , m } , m is the number of particles in the swarm, j = { 1 , 2 , , n } , and n is the dimension of the particle swarm.
If p i j is greater than p z j , particles p z j learn from particles p i j ; if p z j is greater than p i j , particles x i j learn from particles x z j ; if neither is dominated by the other, the particles do not learn.
(2) Local search strategy based on the BM-MOPSO algorithm.
When the particle swarm performs a local search for the Baldwin effect, one of the following situations will occur:
(1) When searching the initial points, most particles are far away from the Pareto frontier in the space, and it is easy to find the dominant solution in the region, which leads to high search efficiency. At a certain time, the optimal solution set of the population is sorted by Pareto dominance, and the particles are dynamically updated to the Pareto solution set. At the same time, the method learns from the dominant solution so that the particle swarm can more quickly approach the Pareto frontier.
(2) If the local search times of PSO converge to the threshold flag, the obtained solutions are all nondominated solutions, which indicates that the population has fallen into a local extreme point at this time. The direction of the nondominated solutions is
d = i = 1 flag x i x 0 | | x i x 0 | | ,
where x i ( i = 1 , 2 , , flag ) is the i th nondominated solution in the local search and x i x 0 is the distance from the local extreme point to the initial value in the local search.
(3) If the above two conditions are satisfied, then the number d of x i and the nondominated solution q (q = 3) satisfy one of the following conditions:
If d < q , the method learns from the dominating solution so that the particle swarm can more quickly approach the Pareto frontier.
If d q , the search direction changes, and the PSO can be guided and diffused in the optimization direction.
The local search method of the BM-MOPSO algorithm is as follows:
Step 1: The parameters are initialized; that is, the initial position of PSO is ( x 0 , y 0 ) , its iteration number is i = 0 , and the direction of the nondominated solution is d = 0 .
Step 2: In the initial position of PSO ( x 0 , y 0 ) , a dominating solution is randomly selected and marked as y 1 .
Step 3: If y 1 < y 0 , the initial position of PSO is locally searched by the Baldwin effect learning strategy and is set to x new = y 0 + k × ( y 1 y 0 ) ( k is the number of executions of the Baldwin effect learning strategy); then, the method goes to step 5. If y 0 < y 1 , then the Pareto solution set is dynamically updated. At the same time, the method learns from the dominant solution. The particle swarm can move closer to the Pareto frontier, ensuring that the individual particle swarm will have this learning ability. If x new = y 1 + k × ( y 0 y 1 ) , the method goes to step 5; otherwise, it goes to step 4.
Step 4: If d + + , then PSO calculates the direction d of the nondominated solution. If d q , it reinforces learning to calculate the optimal step size s and explore the optimal position x new = y 0 + d × s and then goes to step 5; otherwise, it goes to step 2.
Step 5: If i flag , the method goes to step 2; otherwise, this particle swarm population falls into a local extreme point, and is the next initial value of particle swarm x 0 .

5.2. The Basic Process of the Multi-Objective Particle Swarm Optimization Algorithm Based on the Baldwin Effect

The basic process of multi-objective BM-MOPSO is shown in Figure 2. The algorithm steps are analyzed below.
(1) Particle initialization.
First, the maximum number of iterations, number of independent variables of the objective function, maximum velocity of particles, and position information are set randomly in the velocity interval and search space to obtain a one-to-one mapping between the service resources and particles, that is, to perform particle initialization.
Assume that each particle is an optional combination of service resources, where a service resource is selected from each candidate service set list C R S i to form a combination of service resources [35]. The initial scale of the particles is the n-population of the feasible dimension space S = { p 1 , p 2 , , p n } , and each scale of particles corresponds to an optional combination of the number of service resources.
In the same way, the order of particles mapped one by one is. A mapping example of the composition of service resources and particles is shown in Figure 3.
For example, let the combination of C3DP optional service resources be { C R S 1 5 , C R S 2 4 , C R S 3 1 , C R S 4 5 } , which means that a C3DPS resource is selected by the serial atomic task sequence O R 1 to be the fifth service resource C R S 1 5 in the set of service resources C R S 1 . The serial atomic task sequence O R 2 is the fourth service resource C R S 2 4 in the set of service resources C R S 2 [36].
(2) Set of fitness functions.
From the above objective function M O pt ( OR ) , the optimal value of each index set is selected by the evaluation of the C3DP order task in the feasible solution domain of the service resource candidate service subset. Therefore, the moderation function of the BM-MOPSO algorithm is calculated as follows:
f i t n e s s ( p i ) = w 1 ( j = 1 n T ( C R S j p j i ) ) 1 + w 2 ( j = 1 n Q ( C R S j p j i ) ) 1 + w 3 i = 1 n ( j = 1 n M a t ( C R S j p j i ) ) 1 + w 4 i = 1 n ( j = 1 n R ( C R S j p j i ) ) 1 w 5 i = 1 n ( j = 1 n F l e x ( C R S j p j i ) ) 1 + w 6 ( j = 1 n C ( C R S j p j i ) ) 1 + w 7 i = 1 n ( j = 1 n F t ( C R S j p j i ) ) 1 + w 8 i = 1 n ( j = 1 n S a ( C R S j p j i ) ) 1
In the formula, w 1 , w 2 , w 3 , w 4 , w 5 , w 6 , w 7 and w 8 are the weight ratios of C3DP equipment service resources. The larger the fitness function f i t n e s s ( p i ) is, the better the particle p i is.
(3) The range of the particle dimension and moderate function.
In the BM-MOPSO algorithm, the learning strategy based on the Baldwin effect is an iterative process of learning and evolution that balances the relationship between global search and local search. To improve the nondominated solution of particle swarms, an individual particle learns in the same environment to achieve stronger survival adaptability [37].
After completing the Baldwin effect learning operation, the value of a certain dimension of the particle is larger than the value range, so it is necessary to calculate the particle swarm with the extreme value standardization method, which is used if the value is larger than the value range.
Therefore, the value range of each dimension of the particles is a discrete value range { p k j : 1 p k j K j } . After the Baldwin effect learning operation is completed, the dimension value of the particles is greater than the value range, and the extreme value normalization method is used for the particle swarm in this excessive value range. If p k j > K j , the value is set to p k j = K j .
The optimal position of each individual p i of the particle swarm is now x i = ( x i 1 , x i 2 , , x i j ) . The functional generalized derivative f ( x ) representing the feedback information in the learning and evolution process of the BM-MOPSO algorithm is calculated, which is defined as follows:
D f ( p i ) D x i k = f ( x i 1 , , x i k + x i k , , x i j ) f ( x i 1 , , x i k , , x i j ) x i k ,
Here, the individual particle p i yields p i when the learning strategy based on the Baldwin effect is carried out. For each dimension x i k ( k = 1 , 2 , , j ) of the position vector of individual, the formula is as follows:
min | p i p i | ε f ( p i ) = min | p i p i | ε f ( x i 1 , x i 2 , , x i j )                 = min N i 1 , N i 2 , N i j Z + f ( x i 1 + λ i 1 N i 1 D f ( p i ) D x i 1 , , x i j + λ i j N i j D f ( p i ) D x i j )
where λ i k and N i k are the parameters for adjusting the step length. The integer N i k is the number of individuals in the population particle swarm p i , and their initial values are all the same; N i k + 1 N i k and λ i k is the value of the updated particle position within the feasible region [ x L , x R ] [38]. Then, λ i is the quotient of the range distance and the maximum iteration number:
λ i = | x L x R | T max ,
The steps of the algorithm are as follows:
Step 1: Initialization. This includes all parameters of the particles, such as the initial position and velocity. The individual optimal position of a particle is defined as the current position, and the global optimal position is the optimal position of all particles. The initial position of each particle i is x i = ( x 1 , x 2 ) , the speed is v i G = ( v 1 , v 2 ) ( i = 1 , , N ), the number of particles is n p , the number of iterations is n max , and the initial solutions p 0 are randomly generated by n p .
Step 2: Calculate the fitness function value of each particle. That is, the fitness function value f i t n e s s ( p i ) is calculated by the functional generalized derivative f ( x ) when the coordinate x i takes the coordinates into the generalized derivative.
Step 3: Determine whether to perform local search. If the current iteration number meets the conditions for iteration termination ( < 1 3 of the total number of iterations), then local search is not performed and the algorithm moves to step 4; otherwise, step 5 is performed.
Step 4: Perform a global search. If the number of consecutive iterations with no updates meets the preset threshold, each particle will calculate the update speed and position in the global optimal solution according to the basic formula; otherwise, the algorithm moves to step 6;
Step 5: Perform a local search. The particles perform the local search based on the Baldwin effect with probability P k , and the algorithm moves to step 6;
Step 6: Non-uniform mutation. Assuming that there is a particle p t = ( v 1 , , v j ) in the t th generation, a random variable v k is selected and outputs a number between 0 and 1, and a non-uniform mutation operator is locally mutated by the particle for the next generation p t + 1 = ( v 1 , , v k , , v j ) . The particle swarm will become more stable.
Step 7: Repeat step 2 to step 6 until the current iteration number meets the conditions for termination; otherwise, continue to update x i and v i .
Step 8: After execution, output the calculation result.

6. Example Simulation

To ensure the preciseness of the data, the case data were sampled from a CMfg platform developed by a 3D printing company in Wuhan city, which is a smart C3DPS platform that integrates modelling design and 3D printing. It integrates various kinds of C3DPSs of multiple fields and types and provides network access to different types of 3D printing equipment. It also performs online real-time data collection. Here, the real data of different 3D printing devices in the platform were selected for the example simulation.

6.1. C3DPSs Modeling

Taking the construction of C3DPS network as an example, a computer program is clustered to the structure of the C3DPS network. The relevant basic data of the network must be described as a graph in the form of a database table, and the database table is the structure of 3D printing service network. The simulation process of the C3DPS network is divided into the following steps:
Firstly, it expresses the information of C3DPSs. As the assembly resources of personalized design products have many types and different capabilities, the C3DPS template is expressed in XML format.
For example, some of the attributes of Line 3 of No.6 Studio of Urui 3D Printing include 3D printing service mode, 3D printing specifications and dimensions, printing accuracy, maximum printing speed, 3D printing equipment model, and printing type. The node relationship is shown in Figure 4. As shown in Figure 5, it is a C3DPS expressed in XML format. At the same time, the XML format is expressed to the logical relationships among resources, attributes, and attribute values.
For example, some attributes of Line 3 of No.6 Studio of Yourui 3D Printing include 3D printing service mode, 3D printing specifications and dimensions, printing accuracy, maximum printing speed, 3D printing equipment model, and printing type, and their node relationships are shown in Figure 4. As shown in Figure 5 and Figure 6, it is a C3DPS in XML format. At the same time, the C3DPS also uses XML format to express the logical relationship among resources, attributes, and attribute values.
Secondly, it determines the service node. Here, as the research object, 10 C3DPS providers is selected and established to the basic information table of the node, as shown in Table 2.
Here, a program is incorporated in the software to cluster the structure of the C3DPS network, so as to describe the triple data of the cultural creative product personalized design product as a graph in the form of a database table and generate a complex network topology of the cultural creative product personalized design. Among them, the basic information of portrait 3D printing node and data attributes of C3DPS node are shown in Table 3.
Thirdly, it determines the alternative relationship between the candidate C3DPS providers (that is, there is a competitive relationship by the same or similar services providers) and creates a node connection relationship table.
Because of the different comprehensive service capabilities of each C3DPS provider, it is necessary to investigate each service node. For example, Jiayi Hi-Tech 3D scanner service provider provides 3D scanner outsourcing services, and the Worry 3D scanner service providers also provide 3D scanner outsourcing services. The output service capabilities of both are the same, so their similar weight is 1. After the data of all C3DPSs are summarized and sorted, the node connection relationship table is generated according to the form of a database table, as shown in Table 4.
Finally, the basic information table and connection relationship table of C3DPS node are imported into gephi-0.8.1 software for data analysis. The topology of C3DPS ( S _ N e t ) network is shown in Figure 7. It is composed of 33 nodes and 46 edges. Among them, the same color of nodes means that the types of C3DPSs are similar. The thickness of the edge represents the similarity between the two C3DPSs, that is, the service capability is replaceable.

6.2. Simulation Environment

To verify the reliability and universality of the cloud service evaluation based on the hybrid multi-objective BM-MOPSO evaluation model, a verification application in C3DPS creative and innovation product development was used. Its simulation environment was as follows:
  • Windows 7 operating system;
  • Intel (R) Core (tm) i5-4210H 2.90 GHZ CPU;
  • 8G memory.
The experimental environment was as follows:
At run time, the population size was 10, the maximum number of iterations was 100, the target weights were Q 1 and Q 2 , and the values were 0.7 and 0.3, respectively. The range of moderate function deterioration was set to −0.1 for the moderate functions w 1 = 0.20 , w 2 = 0.20 , w 3 = 0.10 , w 4 = 0.05 , w 5 = 0.10 , w 6 = 0.20 , w 7 = 0.05 , and w 8 = 0.10 .
According to the above multi-objective BM-MOPSO, the evaluation data and weight value of each candidate 3D printing device were inputted and simulated in MATLAB. Each algorithm was tested independently for each test function fewer than 30 times. The convergence characteristics of the algorithms are shown in Figure 8 and Figure 9:
Here, it is a flow chart of cloud service evaluation implementation based on the hybrid multi-objective BM-MOPSO evaluation model.

6.3. Analysis of Hybrid Multi-Objective BM-MOPSO

The Knowledge module is responsible for the management and scheduling of various types of knowledge on service resources; the Coordination module creates links between various coordination methods and performs operation monitoring and coordination management in the cloud service evaluation model based on the hybrid multi-objective BM-MOPSO evaluation model. Figure 10 shows the evaluation process for cloud services based on the hybrid multi-objective BM-MOPSO evaluation model.
A user proposes a complex manufacturing task that decomposes into six sub-task nodes T = ( T 1 , T 2 , T 3 , T 4 , T 5 , T 6 ) . After entering the task information into the platform, it is preliminarily matched to several 3D printing equipment resources that meet the user’s needs. A candidate set of 3D printing equipment resources is established, that is, S = ( S 1 , S 2 , S 3 , S 4 , S 5 , S 6 ) . Among them, task node S 1 is matched with three candidate resources, which can be expressed by S 1 = ( S 11 , S 12 , S 13 ) . Each candidate resource is shown in Table 5.
Where S = ( S 1 , S 2 , S 3 , S 4 , S 5 , S 6 ) means Services = (Service 1 (Point cloud data processing), Service 2 (3D Reverse Engineering), Service 3 (3D CAD modeling and design), Service 4 (Slicing procedures for layered manufacturing technique), Service 5 (simulation analysis), and Service 6 (3D printing and post-processing)}. S 1 = ( S 11 , S 12 , S 13 ) means three Candidates for 3D printing device resources.
According to the different 3D printing equipment resources, the evaluation indicators are optimized and quantified from the original data of each 3D printing equipment resource. The form of the data set is as follows:
S 1 = [ q c ( S i 1 ) q t ( S i 1 ) q r ( S i 1 ) q u ( S i 1 ) q e ( S i 1 ) q c ( S i 2 ) q t ( S i 2 ) q r ( S i 2 ) q u ( S i 2 ) q e ( S i 2 ) q c ( S i 3 ) q t ( S i 3 ) q r ( S i 3 ) q u ( S i 3 ) q e ( S i 3 ) q c ( S i j ) q t ( S i j ) q r ( S i j ) q u ( S i j ) q e ( S i j ) ]
Here, the unit of q c ( S i j ) is “Yuan”, and the unit of q t ( S i j ) is “hour”.
The above evaluation data of each candidate 3D printing equipment resource are specifically expressed as follows:
S 1 = [ 1240 35 0.76 0.79 78 1300 55 0.18 0.48 94 1750 54 0.44 0.64 98 ] S 2 = [ 980 42 0.79 0.75 20 2200 59 0.27 0.65 34 1700 61 0.65 0.16 70 3200 59 0.11 0.49 30 ] S 3 = [ 1400 49 0.95 0.34 53 1300 58 0.58 0.22 64 2700 67 0.75 0.25 87 ] S 4 = [ 1850 30 0.50 0.69 49 3300 34 0.89 0.95 40 2000 24 0.54 0.13 68 ] S 5 = [ 1600 55 0.14 0.25 49 2600 36 0.84 0.25 68 1850 37 0.81 0.24 75 900 44 0.92 0.35 53 ] S 6 = [ 750 43 0.19 0.25 35 2700 37 0.61 0.47 67 1800 38 0.35 0.83 59 1000 68 0.58 0.54 76 ] ,
Assuming that the particle population is 14 and the number of iterations is 50, the weight values of the evaluation parameters in the fitness function are, respectively, α 1 = 0.2 , α 2 = 0.3 , α 3 = 0.2 , α 4 = 0.1 , and α 5 = 0.2 .
According to the above multi-objective BM-MOPSO, the evaluation data and weight value of each candidate 3D printing device are inputted and simulated in MATLAB. The simulation results are shown in Figure 11, in which the vertical axis is the fitness function value and the horizontal axis is the population number.
It can be seen from Figure 6 that when the population number is 14, the fitness is lowest, the optimal combination scheme is S = ( S 11 , S 21 , S 31 , S 42 , S 54 , S 63 ) , and the fitness function value is 3.9370. That is, the optimal combination of 3D printing equipment resources is Yourui 3D printing hw-602, WINBO WB-SH105 and Qiaoyi workshop SLM 500, Flashcast Technology Explorer, Jiayi Hi-Tech JOYE-1212E, and Yunle Design Studio.

7. Conclusions

The global economy is transforming from a product economy to a service economy. Manufacturing and services are gradually merging. Individual enterprises pay close attention to their own core business. By providing manufacturing services, they can increase the value for stakeholders in the manufacturing value chain so that these individual enterprises will be more closely coordinated with each other. The C3DP model is a new service model that supports multi-user collaboration, and it is also an application of cloud manufacturing in the field of 3D printing. C3DPS modelling is the basis of C3DPS supply-and-demand matching; that is, this relationship between order tasks and services provides an effective way to match similar elements in the C3DPS platform. A large number of C3DPSs are aggregated according to certain rules and form a multi-complex 3D printing service network.
This paper formally describes C3DPSs, proposes a QoS acquisition and calculation method based on a mutual evaluation mechanism under the CMfg model, and establishes a C3DPS trust evaluation model based on service matching and global trust. The genetic algorithm optimizes the combination of C3DPSs so that it meets the requirements. Active intelligent rent-seeking for C3DPSs will be the next research direction.

Author Contributions

Conceptualization, C.Z. (Chenglei Zhang) and C.Z. (Cunshan Zhang); methodology, C.Z. (Chenglei Zhang); software, J.L. and J.Z.; validation, H.H. and B.Y.; formal analysis, C.Z. (Chenglei Zhang) and J.L.; investigation, C.Z. (Chenglei Zhang), J.L. and S.Z.; resources, R.L.; data curation, C.Z. (Chenglei Zhang); writing—original draft preparation, C.Z. (Chenglei Zhang), J.L. and K.Y.; writing—review and editing, C.Z. (Cunshan Zhang); visualization, J.Z. and H.H.; supervision, C.Z. (Cunshan Zhang) and R.L.; project administration, R.L.; funding acquisition, H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Natural Science Foundation of Shan dong Province of China under Grant (NO. ZR2019PEE019) and High-level talents (high-level doctorate) research project of Linyi University (NO. LYDX2019BS009) and Anyang medical rehabilitation project (NO. 2020-11).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the first author.

Acknowledgments

We sincerely thank Garrett Chapman Goble for his linguistic assistance during the preparation of this manuscript. We are grateful to the editors and anonymous reviewers for the valuable comments and suggestions on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ouelhadj, D.; Petrovic, S. A survey of dynamic scheduling in manufacturing systems. J. Sched. 2008, 12, 417–431. [Google Scholar] [CrossRef] [Green Version]
  2. Wong, T.; Leung, C.; Mak, K.; Fung, R. Dynamic shopfloor scheduling in multi-agent manufacturing systems. Expert Syst. Appl. 2006, 31, 486–494. [Google Scholar] [CrossRef]
  3. Xiang, W.; Lee, H. Ant colony intelligence in multi-agent dynamic manufacturing scheduling. Eng. Appl. Artif. Intell. 2008, 21, 73–85. [Google Scholar] [CrossRef]
  4. Rahman, I.U.; Zakarya, M.; Raza, M.; Khan, R. An n-state switching PSO algorithm for scalable optimization. Soft Comput. 2020, 11, 11297–11314. [Google Scholar] [CrossRef]
  5. Nie, L.; Gao, L.; Li, P.; Shao, X. Reactive scheduling in a job shop where jobs arrive over time. Comput. Ind. Eng. 2013, 66, 389–405. [Google Scholar] [CrossRef]
  6. Zeng, N.; Wang, Z.; Li, Y.; Du, M.; Liu, X. A Hybrid EKF and Switching PSO Algorithm for Joint State and Parameter Estimation of Lateral Flow Immunoassay Models. Trans. Comput. Biol. Bioinform. 2012, 9, 321–329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Cowling, P.; Johansson, M. Using real time information for effective dynamic scheduling. Eur. J. Oper. Res. 2002, 139, 230–244. [Google Scholar] [CrossRef]
  8. Hwang, H.C.; Choi, B.K. Workflow-based dynamic scheduling of job shop operations. Int. J. Comput. Integr. Manuf. 2007, 20, 557–566. [Google Scholar] [CrossRef]
  9. Kim, S.; Lee, J.; Choi, B. 3D printed fluidic valves for remote operation via external magnetic field. Int. J. Precis. Eng. Manuf. 2016, 17, 937–942. [Google Scholar] [CrossRef]
  10. Guerrero, G.; Langa, J.A.; Suárez, A. Architecture of attractor determines dynamics on mutualistic complex networks. Nonlinear Anal. Real World Appl. 2017, 34, 17–40. [Google Scholar] [CrossRef]
  11. Ikram, M.; Sroufe, R.; Zhang, Q. Prioritizing and overcoming barriers to integrated management system (IMS) implementation using AHP and G-TOPSIS. J. Clean. Prod. 2020, 254, 120121. [Google Scholar] [CrossRef]
  12. Fatima, A.; Cyril, G.; Vincent, V.; Stéphane, J.; Olivier, P. Towards normalization selection of Raman data in the context of protein glycation: Application of validity indices to PCA processed spectra. Analyst 2020, 145, 2945–2957. [Google Scholar] [CrossRef]
  13. Cheng, Y.; Tao, F.; Zhang, L.; Zhao, D. Dynamic Supply-Demand Matching for Manufacturing Resource Services in Ser-vice-Oriented Manufacturing Systems: A Hypernetwork-Based Solution Framework. In Proceedings of the ASME 2015 International Manufacturing Science and Engineering Conference, Charlotte, NC, USA, 8–12 June 2015; V002T04A017-1. p. 7. [Google Scholar] [CrossRef]
  14. Aslam, M.; Fahmi, A.; Almahdi, F.A.A.; Yaqoob, N. Extension of TOPSIS method for group decision-making under triangular linguistic neutrosophic cubic sets. Soft Comput. 2021, 25, 3359–3376. [Google Scholar] [CrossRef]
  15. Tao, F.; Cheng, J.; Cheng, Y.; Gu, S.; Zheng, T.; Yang, H. SDMSim: A manufacturing service supply–demand matching simulator under cloud environment. Robot. Comput. Manuf. 2017, 45, 34–46. [Google Scholar] [CrossRef]
  16. Xia, B.; Liu, T.; Ding, T.; Wang, Z. An Improved PSO Localization Algorithm for UWB Sensor Networks. Wirel. Pers. Commun. 2021, 117, 2207–2223. [Google Scholar] [CrossRef]
  17. Cheng, Y.; Tao, F.; Xu, L.; Zhao, D. Advanced manufacturing systems: Supply–demand matching of manufacturing resource based on complex networks and Internet of Things. Enterp. Inf. Syst. 2016, 12, 780–797. [Google Scholar] [CrossRef]
  18. Cheng, Y.; Tao, F.; Zhao, D.; Zhang, L. Modeling of manufacturing service supply-demand matching hypernetwork in ser-vice-oriented manufacturing systems. Robot. Comput. Integr. Manuf. 2016, 45, 59–72. [Google Scholar] [CrossRef]
  19. Raman, M.G.; Somu, N.; Kirthivasan, K.; Sriram, V.S. A Hypergraph and Arithmetic Residue-based Probabilistic Neural Network for classification in Intrusion Detection Systems. Neural Netw. 2017, 92, 89–97. [Google Scholar] [CrossRef] [PubMed]
  20. Sheng, B.Y.; Zhao, F.Y.; Zhang, C.L.; Yin, Y.; Shu, Y. 3D Rubik’s Cube—online 3D modeling system based on Web GL. In Proceedings of the Technology, Networking, Electronic & Automation Control Conference, Chengdu, China, 12–15 December 2017. [Google Scholar]
  21. Kim, J.-H.; Song, H.-Y. Hypergraph-Based Binary Locally Repairable Codes with Availability. IEEE Commun. Lett. 2017, 21, 2332–2335. [Google Scholar] [CrossRef]
  22. Zhang, C.L.; Sheng, B.Y.; Yin, X.Y.; Zhao, F.Y.; Shu, Y. Research and development of off-line services for the 3D automatic printing machine based on cloud manufacturing. J. Ambient. Intell. Humaniz. Comput. 2017, 10, 1109–1128. [Google Scholar] [CrossRef]
  23. Navabi, M.; Davoodi, A.; Reyhanoglu, M. Optimum fuzzy sliding mode control of fuel sloshing in a spacecraft using PSO algorithm. Acta Astronaut. 2020, 167, 331–342. [Google Scholar] [CrossRef]
  24. Minguella, J.; Villegas, M.; Poll, B.; Tena, G.; Calero, J.; Ginebra, M.-P.; Korkusuz, F. Automatic Casting of Advanced Technical Ceramic Parts via Open Source High Resolution 3D Printing Machines. Key Eng. Mater. 2015, 631, 269–274. [Google Scholar] [CrossRef]
  25. Im, S.; Lee, Y.; Kim, J.; Chang, M. A solution for camera occlusion using a repaired pattern from a projector. Int. J. Precis. Eng. Manuf. 2016, 17, 1443–1450. [Google Scholar] [CrossRef]
  26. Yasser, Y. A new knowledge-based link recommendation approach using a non-parametric multilayer model of dynamic complex networks. Knowl. Based Syst. 2018, 143, 81–92. [Google Scholar]
  27. Wang, X.; Sheng, B.; Zhang, C.; Xiao, Z.; Wang, H.; Zhao, F. An effective application of 3D cloud printing service quality evaluation in BM-MOPSO. Concurr. Comput. Pr. Exp. 2018, 30, e4977. [Google Scholar] [CrossRef]
  28. Yasser, Y.; Farshad, S. A statistical infinite feature cascade-based approach to anomaly detection for dynamic social networks. Knowl. Based Syst. 2017, 100, 52–64. [Google Scholar]
  29. Tao, F.; Zhang, L.; Lu, K.; Zhao, D. Study on manufacturing grid resource service optimal-selection and composition framework. Enterp. Inf. Syst. 2012, 6, 237–264. [Google Scholar] [CrossRef]
  30. Tao, F.; Li, C.; Liao, T.W.; Laili, Y.J. BGM-BLA: A new algorithm for dynamic migration of virtual machines in cloud computing. IEEE Trans. Serv. Comput. 2016, 99, 910–925. [Google Scholar] [CrossRef]
  31. Zhang, C.; Liu, J.; Xu, B.; Yuab, B.; Zhuang, S.; Zhao, F. Architecture of Cloud 3D Printing Task Modeling for Nodes Dynamic Scheduling and Coupling Based on Complex Networks. IEEE Access 2020, 8, 135208–135222. [Google Scholar] [CrossRef]
  32. Chen, Z.; Hendrix, W.; Samatova, N.F. Community-based anomaly detection in evolutionary networks. J. Intell. Inf. Syst. 2012, 39, 59–85. [Google Scholar] [CrossRef]
  33. Heard, N.A.; Weston, D.J.; Platanioti, K.; Hand, D. Bayesian anomaly detection methods for social networks. Ann. Appl. Stat. 2016, 4, 645–662. [Google Scholar] [CrossRef] [Green Version]
  34. Anwar, H.; Din, I.; Park, K. Projector calibration for 3D scanning using virtual target images. Int. J. Precis. Eng. Manuf. 2012, 13, 125–131. [Google Scholar] [CrossRef]
  35. Zhang, C.L.; Zhao, F.Y.; Wang, Z.Q. Modeling of Cloud 3D printing service hyper-network in service-oriented manufacturing systems. IEEE Access 2019, 8, 16225–16235. [Google Scholar] [CrossRef]
  36. Zhou, L.; Zhang, L.; Laili, Y.; Zhao, C.; Xiao, Y. Multi-task scheduling of distributed 3D printing services in cloud manufacturing. Int. J. Adv. Manuf. Technol. 2018, 96, 3003–3017. [Google Scholar] [CrossRef]
  37. Seo, Y.B.; Yun, Y.H.; Joo, K.-N. 3D multi-layered film thickness profile measurements based on photometric type imaging ellipsometry. Int. J. Precis. Eng. Manuf. 2016, 17, 989–993. [Google Scholar] [CrossRef]
  38. Liu, C.; Cao, G.-H.; Qu, Y.-Y.; Cheng, Y.-M. An improved PSO algorithm for time-optimal trajectory planning of Delta robot in intelligent packaging. Int. J. Adv. Manuf. Technol. 2020, 107, 1091–1099. [Google Scholar] [CrossRef]
Figure 1. Framework of the evaluation method of the C3DPS order task execution process based on the AHP-TOPSIS optimal set algorithm and Baldwin effect.
Figure 1. Framework of the evaluation method of the C3DPS order task execution process based on the AHP-TOPSIS optimal set algorithm and Baldwin effect.
Micromachines 12 00801 g001
Figure 2. The basic process of the multi-objective BM-MOPSO algorithm.
Figure 2. The basic process of the multi-objective BM-MOPSO algorithm.
Micromachines 12 00801 g002
Figure 3. Mapping example of the composition of service resources and particles.
Figure 3. Mapping example of the composition of service resources and particles.
Micromachines 12 00801 g003
Figure 4. The relationship diagram of C3DPS node.
Figure 4. The relationship diagram of C3DPS node.
Micromachines 12 00801 g004
Figure 5. The description diagram of C3DPS.
Figure 5. The description diagram of C3DPS.
Micromachines 12 00801 g005
Figure 6. The Partial XML information graph of C3DPS template.
Figure 6. The Partial XML information graph of C3DPS template.
Micromachines 12 00801 g006
Figure 7. The network graph of C3DPS Network ( S _ N e t ).
Figure 7. The network graph of C3DPS Network ( S _ N e t ).
Micromachines 12 00801 g007
Figure 8. Convergence characteristics.
Figure 8. Convergence characteristics.
Micromachines 12 00801 g008
Figure 9. Simulation results of the algorithm. (a) Iterative simulation of the objective function; (b) The final Pareto frontier based on the improved BM-MOPSO algorithm. A caption on a single line should be centered.
Figure 9. Simulation results of the algorithm. (a) Iterative simulation of the objective function; (b) The final Pareto frontier based on the improved BM-MOPSO algorithm. A caption on a single line should be centered.
Micromachines 12 00801 g009
Figure 10. Flow chart of cloud service evaluation implementation based on the hybrid multi-objective BM-MOPSO evaluation model.
Figure 10. Flow chart of cloud service evaluation implementation based on the hybrid multi-objective BM-MOPSO evaluation model.
Micromachines 12 00801 g010
Figure 11. Convergence curve of the fitness function.
Figure 11. Convergence curve of the fitness function.
Micromachines 12 00801 g011
Table 1. Standardized decision matrix for multi-attribute decision-making.
Table 1. Standardized decision matrix for multi-attribute decision-making.
r 1 r 2 r i r n
T ( μ 1 T , ν 1 T ) ( μ 2 T , ν 2 T ) ( μ i T , ν i T ) ( μ n T , ν n T )
Q ( μ 1 Q , ν 1 Q ) ( μ 2 Q , ν 2 Q ) ( μ i Q , ν i Q ) ( μ n Q , ν n Q )
Mat ( μ 1 M a t , ν 1 M a t ) ( μ 2 M a t , ν 2 M a t ) ( μ i M a t , ν i M a t ) ( μ n M a t , ν n M a t )
R ( μ 1 R , ν 1 R ) ( μ 2 R , ν 2 R ) ( μ i R , ν i R ) ( μ n R , ν n R )
Rl ( μ 1 R l , ν 1 R l ) ( μ 2 R l , ν 2 R l ) ( μ i R l , ν i R l ) ( μ n R l , ν n R l )
Flex ( μ 1 F l e x , ν 1 F l e x ) ( μ 2 F l e x , ν 2 F l e x ) ( μ i F l e x , ν i F l e x ) ( μ n F l e x , ν n F l e x )
C ( μ 1 C , ν 1 C ) ( μ 2 C , ν 2 C ) ( μ i C , ν i C ) ( μ n C , ν n C )
Table 2. Basic information table of service node.
Table 2. Basic information table of service node.
--Create table
create table C3DS_NODES {
ID INTERGER not null,
NAME VARCHAR(100), //Order Name
ORDERCATEGORY VARCHAR(50), //Order classification
SERVICECATEGORY VARCHAR(50), //Service type
PRINTMATERIAL NUMBER, //Printing material
PROCESSINGTECHNOLOGY VARCHAR(50), //Processing technology
STATUS VARCHAR(1), //Access status
REMARK VARCHAR(200)
}
Table 3. The attribute table of C3DPS node data.
Table 3. The attribute table of C3DPS node data.
RowKeyTimeStampColumns
OrnameOrcategoryPrtechnologyPrmaterial
000010Vatican gypsum relief1Gypsum 3D printing (PP)Gypsum
000020R2D2 robot3Selective laser sintering (SLS)Metal powder
000031Void cube model1Light curing (SLA)Photosensitive polymer
000040Eiffel Tower7Melt extrusion (FDM)Thermoplastic material
Table 4. The relation table of C3DPS node.
Table 4. The relation table of C3DPS node.
--Create table
create table C3DS_NODES
{
ID NUMBER not null,
SOURCE NUMBER, //Source node
TARGET NUMBER, //Target node
TYPE VARCHAR2(20), //Type
WEIGHT NUMBER //Weight
}
Table 5. Candidates for 3D printing device resources.
Table 5. Candidates for 3D printing device resources.
Candidate SetAtomic ServiceWorkshop NameEquipment Model
S 1 S 11 Yourui 3D printingHW-602
S 12 Jiayi Hi-TechJOYE-4035
S 13 Campus storeAurora LVO A8
S 2 S 21 WINBOWB-SH105
S 22 Beien 3DBANSOT M2
S 23 3D printing workshopDimension SST 1200es
S 24 The third brother of HanbangCorporate T1
S 3 S 31 Artful design workshopSLM 500
S 32 E-Plus-3DEP-M100T
S 33 ManhengEOS-M290
S 4 S 41 Tongchuang 3DMOONRAY
S 42 Flashcast Technology StudioExplorer
S 43 Yourui 3D printingDLP-1
S 5 S 51 Flashcast Technology StudioCreator Pro
S 52 Wuhan storeSecond-generation 3D printing
S 53 Jiayi Hi-TechJOYE-1010K
S 54 Jiayi Hi-TechJOYE-1212E
S 6 S 61 Campus Station of College of CultureFORTUS 200 mc
S 62 3D printing workshopProJet 6000
S 63 Yunle Design StudioULTRA
S 64 High-precision printingHOFTX2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, C.; Zhang, C.; Zhuang, J.; Han, H.; Yuan, B.; Liu, J.; Yang, K.; Zhuang, S.; Li, R. Evaluation of Cloud 3D Printing Order Task Execution Based on the AHP-TOPSIS Optimal Set Algorithm and the Baldwin Effect. Micromachines 2021, 12, 801. https://doi.org/10.3390/mi12070801

AMA Style

Zhang C, Zhang C, Zhuang J, Han H, Yuan B, Liu J, Yang K, Zhuang S, Li R. Evaluation of Cloud 3D Printing Order Task Execution Based on the AHP-TOPSIS Optimal Set Algorithm and the Baldwin Effect. Micromachines. 2021; 12(7):801. https://doi.org/10.3390/mi12070801

Chicago/Turabian Style

Zhang, Chenglei, Cunshan Zhang, Jiaojiao Zhuang, Hu Han, Bo Yuan, Jiajia Liu, Kang Yang, Shenle Zhuang, and Ronglan Li. 2021. "Evaluation of Cloud 3D Printing Order Task Execution Based on the AHP-TOPSIS Optimal Set Algorithm and the Baldwin Effect" Micromachines 12, no. 7: 801. https://doi.org/10.3390/mi12070801

APA Style

Zhang, C., Zhang, C., Zhuang, J., Han, H., Yuan, B., Liu, J., Yang, K., Zhuang, S., & Li, R. (2021). Evaluation of Cloud 3D Printing Order Task Execution Based on the AHP-TOPSIS Optimal Set Algorithm and the Baldwin Effect. Micromachines, 12(7), 801. https://doi.org/10.3390/mi12070801

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop