Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 7343794, 12 pages

http://dx.doi.org/10.1155/2016/7343794

## A New Manufacturing Service Selection and Composition Method Using Improved Flower Pollination Algorithm

School of Information, Zhejiang University of Finance and Economics, Hangzhou 310018, China

Received 12 August 2016; Revised 16 November 2016; Accepted 22 November 2016

Academic Editor: Eusebio Valero

Copyright © 2016 Wenyu Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

With an increasing number of manufacturing services, the means by which to select and compose these manufacturing services have become a challenging problem. It can be regarded as a multiobjective optimization problem that involves a variety of conflicting quality of service (QoS) attributes. In this study, a multiobjective optimization model of manufacturing service composition is presented that is based on QoS and an environmental index. Next, the skyline operator is applied to reduce the solution space. And then a new method called improved Flower Pollination Algorithm (FPA) is proposed for solving the problem of manufacturing service selection and composition. The improved FPA enhances the performance of basic FPA by combining the latter with crossover and mutation operators of the Differential Evolution (DE) algorithm. Finally, a case study is conducted to compare the proposed method with other evolutionary algorithms, including the Genetic Algorithm, DE, basic FPA, and extended FPA. The experimental results reveal that the proposed method performs best at solving the problem of manufacturing service selection and composition.

#### 1. Introduction

With the rapid development of the Internet and associated information technology, many traditional manufacturing industries have changed the ways by which they provide products and services. Suppliers increasingly encapsulate their production capacities into services and upload them to shared network platforms. In general, user requirements are extremely complex and cannot be satisfied by a single manufacturing service. Therefore, to meet such increasingly complex requirements, appropriate candidate services must be chosen from a service set for each subtask and combined to form a composite service. The optimal composite manufacturing service is then identified from all the possible solutions and used to perform the multiobjective task [1]. As an increasing number of manufacturing services become available, the manner in which to select and compose them appropriately is a challenging problem. This is referred to as a multiobjective manufacturing grid resource service composition and optimal-selection problem [1].

As for the background of manufacturing service selection and composition, the concept of cloud manufacturing, as a new manufacturing paradigm, was proposed by Zhang et al. [2] and has attracted considerable attention from experts and scholars in the field of manufacturing service. Tao et al. [3] presented the classification of manufacturing resource and service and proposed a five-layered structure (i.e., resource layer, perception layer, network layer, service layer, and application layer) toward cloud manufacturing.

With the rapid increase in the number of web services having the same functional attributes, selecting an appropriate candidate service based on nonfunctional attributes such as quality of service (QoS) has become an equally challenging task that has also attracted considerable attention [4]. As with web services, each manufacturing service has functional and nonfunctional attributes [5]. Regarding nonfunctional attributes, QoS has become crucial in evaluating the overall performance of manufacturing service composition. Time, cost, reliability, and availability are four representative QoS properties that are frequently employed by researchers. Tao et al. [6] presented a QoS description supporting resource service correlation. They also described the basic resource service composite modes (RSCM) for composite resource service and the methods to translate a complicated RSCM into a simple sequence RSCM.

However, it is worth noting that, unlike typical web service composition, additional restrictions on manufacturing service composition exist [7]. For example, manufacturing enterprises belong to the entity enterprises that often produce considerable amount of contamination in the course of providing services. Therefore, from the perspective of green manufacturing [7], an environmental index (i.e., carbon emission) must be factored into the evaluation system to measure the overall performance of a manufacturing service composition. In addition, as the number of candidate services increases, the solution space becomes so large that finding the optimal solution is difficult. Thus, some researchers have proposed using the skyline operator [8] for service composition and have advanced the concept of skyline service [9, 10] to limit the size of the solution space and improve the combined utility of composite solutions.

Even after applying of the skyline operator, selecting an optimal composite solution remains time-consuming. Although many evolutionary algorithms, like Genetic Algorithm (GA), Particle Swarm Optimization (PSO) algorithm, and Differential Evolution (DE) algorithm, have been proposed to solve this problem, room remains for improvement. The Flower Pollination Algorithm (FPA) was first proposed by Yang [11] in 2012 as a new metaheuristic algorithm. It has been applied to many fields and has been shown to perform well, especially at solving optimization problems.

In this study, a novel and efficient method is proposed that employs not only the skyline operator but also the improved FPA (IFPA) to solve the problem of QoS-based manufacturing service selection and composition. The use of the skyline operator means that only those services that belong to the skyline service set are valid candidate services [9]. In this study, the differences between our proposed IFPA and basic FPA are the following three points. First, a grouping strategy and elite replacement strategy are used to improve the convergence speed of FPA. Second, we not only use the adaptive parameter to adjust dynamically the switch probability between global and local pollination but also employ an adaptive scaling factor to control the step size of the global search. Third, we introduce crossover and mutation operators of the DE algorithm into the local search of FPA to improve its performance.

The rest of this paper is organized as follows. Related studies about manufacturing service composition and the development of relevant evolutionary algorithms are discussed in Section 2. Section 3 provides a description of the model that is established in this study and definitions of skyline service. The basic FPA and its improvements are introduced in Section 4. Section 5 offers a case study and several comparative experiments used to illustrate the practicality and effectiveness of the proposed method. A conclusion that outlines the authors’ main contributions and that describes future research is provided in Section 6.

#### 2. Related Work

##### 2.1. Manufacturing Service Selection and Composition

With the rapid development of cloud computing and the advent of big data, the problem of optimizing composite services is causing increasing concern in the manufacturing industry. A composite service is the result of aggregating multiple services into one, and its aim is to execute functions of greater complexity than those that can be performed by a single service [12].

To evaluate the overall performance of service composition, Zeng et al. [13] set QoS attributes as the main measurement, including execution price, execution time, reputation, reliability, and availability, which were widely accepted and applied. To obtain high quality of service composition, they also introduced a middleware platform to express user satisfaction as utility functions. Tao et al. [14] presented different aggregation formulations for four different structures of service composition, including sequence, parallel, selective and cycle. In addition, the correlation of service composition was taken into consideration. To minimize implementation time and cost and to maximize reliability, an improved Particle Swarm Optimization (PSO) algorithm was introduced as the key technology. Considering the heterogeneous QoS values and homogeneous resources in manufacturing service, Zheng et al. [15] took the dynamic customer requirements into account by introducing the fuzzy theory and presented the way of calculating fuzzy QoS value and selecting the best services considering design preference. Hao et al. [16] pointed out that, in order to obtain the best service composition solution, user QoS constraints need to be relaxed rather than overlooked, so they proposed a novel QoS model and computational framework. Tao et al. [17] made an overview about manufacturing service management and analyzed the future research directions of manufacturing service composition and optimization.

However, most of the aforementioned studies only focus on the general QoS attributes and ignore the importance of environmental protection, which has become one of the most important attributes in the process of selecting manufacturing service [7, 17]. Therefore, an environmental index such as carbon emission is added in this paper as a soft constraint to measure the overall performance of manufacturing service composition.

As a multiobjective optimization problem, its solution space is very large and the objectives of manufacturing service composition are in conflict with each other. Thus, Alrifai et al. [18] first introduced the skyline operator [8] into service composition to reduce the number of candidate services. They also proposed the definition of skyline service to improve the efficiency of service selection. Skyline service has been widely used by scholars and become an active research area [10, 12]. In view of its good performance, the proposed method in this paper incorporates the skyline service concept in the model to limit the size of the solution space and improve the overall quality of the optimal manufacturing service composition solution.

##### 2.2. The Development of Evolutionary Algorithms in Manufacturing Service Composition

In order to solve the problem of service composition, experts and scholars have used a variety of evolutionary algorithms. For example, in the authors’ previous work [19], the extended GA was employed to solve the problem of manufacturing resource composition, treating the overall utility of multiple objectives as the metric of the composite solution. Li et al. [4] employed a multipopulation GA to solve the problem of QoS-oriented web service composition for the Internet of Things. Tao et al. [1] presented the concept of multiobjective manufacturing grid resource service composition and optimization problem and used the PSO algorithm to solve it. They not only introduced the basic resource service composite modes (RSCM) and the principles to translate a complicated RSCM into a simple sequence but also used the extended PSO algorithm to improve the quality of the optimal solution. Pop et al. [20] introduced a modified DE algorithm with a new genome-encoding method to solve the problem of service composition. Wei and Liu [21] introduced the Ant Colony Optimization (ACO) algorithm into the cloud manufacturing resource allocation model. Wang et al. [22] proposed a new optimization algorithm called the culture max-min ant system, which not only established a novel comprehensive model including both generic and domain QoS but also combined the Culture Algorithm (CA) and the ACO algorithm into a new approach. Lartigau et al. [7] applied an adapted Artificial Bee Colony (ABC) optimization algorithm to the cloud manufacturing service composition model, which not only included QoS attributes but also took the geoperspective transportation into consideration.

Although many conventional evolutionary algorithms have been used in service composition, each has its own deficiencies, for example, converging too easily to local optimum or premature convergence. In order to obtain better optimization solutions, a variety of emerging algorithms have been proposed. For example, inspired by the biological process of flower pollination, Yang [11] first proposed FPA in 2012 and proved that it was more efficient in solving optimization problems than both the GA and PSO.

As a novel and effective metaheuristic algorithm, FPA has been used in many fields. Abdelaziz and Ali [23] used it in a multimachine power system to seek the optimal parameters of static VAR compensator damping controller. Abdelaziz et al. [24] also used FPA to solve the problem of economic load dispatch and combined economic emission dispatch problems. To solve sizing optimization problem of truss structures, Bekdaş et al. [25] used FPA to minimize the weight of truss structures. To obtain feasible optimized designs, they added an iterative constraint handling strategy into the FPA search engine. El-henawy and Ismail [26] extended the basic FPA by combining it with chaotic search to solve large integer programming problems. Wang et al. [27] applied FPA with bee pollinator in the field of cluster analysis to make up for deficiencies in the -means method. By comparing it with other algorithms such as GA, they proved the better accuracy and higher level of stability of their proposed method. In order to solve the problem of multilevel image thresholding, Ouadfel and Taleb-Ahmed [28] applied FPA to search for the best values from a given number of thresholds with the purpose of maximizing the objective function. To apply the FPA to solve the antenna positioning problem, Dahi et al. [29] proposed four binary variants of the FPA by applying the principal mapping techniques. In addition, to study the applicability of this new metaheuristic algorithm in real-world applications, Draa [30] reviewed the performance of FPA in terms of both qualitative and quantitative aspects, and then he summarized the defects of the basic FPA by comparing it with state-of-the-art algorithms and pointed out the extension and application directions of FPA.

However, even though FPA research is expanding, to the best of our knowledge, it has seldom been used to solve the optimization problem of service composition, let alone be applied in the field of manufacturing service selection and composition. Therefore, in this study, FPA is applied to solve the manufacturing service selection and composition problem by improving the standard FPA in combination with DE. The experimental results show that the proposed IFPA can obtain better optimal solutions than other algorithms such as GA [4], DE [31], basic FPA [11], and even extended FPA (EFPA) [32].

#### 3. Model of Manufacturing Service Composition Problem

##### 3.1. Notations in Manufacturing Service Composition Model

To solve the manufacturing service composition problem and translate it into a mathematical model, we use the following notations to better illustrate the model of manufacturing service composition: : execution time of the candidate service : execution cost of the candidate service Ava: availability of the candidate service Rel: reliability of the candidate service : the number of subtasks in each basic structure : the selected manufacturing service realizing th subtask, : the total execution time in the sequence structure : the total execution cost in the sequence structure : the total availability in the sequence structure : the total reliability in the sequence structure : the processed carbon emission value of th manufacturing service : a soft constraint about carbon emission that is provided by users : real carbon emission value of th candidate service.

##### 3.2. The Calculation of QoS Properties

In this study, the four representative QoS properties of execution cost , execution time , reliability , and availability are employed as criteria for manufacturing service performance evaluation [13, 33].

As there are various metrics for different QoS attributes, in order to ensure the accuracy of assessment results, all the QoS attributes need to be normalized and converted into the same order of magnitude before they are employed in the final QoS aggregation formula. It has been well recognized that service quality is negatively correlated with cost and time, but it has a positive correlation with reliability and availability. Therefore, each QoS attribute is classified as being either positive attribute or negative attribute. The corresponding calculation formulas are as follows [5]:

Here, is th true attribute value of a manufacturing service and represents the corresponding normalized value. represents the maximal value of th QoS attribute for manufacturing services in the same candidate service set, that is, ; stands for the minimal value of th QoS attribute for manufacturing service, that is, .

In general, manufacturing service composition includes four basic structures: sequence, parallel, selective, and cycle structures. There are different aggregation methods for each structure. The aggregation formulas are defined in [13]. For example, the following are the calculation formulas for the sequence structure:

In order to use a unified standard to measure the overall performance of manufacturing service composition, the multiple objectives need to be transformed into a single one by assigning different weights to each attribute according to the different preferences of users. The final QoS evaluation value of manufacturing service composition can be calculated according to the following equation: where , , , and sum to 1 and represent the respective weights of the four attributes, including total time (), total cost (), total reliability (Rel), and total availability (Ava), of manufacturing service composition.

##### 3.3. The Calculation of Environmental Index

In addition to the basic QoS indicators, the proposed model contains an environmental indicator. Customers generally consider the QoS performance of manufacturing service composition to be more important than its environmental performance, and therefore they have different expectations of them. For that reason, the constraints proposed by customers can be divided into two types: rigid constraints and soft constraints. QoS constraints are defined as being rigid ones that must be satisfied, whereas environmental indices are defined as soft constraints that do not have to be satisfied strictly [22].

In the present paper, the authors regard carbon emission as an environmental protection index. As with manufacturing tasks, customers impose their own soft constraints that should be satisfied as far as possible. For each candidate manufacturing service , the deviation value between its carbon emission and soft constraint is calculated according to the following equation [22]:where is the deviation carbon emission value, is the real value (i.e., carbon emission value) provided by suppliers, and is a soft constraint about carbon emission that is provided by users. From the above equation, it is clear that is a binary index. If carbon emission meets the user requirement for environmental protection, then ; otherwise, . It is well known that the smaller the deviation value, the higher the environmental indicator. Thus, the deviation value is a negative attribute, and its standardized formula is defined in (2) too. The aggregation formulas of deviation value about carbon emission are similar to the aggregation formula of execution cost, so we omitted them here.

##### 3.4. The Comprehensive Utility of a Manufacturing Service Composition

Based on the above description, the final objective function of manufacturing service composition is defined aswhere is the comprehensive utility of a manufacturing service composition, represents the final QoS evaluation, and represents the final carbon emission value of manufacturing service composition. Terms and sum to 1 and represent the weights of QoS and environmental index in a manufacturing service composition solution, respectively. In this model, the goal is to obtain the maximum value of the comprehensive utility, and this comprehensive utility be regarded as the fitness value in the experiment section of this paper.

##### 3.5. The Introduction of Skyline Service

Before the new evolutionary algorithm is carried out to solve the problem of manufacturing service composition, the skyline service is introduced to reduce the number of candidate manufacturing services and improve the quality of optimal solution. The concept of skyline service, which was based on dominance relationship, was first proposed by Alrifai et al. [18]. The following is a detailed description of the related definitions [9].

*Definition 1 (dominance). *One assumes that there are two vectors and in -dimensional space. If the following two conditions are satisfied, one says that dominates . Firstly, each dimension in is greater than or equal to the corresponding in . Secondly, at least one dimension in is strictly greater than the corresponding in . The two conditions can be expressed mathematically as follows [34]:Otherwise, if either of the two conditions is not satisfied, then does not dominate .

*Definition 2 (skyline service). *For two manufacturing services and , if all the QoS attributes of are better than or equal to the corresponding ones of , and at least one attribute of performs strictly better than the corresponding attribute of , then dominates , which is denoted as [34]. That is to say, is better than . A skyline service refers to such a service that is not dominated by any other services.

The application of skyline operator can limit the size of the solution space of manufacturing service composition because the process of choosing skyline service eliminates those other services that could not possibly be selected as the final manufacturing service composition.

#### 4. Improved Flower Pollination Algorithm

##### 4.1. Introduction to the Standard FPA

FPA is a novel bionic evolutionary algorithm that was first proposed by Yang in 2012 [11]. The inspiration for FPA comes from the natural pollination process that takes place in flowering plants. In FPA, each flower stands for a feasible solution and the objective function value in Section 3.3 is regarded as the fitness value. To mimic the pollination process, there are two different pollination methods for each flower to choose: global or local pollination with switch probability . The specific process is described as follows.

For each individual, there is a random number rand. If , then global pollination should be carried out. In the global pollination process, each flower updates its position according to the following equation [11]: where and are the old and new positions of th flower, respectively, is the best flower at current iteration , which has the best fitness value in the whole population, and is the scaling factor that controls the step size of global pollination. Parameter , the Lévy flight, is used as the strength of pollination in the basic FPA, and the step size obeys the Lévy distribution: where is the standard gamma function and is from 0.3 to 1.99 [35]. In the basic FPA, Yang [11] set the distribution factor as 1.5. This distribution is valid when large step is bigger than 0.

If , the local pollination process should be carried out. Flower obtains its new position according to the difference between its old position and the position of two neighboring flowers and . This process is considered as local search, and the updating equation [11] is defined aswhere is drawn from a uniform distribution , and it is considered as a local random walk.

After pollination is completed, the new individuals update their positions by comparing fitness values. If the fitness of is better than that of , the new position of th flower will be replaced by . Otherwise, th flower remains at .

##### 4.2. Improvements of IFPA Compared with the Standard FPA

Like many other evolutionary algorithms, basic FPA suffers from premature convergence and can easily fall into local optimum. In order to overcome these shortcomings, the authors propose an IFPA to extend the standard FPA in three aspects. Firstly, the whole population is divided into different groups on average, and an elite replacement strategy is introduced to replace the worst individual in each group. Secondly, the switch probability and scaling parameter are adapted dynamically with the iterative process rather than remaining fixed like the basic FPA. This modification increases not only the global search ability but also the local search ability of basic FPA. Thirdly, in order to enhance its local search ability and obtain better optimal solution, the local search method is combined with the mutation and crossover operators of DE.

###### 4.2.1. Grouping Strategy and Elite Replacement Strategy

In the original FPA, flowers search for the optimal solution in a single population only. This method may take a relatively long time and its global search ability is insufficient. Based on this idea and inspired by [32], the authors introduce a parallel grouping strategy to solve this problem effectively. Before starting the iteration, the whole population is divided into three different groups on average according to their fitness values. The detailed operation steps are described as follows. Firstly, the fitness is calculated and ranked. Secondly, the whole population is divided into three groups on average according to their fitness. The first group consists of individuals with the best fitness, the third group consists of individuals with the worst fitness, and the remaining individuals belong to the second group. The iterative process is carried out in each group concurrently.

In addition, in order to improve the FPA convergence and the optimal solution quality, an elite replacement strategy is introduced. From each group, the best individual is selected as an elite individual. That is to say, after global and local pollination, the worst individual is replaced with the elite individual in each group. If there are duplicate individuals, then the duplicate one can be modified by selecting one dimension randomly to carry out mutation before starting the next iteration. In each group, FPA is carried out concurrently in the same way. Before executing the next generation, these groups are combined into a single population.

###### 4.2.2. Adaptive Switch Probability and Scaling Factor

In standard FPA, Yang [11] set the switch probability as 0.8. However, a fixed switch probability is insufficient. In order to retain a balance between global and local pollination in FPA, the switch probability is generated dynamically according to the fitness value of the current generation, which is represented by the objective function. The adaptive switch probability is obtained according to the following equation: where and are the maximum and minimum fitness values, respectively, of an individual in the current iteration.

For the same reason, in order to escape local optimum and avoid premature convergence in global pollination, the scaling factor is also generated dynamically according to the number of iterations. The calculation of scaling factor depends on the following formula: where represents the maximum iterations and is the current iterations. According to this formula, the value of scaling factor ranges from 0.2 to 0.5, and its value decreases with the increase of iterations.

###### 4.2.3. Modified Local Pollination Process by Combining with DE Algorithm

In the local pollination process of , the individuals and are selected randomly except for the rule that and should differ from the objective individual . However, this method is not sufficiently far-sighted and there is no direction of evolution to follow, which may lead the algorithm to fall into local convergence. In order to avoid this, the local search method is modified by combining it with DE. The detailed steps are as follows.

Firstly, the DE mutation operator [31] is applied in local search. For , a new mutated individual will be generated according to the following equation: where is the best individual in current generation in the whole population, and are different from and belong to the same group with , and is a random number that is drawn from a uniform distribution .

Secondly, the crossover operator [31] is also introduced after mutation. The two-point crossover operator can be carried out as: Here, is th dimension of th individual ; and are generated randomly from . The mutated individual obtained by mutation operator changes its location according the above equation. This modification not only increases the local search ability but also improves the population diversity to avoid local optimal.

The pseudocode of improved FPA is summarized in Pseudocode 1.