Next Article in Journal
Analysis of Tail Dependence between Sovereign Debt Distress and Bank Non-Performing Loans
Next Article in Special Issue
A Fuzzy ANP-QFD Methodology for Determining Stakeholders in Product-Service Systems Development from Ecosystem Perspective
Previous Article in Journal
Complex Network Analysis of Transmission Mechanism for Sustainable Incentive Policies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of the Product–Service System Configuration Based on a Multilayer Network

1
School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444, China
2
Université de Bourgogne Franche-Comté, UTBM, Pôle Industrie 4.0, Pôle ERgonomie et COnception des Systèmes ERCOS/ELLIADD EA4661, 25000 Belfort, France
3
Shanghai Aerospace Equipments Manufacturer Co., Ltd., Shanghai 200245, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(2), 746; https://0-doi-org.brum.beds.ac.uk/10.3390/su12020746
Submission received: 18 November 2019 / Revised: 8 January 2020 / Accepted: 10 January 2020 / Published: 20 January 2020

Abstract

:
Product–service systems (PSS) accelerate the transition of value creation patterns for manufacturing industries, from product design and production to the delivery of overall solution integrating products and services. Existing PSS configuration solutions provide customers with preferable product modules and service modules characterized by the module granularity. Every service module is essentially a whole service flow. However, the performance of the PSS configuration solution is greatly influenced by service details. In summary, this paper studied the configuration optimization of product-oriented PSS using a fine-grained perspective. A multilayer network composed of (i) a product layer, (ii) a service layer, and (iii) a resource layer was constructed to represent the elements (product parts, service activities, resources) and relationships in PSS. Service activities selection and resource allocation were considered jointly to construct the mathematical model of PSS configuration optimization, thus enabling the calculation of optimizing objectives (time, cost, and reliability) under constraints closer to the actual implementation. The importance degree of service activity was considered to improve the performance of service activities with higher importance. Corresponding algorithms were improved and applied for obtaining the optimal solutions. The case study in the automotive industry shows the various advantages of the proposed method.

1. Introduction

Manufacturing industries are experiencing a transition of value creation patterns from the design and production of products to the design and the delivery of services based on products, to accommodate increasingly global competition and customer-centered business settings [1]. As an overall solution that effectively integrates products and services [2], product–service systems (PSS) [3,4,5] play a vital role in accelerating this transition for manufacturers. By prolonging the service life of products and commoditizing the related services, PSS can change the revenue models and pricing capabilities for manufacturers [6,7]. This exemplifies a new source of competitive advantage and differentiation. Furthermore, compared to traditional business models, PSS has the potential to reduce environmental impact, thus receiving much attention as one strategic alternative for the sustainable development of firms [8,9,10].
PSS can be divided into three types: (a) product-oriented PSS, (b) use-oriented PSS, and (c) result-oriented PSS [5]. In product-oriented PSS, customers buy the product as usual; a service is also sold to customers additionally, which then will be operated on products or provided to customers [11]. This means that the product is the principal purpose, so the service is rendered to help customers achieve a better user experience. This paper focuses extensively on product-oriented PSS.
Services help to attain higher margins, but performing the transition or maintaining the success of service strategy is still a challenge for companies [12]. With the increase in service offers, higher costs or failures in achieving expected returns sometimes occur [13,14]. To overcome this glitch in the profit generation and commercial success, companies should target the role of PSS design [15,16]. Many studies have focused on PSS design to provide support, such as design frames [17,18,19], design methods [20], modularization approaches [21,22], configuration and evaluation methods [23,24,25].
Product–service integration solutions provided to customers include product modules and service modules [26]. There are usually many possible PSS solutions. A PSS configuration is the construction of technical systems by selecting and assembling preferable modules from a predefined product and service library according to customers’ needs under certain constraints [27]. Thus, PSS configuration is an imperative part of PSS design. Studies on PSS configuration have mainly focused on the PSS configuration framework [24,28], customer perception in PSS configuration [29], the extraction of PSS configuration rules [30,31,32,33], and PSS configuration optimization methods [23,34]. Most studies on PSS configuration, especially research on PSS configuration optimization, commonly deal with some predefined level of module granularity. Overall, the product module granularity is introduced to achieve flexible manufacturing systems for product functions realization. With regard to service module granularity, every service module is essentially a whole service flow, depending on the specific module’s granularity level. However, it presently lacks the PSS configuration on a finer granularity level, which can inspire the final configuration’s optimization.
Considering the balance between differentiated customization and mass production, inappropriate product module granularity will increase managerial difficulty and manufacturing cost [35]. Services are generally under-designed in comparison to physical products [36]. Also, service modules are not restricted to physical artefacts. Thus, service modules should be designed considering a finer granularity level, so they can integrate tangible and non-tangible artefacts. Recent studies on service modules in PSS have mainly focused on the identification and evaluation of customer requirements [37,38], modularization and related quantification [39,40], and the evaluation of internal and external performance (such as cost, customer satisfaction, service cycle time, availability, and reliability) [11,41]. Modularization can normalize the PSS configuration and improve its suitability. However, a service module in PSS is essentially a composition of service activities. It is coarse-grained to use a whole service flow as a service module unit in PSS configuration. Therefore, configuring a PSS in finer granularity by using service activity as the minimum configuration unit is proposed and developed in this paper.
There are various methods to optimize the composition of service activities, which can be embodied in the following three approaches [42]:
(1) Service composition solver based on artificial intelligence planning, which can automatically generate service activities composition according to customer demands;
(2) The method based on semantics, studying the service matching and composability;
(3) The service workflow-based semiautomatic service composition method, which is the most widely used one.
In this context, Zhu et al. [43] proposed an ontology-based service process model to support decision-making for maintenance, repair, and overhaul services in the aerospace industry. Cao et al. [44] transformed the service process model into a process structure tree to consider service process structuration and service activities selection simultaneously. Wang et al. [45] advocated a model of service activities composition based on stochastic Petri net. Hu et al. [46] modeled the service process by logical Petri net to study the structural transformation of the service process.
The above-mentioned studies only focused on service activities selection. The impact of resource allocation on service activity has not been considered during the service activities selection. Service activities selection and resource allocation were consequently divided into two independent processes. In a practical situation, the performance of service activity operation will be altered based on different resource allocation instances. It is necessary to combine service activities selection and resource allocation into a whole process during the PSS configuration optimization. Furthermore, most of the previous studies highlighted multi-objective optimization of the overall solution, but high-importance service activities usually exert a greater impact on the solution, compared to low-importance services. Therefore, the importance degree of service activity is also taken into consideration in this paper, to improve the performance of service activities with higher importance.
Based on the above analysis, it was found that: (a) PSS configuration is an imperative part of PSS design; (b) existing PSS configuration solutions are obtained in the module granularity; and (c) the significant impact of service details on the performance of a PSS configuration solution has not been given sufficient attention in recent research studies. PSS configuration optimization is, therefore, studied from a fine-grained perspective in this paper. The optimization is performed by combining service activities selection and resource allocation. As a result, the PSS configuration solutions provided to the customer will be more refined and precise. In particular, the optimizing objectives, which present the performance of the solutions, will be calculated closer to the actual implementation.
The complex network is chosen as the framework for the PSS modeling because its structural properties enable the exploration of the interactions between products, resources, and service activities. The mathematical characteristics of complex networks are very useful for analyzing the characteristics of complex systems. However, it is difficult for a general complex network to effectively reveal the multilayered or multidimensional relationships in a complex system when the scale of the system increases substantially. Many approaches have been proposed to enrich the representation of networks. Hypernetworks [47,48,49] are developed based on the definition of the hypergraphs to generalize the multidimensional relationships, in which multi-elements are connected rather than only two elements. Multilayer networks [50] are applied to describe the complex structures of systems, which include different types of relationships between their components. Originally, multilayer networks were proposed to express different types of relationships between one single type of elements in different layers [51,52]. More recently, they have been used to model other systems, including those with different types of elements. Omodei et al. [53] proposed a method based on the multilayer networks of citations and disciplines to assess the interdisciplinary importance of scholars, institutions, and countries. Additionally, Pasqual and Weck [54] introduced a multilayer network model integrating three coupled layers (product layer, change layer, and social layer) to analyze and manage engineering change propagation. Leng and Jiang [55] also offered a multilayer network model comprising processes, machines, and work in processes to assess dynamic scheduling in a radio frequency identification-driven discrete manufacturing system.
Based upon the comparison of the general complex network, hypernetwork, and multilayer network, it is obvious that a multilayer network can perfectly accommodate multiple heterogeneous elements and provide a better description of relationships between heterogeneous elements than a general complex network or hypernetwork. Therefore, a multilayer network, as a special type of complex network, is applied for the modeling of PSS in this paper.
The organization of this paper is as follows. In Section 2, this paper’s research problem is described and the framework is then introduced. The PSS multilayer network model is constructed in Section 3. Section 4 establishes the mathematical model of PSS configuration. This section also delineates the improved optimization algorithms and their applications for obtaining the optimal solutions. Moreover, a case is studied in Section 5 to discuss the feasibility and effectiveness of this method. Section 6 highlights the benefits and the contributions to the sustainability of the proposed approach. Finally, conclusions are stated in Section 7.

2. Problem Description and Framework

In product-oriented PSS, products are delivered to customers as tangible objects in a definite time, just after the PSS configuration solution is accepted. However, that is not the case with services. Services are provided to customers only when the customer requests or the service premises are satisfied and operation of services is affected by the real-time conditions of PSS. Service details, which include service activities selection and resource allocation, are consequently not strictly determined when the PSS configuration solution is accepted. Contracts to implement corresponding service functions are provided to the customer. But service details definitely impact the performance of the PSS configuration solution. Service time depends on service activities selection, cost and reliability are accurately influenced by resource allocation. Thus, the PSS configuration should be optimized by combining service activities selection and resource allocation. In this way, the PSS solutions obtained after optimization will be more precise, and the results based on the optimizing objectives (executing time, executing cost, and executing reliability) will be closer to the actual implementation.
The framework of this fine-grained configuration optimization of the product-oriented PSS is shown in Figure 1. A PSS model is constructed based on a three-layer network, which contains product and services selected by the customer, and the resources in the system. A node in every layer (product layer, service layer, resource layer) of the multilayer network respectively represents an element (product, service, resource). The relationships between homogeneous elements in a layer and the relationships between heterogeneous elements in different layers are represented as edges. With appropriate decomposition, PSS can be clearly described by the proposed multilayer network. Elements and relationships in PSS can be expressed by connection matrixes, which make the description of the mathematical model and realization of optimization more convenient. By using this PSS multilayer network model, the importance degrees of service activities are analyzed based on the attributes of the multilayer network. A mathematical model based on the connection matrixes is constructed to express the objectives and constraints of PSS configuration optimization, and the corresponding solution algorithm is proposed to solve the mathematical model of optimization.

3. PSS Multilayer Network Model

The PSS multilayer network ( N E T P S S ) is composed of the product layer ( N e t P ), service layer ( N e t S ), resource layer ( N e t R ), and set of inter-layer edges ( E I ) which connect layers to each other. Therefore, there are three types of nodes and five kinds of edges considered in N E T P S S . Specifically, nodes in N e t P , nodes in N e t S , and nodes in N e t R are denoted as sets P , S , and R , respectively; edges in N e t P , edges in N e t S , edges in N e t R , inter-layer edges between N e t S and N e t P , and inter-layer edges between N e t S and N e t R are denoted as sets E P , E S , E R , E P , S , and E R , S , respectively. The inter-layer edges between N e t R and N e t P are not considered in this model.
Then, N E T P S S is described by Formula (1):
N E T P S S = V ( P , S , R ) , E ( E P , E S , E R , E P , S , E R , S ) .
A product can be decomposed into parts according to required granularity. N e t P is formed by part nodes (p) connecting to each other by interrelationships between parts. N e t S consists of service flows, each of which is described by a service trigger node (ss), numbers of service activity nodes (sa), and a service end node (se) logically associated with each other. N e t R contains resources required by service activities, including personnel resources (pr) and equipment resources (er), as resource types are denoted as kr. E I contains E P , S , E R , S , and the weights on inter-layer edges. Detailed descriptions of N e t P , N e t S , N e t R , and E I are given in the following subsections.

3.1. Product Layer

A node p i in N e t P is described by Formula (2):
p i = p _ B A S i , p _ R E C i , p _ S T i ,
where p _ B A S i is the set of basic information, including information such as name, quality, and material. p _ R E C i is the set of accepted service records. p _ S T i is the set of working statuses, including various data of physical status, material status, and environmental status. All statuses vary with time t. For any p i whose status data cannot be monitored, p _ S T i = .
There are four types of interactions among parts [56]: spatial-type interaction, energy-type interaction, information-type interaction, and material-type interaction. These interactions could be extracted from the product design manually. If p i provides spatial reference (adjacency or orientation) to p j , then s e i j = 1 . If p i provides or transmits energy to p j , then e e i j = 1 . If p i provides or transfers information to p j , then i e i j = 1 . If p i determines or impacts the material selection of p j , then m e i j = 1 . If there are one or more interactions from p i to p j , then e i j P = 1 , denotes the directed edge from p i to p j ; otherwise, e i j P = 0 . Obviously, e i j P e j i P .
In order to describe the overall strength of the relationship between p i and p j , the weight on e i j P is represented by the value w i j P , following Formula (3):
w i j P = w s s e i j + w e e e i j + w i i e i j + w m m e i j ,
where w s , w e , w i , and w m , respectively, represent the proportion of s e i j , e e i j , i e i j , and m e i j in the edge e i j P . w i j P is defined in [0,1]. If the four types of interactions are equal in importance, then w s = w e = w i = w m = 0.25 .
N e t P is consequently represented by Formula (4):
N e t P = P , E P ; P = { p 1 , p 2 , , p n } ; E P = { e i j P } ; W P = { w i j P } .

3.2. Service Layer

Services in a product lifecycle can be classified into two types: production services and product services. All services in N e t S are product services.

3.2.1. Service Trigger Node

A node s s i in N e t S is as described by Formula (5):
s s i = s s P i , T r i i , s _ s t i ; T r i i = { t r i i ( p j ) | p j s s P i } ; s _ s t i = { o n ( s t a r t _ t i , d d l _ t i ) o f f } ;
where s s P i is the set of parts whose working statuses are monitored by s s i . T r i i is the set of triggering conditions of s s i , t r i i ( p j ) is the triggering condition for working statuses p _ S T i of p i to trigger s s i . s _ s t i denotes the status of s s i : if triggering conditions are satisfied, the status will be o n ( s t a r t _ t i , d d l _ t i ) , where s t a r t _ t i and d d l _ t i are, respectively, the start time and the deadline of the entire service set according to customer requirements; otherwise, the status will be o f f .

3.2.2. Service Activity Node

The function of service activity is to change the status of the service activity receiver [57]. Service activity is operated by the specified resources and executed on the specified component. By decomposing the service flow at specified granularity, the corresponding set of service activities could be obtained. A service activity node s a i in N e t S is as described by Formula (6):
s a i = a _ R E C i , T i , C i , R E i , O i , a _ s t i ; a _ s t i = { o f f o n ( b t i , f t i ) } ;
where a _ R E C i is the set of service execution records of s a i . T i and C i are the theoretical executing time and fixed cost for operating s a i , which are extracted from historical data. R E i denotes the reliability of that s a i which can be completed in time T i . It is represented by the ratio between the number of times that s a i has been completed in time T i and the total number of times that s a i has been executed. R E i ( 0 , 1 ] . O i is set of service objects of s a i . a _ s t i denotes the status of s a i : o f f is idle status, o n ( b t i , f t i ) denotes working status. b t i is the beginning time of s a i , ending time is defined as f t i = b t i + T i .

3.2.3. Service End Node

A service end node s e 0 is put in N e t S to distinguish whether the service flow is finished at the end of s a i . An edge between s a i and s e 0 means that s a i is the last service activity in the service flow.

3.2.4. Edges in Service Layer

Regarding s s j and s e 0 as the activities of “trigger service” and “end service”, respectively, then s s j , s e 0 , and s a i are homogeneous. Edges between s a i and s a j , edges between s a i and s s j , and edges between s a i and s e 0 are all denoted by e S ; e i j S e j i S , e i i S = 0 . w i j S defined in the interval [0,1] denotes weight on e i j S . w i j S = 0 when e i j S = 0 . The description of an edge e S between two nodes considers the following cases:
  • Between s s j and s a i : if s a i is the first one executed in the service flow triggered by s s j , then e j i S = 1 , b t i = s t a r t _ t j . s a i can participate in different service flows triggered by s s j , denoted as e i j S = 1 with w i j S 0 , s a i will not output any weight to s s j ; e i j S = 0 between s s i and s s j ;
  • Between s e 0 and s a i : if s a i is the last one in the service flow being executed, then e i 0 S = 1 ; e 0 i S = 0 for any s a i and e i 0 S = e 0 i S = 0 between s e 0 and s s i ;
  • Between s a j and s a i : if s a i and s a j are in same service flow and are executed in time sequence, then e i j S = 1 . Otherwise, e i j S = 0 .
There are three types of correlations among service activities [58]: correlation of function similarity, correlation of function complementary, and no correlation. If two service activities have a correlation of function similarity, their function might be identical, inclusive, or partially similar. If two service activities have a correlation of function complementary, the composable situation might be accurate composition, inclusive composition, blocked composition, or crossed composition.
According to the above correlations, logical relationships among service activity nodes are divided into serial relationship, parallel relationship, and selective relationship. Two activities in a serial relationship are sequentially executed in chronological order. Two activities in a parallel relationship are executed at the same time, and only when they are fully completed will the execution of subsequent activities begin. In a selective relationship, there are multiple activities for selection: if only one of them is executed, then the execution of subsequent activities will begin.
The in-degree and out-degree of s a i in N e t S are denoted as i n w i and o u t w i which are described by Formulas (7) and (8):
i n w i = j e j i S w j i S ,
o u t w i = j e i j S w i j S .
Then, the description of the three types of logical relationships is expanded as follows:
  • If s a h and s a l are in a serial relationship, and i n w h = x , then e h l S = 1 , w h l S = x , b t l = f t h ;
  • If n nodes ( s a 1 , s a 2 , , s a n ) are in a selective relationship, and i n w h = x , then e h i S = e h j S = e i l S = e j l S = 1 , w h i S = w h j S = w i l S = w j l S = x , o u t w l = x . In this case, o u t w h > i n w h and o u t w l < i n w l . If s a i is selected, then b t i = f t h , b t l = f t i ;
  • If n nodes ( s a 1 , s a 2 , , s a n ) are in a parallel relationship, and i n w h = x , then e h i S = e h j S = e i l S = e j l S = 1 , w h i S = w h j S = w i l S = w j l S = x / n , o u t w l = x , and b t j = b t i = f t h , b t l = max ( f t i , , f t j ) . In this case, o u t w h = i n w h and o u t w l = i n w l .
Connections of these logical relationships are shown in Figure 2. If s a h is replaced by a service trigger node, then x = 1 .
Then, the logical relationship between nodes in N e t S can be distinguished by values of i n w and o u t w . However, if s a h has multi inputs and multi outputs, as depicted in Figure 3a, the logical relationship between s a h and other nodes will be uncertain. The virtual service activity node s a m is applied to prevent that. Then, s a h will become a node with a single input or single output, as shown in Figure 3b,c. There is not any essential service content in a virtual service activity node, so T m = 0 , C m = 0 , R E m = 1 , a _ R E C m = , O m = , and s a m does not require any resources.
If s a h connects the other nodes, which are in selective relationships and parallel relationships, intricately, it will be impractical to distinguish different relationships by values of i n w h and o u t w h . As shown in Figure 4a, when w h i S = w h j S = w h k S = w h g S = x / 2 , relationships among s a i , s a j , s a k , and s a g cannot be distinguished directly. In order to solve this situation, there will be two situations to integrate in the virtual service activity nodes s a m and s a n , as follows:
  • w h m S = w h n S = w m i S = w m j S = w n k S = w n g S = x / 2 means that s a i and s a j are in a selective relationship; s a k and s a g are in a selective relationship. Meanwhile, s a i (or s a j ) and s a k (or s a g ) are in a parallel relationship, as expressed in Figure 4b.
  • w h m S = w h n S = x , w m i S = w m j S = w n k S = w n g S = x / 2 means that s a i and s a j are in a parallel relationship; s a k and s a g are in a parallel relationship. Meanwhile, two parallel clusters are in a selective relationship, as indicated in Figure 4c.
Assuming that the number of s a and s s in N e t S are indicated, respectively, by δ and γ, and there is only one s a , then the description of N e t S is consequently shown by Formula (9):
N e t S = S , E S ; S = { s 0 , s 1 , s 2 , , s δ + γ } = { s e 0 , s a 1 , s a 2 , , s a δ , s s δ + 1 , s s δ + 2 , , s s δ + γ } ; E S = { e i j S } ; W S = { w i j S } .

3.3. Resource Layer

The node p r i or e r i in N e t R is described by Formula (10):
r i = p r i e r i = r _ R E C i , c i , r _ s t i ; r _ s t i = i d l e b u s y ;
where r _ R E C i is the set of records of services operated by r i , c i is the cost for r i working in one unit of time, r _ s t i denotes the status of r i . If r _ s t i = i d l e , then r i could accept a new assignment of the service activity; if r _ s t i = b u s y , it means r i is in working status and will not accept any new assignments until the service activity r i is working on is finished. If there are not any resource contents in k r j , then r _ R E C j = , c j = 0 , r _ s t j = .
If there is a subordinate relationship between p r i e r i and k r j , which means that the resource type of p r i e r i is k r j , then e j i R = 1 . This denotes that there is a directed edge from the father node k r j to the child node p r i e r i . Otherwise, e j i R = 0 . For any p r i e r i , j e j i R = 1 .
Assuming that quantities of p r , e r and k r in N e t R are indicated by α , β , and b , respectively, then the description of N e t R is shown by Formula (11):
N e t R = R , E R ; R = { p r 1 , , p r α , e r α + 1 , , e r α + β , k r α + β + 1 , , k r α + β + b } = { r 1 , r 2 , , r α + β + b } ; E R = { e j i R }   , ( i = 1 , 2 , , α + β ; j = α + β + 1 , , α + β + b ) .

3.4. Inter-Layer Edges

The service content is determined by the working status of products, while service activities are executed on products by required resources. In the PSS multilayer network, some inter-layer edges exist between N e t S and N e t P , the rest exist between N e t S and N e t R .

3.4.1. Inter-Layer Edges Between N e t S and N e t P

For p i and s a j , if p i O j , then e j i S P = 1 , which denotes the directed inter-layer edge from s a j to p i . This means that p i is one of the service objects of s a j . Otherwise, e j i S P = 0 . Meanwhile, e S P = 1 only makes sense between p and s a ; for s s or s e , e S P = 0 .
Therefore, the set of inter-layer edges between N e t S and N e t P is as described by Formula (12):
E P , S = { e j i S P } .

3.4.2. Inter-Layer Edges Between N e t S and N e t R

Inter-layer edges between N e t S and N e t R include directed inter-layer edges from service nodes to resource nodes and directed inter-layer edges from resource nodes to service nodes. Their descriptions are depicted below:
  • The execution of service activities demands the resources be in corresponding specific kinds and quantities. If s a j demands the resources whose type is k r i , then e j i S R = 1 , which denotes a directed inter-layer edge from s a j to k r i and the quantity demand is represented by integer weight w j i S R . Otherwise, e j i S R = 0 , w j i S R = 0 . Meanwhile, e S R = 1 only makes sense between k r and s a ; for other nodes, e S R = 0 ;
  • If e j h S R = 1 , e h i R = 1 , and r _ s t i = i d l e , it means that p r i e r i is assignable to s a j , then e i j R S = 1 , which denotes a directed inter-layer edge from p r i e r i to s a j . Otherwise, e i j R S = 0 . Meanwhile, e R S = 1 only makes sense between p r e r and s a ; for other nodes, e R S = 0 . The reliability of s a j being completed in time T j by p r i e r i is denoted by weight w i j R S in [0,1].
Every p r i has a different knowledge level and proficiency level about different s a . Assuming P A i denotes the number of times that s a j has been completed in the time T j by p r i , P B i denotes the total number of times that p r i has participated in the execution of s a j , then w i j R S between p r i and s a j is represented by P A i / P B i . Every e r i has a constant function in a different s a . Assuming E A i denotes the number of times that any s a has been completed in time by e r i and E B i denotes the total number of times that e r i has participated in the execution of any s a , then w i j R S between e r i and s a j is represented by E A i / E B i .
The set of inter-layer edges between N e t S and N e t R , and the set of weights on these inter-layer edges are described by Formula (13):
E R , S = { e i j R S , e j i S R } ; W R , S = { w i j R S , w j i S R } .
Therefore, E I is as given by Formula (14):
E I = E P , S , E R , S .

4. PSS Configuration Optimization

In most practical cases, many services are usually triggered at the same time. Consequently, the configuration optimization should be designed by considering multi-services. After any s s is triggered, service activities selection and resource allocation should be accomplished under given constraints to obtain the PSS configuration under constraint conditions.
Let us assume that the number of triggered s s is σ up to the point of time t . If s s φ is triggered and is represented by u φ following Formula (15), then the situation of all s s is represented by trigger matrix [ u δ + 1 , u δ + 2 , , u δ + γ ] with σ = φ = δ + 1 δ + γ u φ :
u φ = { 1 , s s φ   i s   t r i g g e r e d 0 , o t h e r w i s e .
As a result, service activities selection is regarded as a selection of the paths with multi s s as start points in N e t S . It is necessary to consider the executing time, cost, and reliability of s a in the selective relationship. The allocation of resources is viewed as resource allocation for multitasks (service activities), in which the cost and reliability of different r i are considered. Therefore, the services solution is optimized with the goal of executing time minimization, cost minimization, and reliability maximization.

4.1. Importance Degree of Service Activity

When a certain type of resource is required by several service activities simultaneously, then those resources with higher reliability should be allocated to more important service activities. Product parts are the objects of service activities; hence, the importance degree of product parts exerts a decisive influence on the importance degrees of service activities. Consequently, the importance degree of product parts is taken as the base of calculation for the importance degree of service activities, which provides a theoretical reference for the resource allocation.
Today, network-based importance measurement is thoroughly developed and widely applied. Various measure methods can be summarized as:
(1) Importance calculation methods based on neighbor relationships between nodes, such as degree centrality or semi-local centrality [59];
(2) Path-based importance calculation methods, such as closeness centrality and betweenness centrality;
(3) Importance calculation methods based on eigenvectors, such as search algorithm LeaderRank [60] and SRank [61] based on algorithm PageRank [62];
(4) Importance calculation methods based on node removal and contraction, in which the importance degree of nodes is measured by the variation in network attributes after node removal or contraction [63].

4.1.1. Importance Degree of Product Parts

Importance calculation methods based on eigenvector are applied in this paper. In the concept of PageRank and LeaderRank, the importance of a page in a network is determined by the number and the quality of other pages pointing to it. However, the importance of p i in N e t P is ascertained by the number and quality of other nodes it points to. This makes PageRank and LeaderRank unsuitable for calculating the importance degree of product parts. Referring to the weighted LeaderRank [64], PartRank is proposed to calculate the importance degree of product parts. Details are illustrated as follows.
A ground node p g , which connects with every part node by a bidirectional edge, is introduced. Then e g i P = 1 and e i g P = 1 , i g . Based on the concept of weighted LeaderRank, nodes with higher input intensity embody a higher probability of being visited by p g . Inversely, p i with a higher output intensity is more likely to visit p g in N e t P . w i g P is consequently described by formula (16), where θ is a free parameter, n is total number of nodes in N e t P (excluding p g ). w g i P is described as Formula (17):
w i g P = ( h = 1 n + 1 e i h P w i h P ) θ ;
w g i P = h = 1 n e h i P w h i P h = 1 n e h i P .
The importance degrees are calculated iteratively until they reach a steady state. Assuming that calculation reaches a steady state at step k , then the importance degree of p g will be allotted to every p i . The importance degree of p i is denoted as P R i , calculated as Formula (18),
{ P R g ( 0 ) = 0 P R i ( 0 ) = 1 , i g P R i ( t ) = j = 1 n + 1 e i j P w i j P k j i n P R j ( t 1 ) , t = 1 , 2 , , k P R i = P R i ( k ) + w i g P k g i n P R g ( k ) ,
where k j i n = h = 1 n + 1 ( e h j P w h j P ) is the input intensity of p i in N e t P .

4.1.2. Importance Degree of Service Activities

Knowing that the service object of sa is p and if s a is operated on p with a high importance degree, a designer should pay more attention to that s a . C N I j denotes the cross layer importance degree of s a j , which is related to the product parts on which s a j operates. The value of C N I j is determined by p i with the highest P R i in the set of service objects of s a j . C N I j is consequently calculated by Formula (19):
C N I j = max 1 i n ( e j i S P P R i ) .
After being normalized, the importance degree of s a j is denoted as D j = C N I j / max , where max is the maximum of all C N I j in N e t S . Therein, the virtual service activity node does not have any service object, so the importance degree of every virtual service activity node is consequently 0.

4.2. Mathematical Model

Define the binary variable y j to indicate whether s a j is selected, and define the binary variable a x j i to indicate whether to allocate r i to s a j , by the Formulas (20) and (21), respectively:
y j = { 1 , s a j   i s   s e l e c t e d 0 , o t h e r w i s e ,
a x j i = { 1 , a l l o c a t e   r i   t o   s a j 0 , o t h e r w i s e .
The mathematical model of PSS configuration optimization is constructed by Formulas (A1)–(A15) in Appendix A. The objective functions are given by Formulas (A1)–(A3) and constraints are given by the Formulas (A4)–(A15).
The objective function of executing time minimization is described by Formula (A1). The executing times of service activities are assumed to be fixed parameters, and they are not impacted by resource allocation. T S denotes the sum of executing times for all selected s a .
The objective function of cost minimization is described by Formula (A2). The executing cost of s a j includes the fixed cost C j and the work costs of all resources allocated to s a j , which are related to the executing times of s a j . C S denotes the sum of executing costs of all selected s a .
The objective function of reliability maximization is described by Formula (A3). The executing reliability of s a j is determined by w i j R S of all r i that are allocated to s a j . Therefore, a logic function z j i , which depends on x j i as the formula (A12) is defined. Then, e x r e j , denoting the executing reliability of s a j , is described as e x r e j = i = 1 α + β z j i . The target of reliability is always to maximize the executing reliability of the whole solution, which is related to the e x r e j of all selected s a j . However, if a certain type of resource is required by several service activities with overlapping service times, the allocation result of that type tends to make a s a j with a higher importance degree for reaching a higher e x r e j . Therefore, A V R E S is defined as the average executing reliability of all selected s a j , which is shown by Formula (22):
A V R E S = j = 1 δ ( y j D j e x r e j ) j = 1 δ y j D j .
The e x r e j of a s a j with higher D j has a greater impact on the value of A V R E S .
When y j = 1 and r i is assigned to s a j , allocating r i to s a j might be feasible. Therefore, x j i is defined by the formula (A13). Here, x j i = 1 means r i is allocated to s a j successfully, x j i = 0 means r i is not allocated to s a j . Constraints for y j and a x j i are shown by Formulas (A14) and (A15). Constraints for the three objective functions are developed as indicated below:
  • All selected service activities must be part of the construction for the triggered service flows. They are represented by Formula (A4);
  • Every service flow starts with the triggered s s φ , and the out-degree of s s φ is 1. This is expressed by Formula (A5);
  • The number of triggered s s is σ . All triggered service flows end at s e 0 . Consequently, the in-degree of s e 0 is σ . This is shown by Formula (A6);
  • If all selected service activities constitute complete service flows, there will not be any selective relationship between selected s a i and s a j . As a result, the sum of weights on edges from triggered s s φ and selected s a j to s a i will equal to the sum of weights on edges from s a i to s e 0 and selected s a j . This is shown by Formula (A7);
  • For every triggered s s φ , s t a r t _ t φ = t . Every triggered flow must be completed within the service time range [ s t a r t _ t i , d d l _ t i ] requested by a customer, as illustrated in Formula (A8);
  • The resources allocated to the selected s a j must satisfy the constraints as the resource type and the quantity demands of s a j . This constraint is denoted by Formula (A9);
  • Every r i can only participate in a maximum of one service activity at the same time. For all service activities that r i  is allocated to, their service time intervals cannot overlap. Defining oljh indicates whether there is any overlap between the service time intervals of s a j and s a h , as shown in Formula (A11). That constraint is described by Formula (A10).

4.3. Obtaining the Optimal Solutions

To obtain the optimal solutions by the trigger matrix [ u δ + 1 , u δ + 2 , , u δ + γ ] and the proposed mathematical model, the solving procedures are as follows.
Step 1:
Obtaining S A S i as the service activities solution set for s s i , the number of solutions is denoted as n i . If u i = 1 , obtain S A S i = { S A 1 i , , S A x i i , , S A n i i } by Depth-First Traversal (DFT), wherein S A x i is the set of all service activities selected in the x i time traversal; if u i = 0 , s s i is not activated, then n i = 1 , S A S i = { S A 1 i } , wherein S A 1 i = . N e t S is different from ordinary complex networks, and as a consequence, the traditional DFT is unsuitable for N e t S . Depth-First Traversal based on service layer (SL-DFT) is designed and detailed in Section 4.3.1.
Step 2:
Solutions of service activities for each triggered s s i are arranged by permutation and combination. Obtaining a solutions set of service activities as S A S a l l = { S A 1 a l l , S A 2 a l l , , S A n a l l a l l } for system, n a l l = v n v , S A x a l l = { S A x δ + 1 δ + 1 , S A x δ + 2 δ + 2 , , S A x δ + γ δ + γ } .
Step 3:
Using the Pareto-improved multiple-objective discrete cuckoo search algorithm (Pareto-IMODCS) [23] to obtain the set of optimal solutions of resource allocation for each S A i a l l as R A S i = { R A 1 i , , R A m i i } , the number of solutions is denoted as m i . Pareto optimality and Pareto-IMODCS are detailed in Section 4.3.2 and Section 4.3.3, respectively.
Step 4:
Merging all R A S i into one as R A S a l l = { R A 1 δ + 1 , , R A m δ + 1 δ + 1 , , R A 1 δ + γ , , R A m δ + γ δ + γ } , filter R A S a l l to obtain the set of non-dominated solutions as R A S p a r e t o based on Pareto optimality. R A S p a r e t o is the final set of optimal solutions.

4.3.1. Obtaining Solutions of Service Activities by SL-DFT

Obtaining S A S v by SL-DFT involves searching all paths from s s v to s e 0 in N e t S . Assuming the current visited node is s x :
  • If i = 0 δ e x i S = 1 and e x y S = 1 , then add node s y into the path and continue searching from s y ;
  • If i = 0 δ e x i S > 1 and o u t w x > i n w x , nodes satisfying e x i S = 1 are in a selective relationship, then select one of these to visit. Assuming that the selected node is s y , continue to search from s y until all paths from s y to s e 0 are obtained, then backtrack to s x , selecting an unvisited node from nodes satisfying e x i S = 1 to visit. This new search forms a new path, so nodes visited in the search from s y can still be visited. Repeat this step until all nodes satisfying e x i S = 1 are visited;
  • If i = 0 δ e x i S > 1 and o u t w x = i n w x , it means nodes satisfying e x i S = 1 are in a parallel relationship. Then, select one to visit and continue searching until the visited node s z satisfies i = 0 δ e i z S > 1 and i n w z = o u t w z , which means s z is the confluence node of these parallel paths. Then, backtracking to s x , select an unvisited node from nodes satisfying e x i S = 1 to visit. Repeat this step until nodes satisfying e x i S = 1 are all visited. Add all nodes visited in this step into the path and continue to search from s z .
After all paths from s x to s e 0 are obtained, if s x is a service trigger node, then end the search. Otherwise, backtrack to the node before s x . Based on the above search, service activity compositions obtained by SL-DFT satisfy the constraints as given by Formulas (A4)–(A7).

4.3.2. Pareto Optimality

In the optimizing process, the solution with smaller executing time, smaller executing cost, and greater executing reliability is considered better. However, objectives in one solution might conflict with each other in reality. Traditional linear weighting methods are unsuitable for solving such multi-objective optimization, which could be solved based on Pareto optimality [65].
If x i and x j satisfy that T S ( x i ) T S ( x j ) , C S ( x i ) C S ( x j ) and A V R E S ( x i ) A V R E S ( x j ) , and at least one of these inequalities is a strict inequality, then x i dominates x j . If x m is not dominated by any x , then x m is a non-dominated solution, which is the Pareto optimal solution.

4.3.3. Obtaining Optimal Solutions of Resource Allocation by Pareto-IMODCS

Let us assume that the number of s a in S A a l l is v regardless of s s , s e , and virtual service activity nodes. The types and quantities of resources that every s a demands can be determined based on the multilayer network model. Then, each possible resource allocation solution is represented by a cuckoo nest: the dimension L ( x K x N x K x ) of nest represents the N x K x th resource instance for the K x th resource type demand of the x th s a in set S A a l l ; various instances in each dimension are all represented by decimal integers. The subscript number of the service activity node corresponding to the x th s a is denoted as n s a ( x K x ) . The subscript number of the resource type node corresponding to the K x th resource type is denoted as   n r ( x K x ) . n s a ( x K x ) and   n r ( x K x ) facilitate in decoding the nest into the corresponding resource allocation solutions, which satisfy the constraint defined by Formula (A9).
The original Pareto-IMODCS algorithm is detailed in reference [23]. If multiple s a in S A i a l l demand same type of resources in the same time interval, the solution obtained by the original Pareto-IMODCS algorithm might allocate the same resource instance to multiple s a simultaneously. This does not comply with the constraint described by Formula (A10). To avoid this, a “nest pre-processing” step is added to Pareto-IMODCS. The specific procedures of the improved Pareto-IMODCS algorithm to obtain Pareto optimal solutions for each S A i a l l are shown in Figure 5.
The procedure “nest pre-processing” is specified as: for any two dimensions as L ( h i j ) and L ( l m n ) , when n r ( h i ) = n r ( l m ) and L ( h i j ) = L ( l m n ) = x , it means that the same resource instance is allocated to multiple dimensions. If n s a ( h i ) = n s a ( l m ) or o l n s a ( h i )   n s a ( l m ) = 1 , it denotes that service activities related to L ( h i j ) and L ( l m n ) are the same one or their working time intervals overlap, then L ( h i j ) = x and L ( l m n ) = x + 1 .

5. Case Study

In this section, a case study for the PSS of an automobile enterprise is provided to validate the proposed multilayer network-based PSS model and the fine-grained PSS configuration optimization method. A product was delivered to a customer with a contract to implement eight selected services (running system maintenance, steering system maintenance, braking system repair, etc.). These services could be decomposed as service flows consisting of service activities. For example, running system maintenance was decomposed as the maintenance of axles, maintenance of suspensions, maintenance of shock absorbers, maintenance of wheels, and maintenance of frames. As indicated, the PSS of the automobile enterprise was modeled as a multilayer network, which included 193 product nodes, 89 service nodes (1 s e , 80 s a , and 8 s s ), and 140 resource nodes (70 p r , 50 e r , and 20 k r ).
When services are triggered, a set of PSS configuration solutions should be provided to the customer for making a final decision. The customer’s desired performances for the final solution were the minimal executing time and cost and maximal executing reliability. The designer had a PSS configuration optimization task at time t , which contained the configuration of services as running system maintenance ( s s 81 ) and steering system maintenance ( s s 82 ). The network status of the case at time t is modeled in Figure 6, in which the sparse matrixes of edges denoting corresponding relationships among nodes in multilayer network are detailed in Figure 7.
According to the PSS configuration optimization task, the trigger matrix at time t was [1 1 0 0 0 0 0 0]. Structures of triggered service flows in service layer are detailed in Figure 8. The executing time, cost, and reliability of related s a j are detailed in Table 1. It is evident that s a 1 and s a 5 were virtual service activity nodes. The cost for all r i working in one unit of time is shown in Figure 9.
In the above-mentioned network status and solution constraints, all paths from triggered s s 81 to s e 0 and all paths from triggered s s 82 to s e 0 were searched by SL-DFT. Then, one from both sets of paths was selected to form service activity compositions. Pareto optimal solutions of resource allocation for every service activity composition were obtained by using the mentioned improved Pareto-IMODCS algorithm. Relevant parameters of the algorithm were as follows: the number of initial nests was 200, the probability of cuckoo eggs being found was 0.25, the maximum number of iterations for a nest was 3000, and the number of times of an iteration run was 10. Then, solutions of all compositions were filtered according to Pareto optimality. As a result, 14 non-dominated solutions were obtained, which were recorded as set A (Table 2).

5.1. Comparison with Method For Selecting Service Activities Before Resource Allocation

To validate the effectiveness of combining service activities selection and resource allocation, the same case was solved with another method as a contrast. In summary, the details were as follows:
Firstly, all paths from triggered s s to s e 0 were searched by SL-DFT and all service activity compositions were obtained. Next, all compositions were filtered by objective functions as executing time, cost, and reliability to obtain Pareto optimal service activity compositions. Then, Pareto optimal solutions of resource allocation for each composition were acquired using the improved Pareto-IMODCS with same parameters as last time. Finally, all solutions were filtered according to Pareto optimality. As a result, six non-dominated solutions were recorded as set B (Table 3).
By analyzing solutions in Table 2 and Table 3, it can be found that set A and set B share four common solutions. Nevertheless, set A also contains other solutions in which service activity compositions are different from solutions in set B. The other two solutions in set B (solutions B5 and B6) are dominated by solution A3. Obviously, combining the service activities selection and resource allocation effectively enriches the diversity of the optimal solution set by retaining more non-dominated solutions. The last ones might have been hindered by using the method where service activities selection finished individually before resource allocation.

5.2. Comparison with the Method Using Different Objective Functions

By calculating the importance degrees of service activities related to triggered s s 81 and s s 82 , it was found that the importance degree of s a 9 is significantly higher than other service activities. Reverting to the actual setting of the case, s a 9 was frame maintenance in the running system maintenance. The frame is the base of the whole automobile assembly, which makes frame maintenance more important compared to other activities in the maintenance of the running system and the steering system.
To validate the maximization of A V R E S , an objective function which improves the executing reliability of service activities with higher importance degrees, the same case was solved in another way, as follows.
The value 1 was given to the weights of all e x r e j in the objective function of executing reliability, regardless of the importance degrees of the service activities. The altered objective function of the executing reliability is shown in Formula (23), where E X R E S denotes the overall executing reliability of all selected activities. As a result, 13 Pareto optimal solutions were obtained, which were recorded as set C.
max E X R E S = j = 1 δ y j e x r e j .
The executing reliabilities of s a 9 ( e x r e 9 ) in different solutions with set A and set C are shown in Figure 10. The values of e x r e 9 in set A are all higher than those in set C. The values of e x r e 9 in solution A3 and C11 are the highest in set A and set C, respectively. Service activity compositions in solution A3 and C11 are the same. Comparing and analyzing these two solutions shows why maximizing A V R E S , as one of the objective functions, could facilitate more beneficial solutions, in which service activities with higher importance degrees have higher executing reliabilities.
The Gantt chart of selected service activity nodes in solution A3 and C11 is depicted in Figure 11. From the perspective of the service time interval, s a 9 overlaps with s a 19 and s a 22 , s a 19 overlaps with s a 8 , s a 8 overlaps with other nodes. Situations of service time interval overlap between service activities are intricate. Consequently, resource allocation of one activity might influence another one and indirectly impact other activities.
The types and quantities of resources demanded by s a 8 , s a 9 , s a 19 , and s a 22 are shown in Table 4. The actual situation of resource allocation about the k r related to s a 9 in solution A3 and C11 is illustrated in Table 5. w i j R S on the edges that connect s a j and r i (all p r i and e r i that are related to k r 124 , k r 126 , k r 134 , and k r 136 ) are detailed in Figure 12 and Figure 13.
The decision process in multi-objective optimization is complicated. In solution C11, r i with the highest w i 9 R S is not allocated to s a 9 in the process of overall executing reliability maximization. In solution A3, p r 28 , p r 46 , e r 90 , and e r 98 are all allocated to s a 9 . w i 9 R S of these four resources is the highest among the same type of resources. In turn, the e x r e 9 in solution A3 is higher than solution C11. For the remaining three s a , for whom the importance degrees are not conspicuously higher than each other, the resource allocation is not only based on w i j R S , but also on the cost (as shown in Figure 9). By including the importance degree of s a j with the objective function of executing reliability as weight, it enables more important s a j to have higher executing reliabilities, by allocating r i with the highest w i j R S to this s a j . If the importance degree of s a j is significantly higher, the impact on the overall executing reliability and executing cost will be tolerable.

6. Discussion

As an imperative part of PSS design, PSS configuration is usually studied in a predefined level of module granularity [21,22]. The service modules in existing studies of PSS [23,34] are essentially whole service flows, which can be decomposed into service activities. In the studies of service activity composition [42,45], the impact of resource allocation on service activity is not considered during the service activities selection. This paper studied the PSS configuration in a finer granularity level to inspire the final configuration’s optimization. By combining service activities selection and resource allocation, the proposed approach enriched the diversity of the optimal solution set by retaining more non-dominated solutions. These solutions might be eliminated in advance within the process of service activities selection, when service activities selection and resource allocation are divided into two independent processes in the PSS configuration optimization. Furthermore, the above-mentioned existing studies highlight the multi-objective optimization of the overall solution. The proposed approach improves the executing reliabilities of service activities with a higher importance degree, by allocating resources with higher reliabilities to them. In terms of modeling, a multilayer network model can provide a better description of relationships among three types of elements (product parts, service activities, and resources), compared to an ontology-based model [43], structure tree model [44], and Petri net model [45,46]. The mathematical characteristics of multilayer networks are useful for analyzing the importance degree of product parts and service activities. The sparse matrixes of edges in multilayer networks are applicable as the input of the proposed optimization algorithm.
This paper focused on product-oriented PSS. Adding services to traditional products is used in product-oriented PSS as an approach to reach sustainability [66]. The ultimate purpose of the proposed approach is to achieve PSS configuration optimization. Even if the purpose of sustainability is not directly mentioned, the PSS configuration optimization is still related to achieving sustainability. Optimization of the executing time means a reduction in energy consumption and decreased work time of workers and equipment. The reduced resource utilization can lead to the concept of sustainability in PSS [5]. The optimization of executing costs is an economic benefit. The optimization of executing reliability can also be regarded as a guarantee of the optimization of executing time and cost. Furthermore, one goal of the proposed multi-objective optimization is to improve customer satisfaction. Economic perspective and customer satisfaction are also important measurement indicators of sustainability [9].

7. Conclusions

PSS enables manufacturing companies to provide customers with solutions by integrating preferable product modules and service modules. How to optimize PSS configuration in a fine-grained perspective when modules are confirmed represents an important issue for PSS management. A PSS multilayer network-based method combining service activities selection and resource allocation was developed to obtain optimal PSS solutions. Product parts, services, and resources were mapped as different layers in a PSS multilayer network. In particular, the mapping rules from logical relationships among service activities to edges among nodes were developed to model service layers appropriately. Service activities selection and resource allocation were combined to evaluate the objective functions of solutions systematically and precisely. The importance degree of a service activity was proposed to ascertain the priority of service activities in resource allocation. The advantages of the proposed method are as follows:
  • The impact of resources on the cost and reliability of service activities is thoroughly considered, which can affect the PSS configuration optimization equally with service activities. A combination of service activities selection and resource allocation can make the calculation of optimizing indexes closer to the actual implementation;
  • The PSS multilayer network model, which accommodates heterogeneous nodes and relationships without distortion, can be transformed into matrices to present information. It helps designers understand the construction of the PSS and build optimization functions easily;
  • Cross layer importance degree can sort the service activities based on the importance of the product parts and the objects of service activities.
Although the approach optimizes PSS configuration in a more fine-grained perspective, it also has some limitations. The combination of service activities selection and resource allocation leads to more complexity within the calculation. The optimization is based on a given product in a situation where customers have explicit targets on product module selection. Further work would consider improving the algorithm to shorten the calculation’s response time. Besides, the multilayer network-based modeling method might be extended to the design of use-oriented PSS and result-oriented PSS to satisfy demands for different types of PSS.

Author Contributions

Conceptualization, Z.Z. and D.X.; methodology, Z.Z. and D.X.; software, D.X.; validation, Z.Z., D.X., E.O. and H.C.; formal analysis, D.X.; investigation, D.X.; resources, Z.Z.; data curation, D.X.; writing—original draft, D.X.; writing—review and editing, Z.Z., D.X., E.O. and H.C.; supervision, Z.Z.; project administration, Z.Z.; funding acquisition, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51205242, and the Shanghai Science and Technology Innovation Action Plan, grant number 16111106402.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

min T S = j = 1 δ y j T j
min C S = j = 1 δ y j ( C j + i = 1 α + β x j i c i T j )
max A V R E S = 1 j = 1 δ y j D j j = 1 δ ( y j D j i = 1 α + β z j i )
y j ( φ = δ + 1 δ + γ e j φ S u φ ) = y j , j = 1 , 2 , , δ
j = 1 δ y j e φ j S w φ j S = u φ , φ = δ + 1 , ,   δ + γ
j = 1 δ y j e j 0 S w j 0 S = σ
[ ( j = 1 δ y j e i j S w i j S ) + e i 0 S w i 0 S ] y i = y i [ ( φ = δ + 1 δ + γ e φ i S w φ i S ) + ( j = 1 δ y j e j i S w j i S ) ] , i = 1 , 2 , , δ
y j e j φ S f t j u φ d d l _ t φ , j = 1 , 2 , , δ ; φ = δ + 1 , , δ + γ
i = 1 α + β e h i R x j i = y j e j h S R w j h S R , j = 1 , 2 , , δ ; h = α + β + 1 , , α + β + b
i = 1 α + β j = 1 δ h = 1 δ ( x j i x h i o l j h ) = 0 , j h
o l j h = { 1 , ( b t j , f t j ) ( b t h , f t h ) 0 , ( b t j , f t j ) ( b t h , f t h ) = , j h
z j i = { w i j R S , x j i = 1 1 , x j i = 0
x j i = a x j i e i j R S y j
a x j i ( a x j i 1 ) = 0
y j ( y j 1 ) = 0

References

  1. Meier, H.; Völker, O.; Funke, B. Industrial product-service systems (IPS 2). Int. J. Adv. Manuf. Technol. 2011, 52, 1175–1191. [Google Scholar] [CrossRef]
  2. Lindahl, M.; Sundin, E.; Sakao, T. Environmental and economic benefits of Integrated Product Service Offerings quantified with real business cases. J. Clean. Prod. 2014, 64, 288–296. [Google Scholar] [CrossRef] [Green Version]
  3. Mont, O.K. Clarifying the concept of product–service system. J. Clean. Prod. 2002, 10, 237–245. [Google Scholar] [CrossRef]
  4. Tukker, A.; Tischner, U. Product-services as a research field: Past, present and future. Reflections from a decade of research. J. Clean. Prod. 2006, 14, 1552–1556. [Google Scholar] [CrossRef]
  5. Baines, T.S.; Lightfoot, H.W.; Evans, S.; Neely, A.; Greenough, R.; Peppard, J.; Roy, R.; Shehab, E.; Braganza, A.; Tiwari, A. State-of-the-art in product-service systems. Proc. Inst. Mech. Eng. Part B 2007, 221, 1543–1552. [Google Scholar] [CrossRef] [Green Version]
  6. Tukker, A. Product services for a resource-efficient and circular economy—A review. J. Clean. Prod. 2015, 97, 76–91. [Google Scholar] [CrossRef]
  7. Rapaccini, M. Pricing strategies of service offerings in manufacturing companies: A literature review and empirical investigation. Prod. Plan. Control. 2015, 26, 1247–1263. [Google Scholar] [CrossRef]
  8. Tukker, A. Eight types of product–service system: Eight ways to sustainability? Experiences from SusProNet. Bus. Strategy Environ. 2004, 13, 246–260. [Google Scholar] [CrossRef]
  9. Lee, S.; Geum, Y.; Lee, H.; Park, Y. Dynamic and multidimensional measurement of product-service system (PSS) sustainability: A triple bottom line (TBL)-based system dynamics approach. J. Clean. Prod. 2012, 32, 173–182. [Google Scholar] [CrossRef]
  10. Sakao, T. Engineering PSS (Product/Service Systems) toward sustainability: Review of research. In Handbook of Sustainable Engineering; Kauffman, J., Lee, K., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 597–613. [Google Scholar]
  11. Pascual, R.; Siña, M.; Santelices, G.; Román, M.; López Droguett, E. Optimal channel coordination in use-based product-service system contracts. Int. J. Prod. Res. 2017, 55, 6946–6956. [Google Scholar] [CrossRef]
  12. Benedettini, O.; Neely, A.; Swink, M. Why do servitized firms fail? A risk-based explanation. Int. J. Oper. Prod. Manag. 2015, 35, 946–979. [Google Scholar] [CrossRef]
  13. Gebauer, H.; Fleisch, E. An investigation of the relationship between behavioral processes, motivation, investments in the service business and service revenue. Ind. Mark. Manag. 2007, 36, 337–348. [Google Scholar] [CrossRef]
  14. Neely, A. Exploring the financial consequences of the servitization of manufacturing. Oper. Manag. Res. 2008, 1, 103–118. [Google Scholar] [CrossRef] [Green Version]
  15. Cavalieri, S.; Pezzotta, G. Product–Service Systems Engineering: State of the art and research challenges. Comput. Ind. 2012, 63, 278–288. [Google Scholar] [CrossRef]
  16. Matschewsky, J.; Kambanou, M.L.; Sakao, T. Designing and providing integrated product-service systems–challenges, opportunities and solutions resulting from prescriptive approaches in two industrial companies. Int. J. Prod. Res. 2018, 56, 2150–2168. [Google Scholar] [CrossRef] [Green Version]
  17. Hussain, R.; Lockett, H.; Vasantha, G.V.A. A framework to inform PSS Conceptual Design by using system-in-use data. Comput. Ind. 2012, 63, 319–327. [Google Scholar] [CrossRef] [Green Version]
  18. Song, W.; Sakao, T. A customization-oriented framework for design of sustainable product/service system. J. Clean. Prod. 2017, 140, 1672–1685. [Google Scholar] [CrossRef]
  19. Mourtzis, D.; Fotia, S.; Vlachou, E.; Koutoupes, A. A Lean PSS design and evaluation framework supported by KPI monitoring and context sensitivity tools. Int. J. Adv. Manuf. Technol. 2018, 94, 1623–1637. [Google Scholar] [CrossRef]
  20. Marilungo, E.; Peruzzini, M.; Germani, M. Review of Product-Service System Design Methods. In Product Lifecycle Management in the Era of Internet of Things, Proceedings of the IFIP International Conference on Product Lifecycle Management, Doha, Qatar, 19–21 October 2015; Springer: Cham, Switzerland, 2016; pp. 271–279. [Google Scholar]
  21. Geng, X.; Xu, S.; Ye, C. Module configuration approach for product service system design driven by customer requirements. Comput. Integr. Manuf. Syst. 2016, 22, 55–61. [Google Scholar]
  22. Sun, J.; Chai, N.; Pi, G.; Zhang, Z.; Fan, B. Modularization of product service system based on functional requirement. Procedia CIRP 2017, 64, 301–305. [Google Scholar] [CrossRef]
  23. Zhang, Z.; Shang, Y.; Sun, J.; Fan, B. Product services system configuration optimization based on multi-objective discrete cuckoo search algorithm. Comput. Integr. Manuf. Syst. 2017, 23, 1774–1786. [Google Scholar]
  24. Shen, J.; Erkoyuncu, J.A.; Roy, R.; Wu, B. A framework for cost evaluation in product service system configuration. Int. J. Prod. Res. 2017, 55, 6120–6144. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, Z.; Xu, D.; Ostrosi, E.; Yu, L.; Fan, B. A systematic decision-making method for evaluating design alternatives of product service system based on variable precision rough set. J. Intell. Manuf. 2019, 30, 1895–1909. [Google Scholar] [CrossRef]
  26. Li, H.; Ji, Y.; Gu, X.; Qi, G.; Tang, R. Module partition process model and method of integrated service product. Comput. Ind. 2012, 63, 298–308. [Google Scholar] [CrossRef]
  27. Long, H.; Wang, L.; Shen, J.; Wu, M.; Jiang, Z. Product service system configuration based on support vector machine considering customer perception. Int. J. Prod. Res. 2013, 51, 5450–5468. [Google Scholar] [CrossRef]
  28. Wu, B.; Fan, S.; Yu, A.J.; Xi, L. Configuration and operation architecture for dynamic cellular manufacturing product–service system. J. Clean. Prod. 2016, 131, 716–727. [Google Scholar] [CrossRef]
  29. Long, H.J.; Wang, L.Y. A rough set based approach to knowledge acquisition for product service system configuration. Appl. Mech. Mater. 2012, 220–223, 2534–2539. [Google Scholar] [CrossRef]
  30. Zhang, Z.; Chai, N.; Ostrosi, E.; Shang, Y. Extraction of association rules in the schematic design of product service system based on Pareto-MODGDFA. Comput. Ind. Eng. 2019, 129, 392–403. [Google Scholar] [CrossRef]
  31. Liu, Y.; Zhang, Z.F.; Yao, D.; Chu, X.N. Extraction of Product Service System Configuration Rules Based on Multi-objective DPSO Algorithm. J. Shanghai Jiaotong Univ. 2015, 49, 1123–1130. [Google Scholar]
  32. Shen, J.; Wu, B.; Yu, L. Personalized configuration rules extraction in product service systems by using Local Cluster Neural Network. Ind. Manag. Data Syst. 2015, 115, 1529–1546. [Google Scholar] [CrossRef]
  33. Long, H.; Wang, L.; Zhao, S.; Jiang, Z. An approach to rule extraction for product service system configuration that considers customer perception. Int. J. Prod. Res. 2016, 54, 5337–5360. [Google Scholar] [CrossRef]
  34. Li, H.; Ji, Y.; Chen, L.; Jiao, R.J. Bi-level coordinated configuration optimization for product-service system modular design. IEEE Trans. Syst. Man Cybern. Syst. 2015, 47, 537–554. [Google Scholar] [CrossRef]
  35. Yan, H.; Zhong, W. Modular decomposition model of complex product system and its application. J. Beijing Univ. Aeronaut. Astronaut. 2017, 43, 654–659. [Google Scholar]
  36. Pezzotta, G.; Pirola, F.; Pinto, R.; Akasaka, F.; Shimomura, Y. A Service Engineering framework to design and assess an integrated product-service. Mechatronics 2015, 31, 169–179. [Google Scholar] [CrossRef]
  37. Shimomura, Y.; Nemoto, Y.; Ishii, T.; Nakamura, T. A method for identifying customer orientations and requirements for product–service systems design. Int. J. Prod. Res. 2018, 56, 2585–2595. [Google Scholar] [CrossRef]
  38. Szwejczewski, M.; Goffin, K.; Anagnostopoulos, Z. Product service systems, after-sales service and new product development. Int. J. Prod. Res. 2015, 53, 5334–5353. [Google Scholar] [CrossRef]
  39. Fargnoli, M.; Haber, N.; Sakao, T. PSS modularisation: A customer-driven integrated approach. Int. J. Prod. Res. 2018, 57, 1–17. [Google Scholar] [CrossRef] [Green Version]
  40. Mourtzis, D.; Fotia, S.; Boli, N.; Vlachou, E. An approach for the modelling and quantification of PSS customisation. Int. J. Prod. Res. 2018, 56, 1137–1153. [Google Scholar] [CrossRef]
  41. Pezzotta, G.; Pirola, F.; Rondini, A.; Pinto, R.; Ouertani, M.-Z. Towards a methodology to engineer industrial product-service system–Evidence from power and automation industry. CIRP J. Manuf. Sci. Technol. 2016, 15, 19–32. [Google Scholar] [CrossRef]
  42. Zhang, M.-W.; Zhang, B.; Liu, Y.; Na, J.; Zhu, Z.-L. Web service composition based on QoS rules. J. Comput. Sci. Technol. 2010, 25, 1143–1156. [Google Scholar] [CrossRef]
  43. Zhu, H.; Gao, J.; Li, D.; Tang, D. A Web-based Product Service System for aerospace maintenance, repair and overhaul services. Comput. Ind. 2012, 63, 338–348. [Google Scholar] [CrossRef]
  44. Cao, J.; Wang, J.; Zhao, H.; Sun, X. A service process optimization method based on model refinement. J. Supercomput. 2013, 63, 72–88. [Google Scholar] [CrossRef]
  45. Wang, Y.; Lin, C.; Ungsunan, P.D.; Huang, X. Modeling and survivability analysis of service composition using Stochastic Petri Nets. J. Supercomput. 2011, 56, 79–105. [Google Scholar] [CrossRef]
  46. Hu, Q.; Du, Y.; Yu, S. Service net algebra based on logic Petri nets. Inf. Sci. 2014, 268, 271–289. [Google Scholar] [CrossRef]
  47. Johnson, J. Hypernetworks in the Science of Complex Systems; Imperial College Press: London, UK, 2013. [Google Scholar]
  48. Johnson, J.H. Hypernetworks: Multidimensional relationships in multilevel systems. Eur. Phys. J. Spec. Top. 2016, 225, 1037–1052. [Google Scholar] [CrossRef] [Green Version]
  49. Johnson, J.H. Hypernetworks for reconstructing the dynamics of multilevel systems. In Proceedings of the European Conference on Complex Systems, Oxford, UK, 22–29 September 2006. [Google Scholar]
  50. Kivelä, M.; Arenas, A.; Barthelemy, M.; Gleeson, J.P.; Moreno, Y.; Porter, M.A. Multilayer networks. J. Complex Netw. 2014, 2, 203–271. [Google Scholar] [CrossRef] [Green Version]
  51. De Domenico, M.; Solé-Ribalta, A.; Cozzo, E.; Kivelä, M.; Moreno, Y.; Porter, M.A.; Gómez, S.; Arenas, A. Mathematical formulation of multilayer networks. Phys. Rev. X 2013, 3, 041022. [Google Scholar] [CrossRef] [Green Version]
  52. De Domenico, M.; Granell, C.; Porter, M.A.; Arenas, A. The physics of spreading processes in multilayer networks. Nat. Phys. 2016, 12, 901–906. [Google Scholar] [CrossRef] [Green Version]
  53. Omodei, E.; De Domenico, M.; Arenas, A. Evaluating the impact of interdisciplinary research: A multilayer network approach. Netw. Sci. 2017, 5, 235–246. [Google Scholar] [CrossRef] [Green Version]
  54. Pasqual, M.C.; de Weck, O.L. Multilayer network model for analysis and management of change propagation. Res. Eng. Des. 2012, 23, 305–328. [Google Scholar] [CrossRef]
  55. Leng, J.; Jiang, P. Dynamic scheduling in RFID-driven discrete manufacturing system by using multi-layer network metrics as heuristic information. J. Intell. Manuf. 2019, 30, 979–994. [Google Scholar] [CrossRef]
  56. Eppinger, S.D.; Browning, T.R. Design Structure Matrix Methods and Applications; MIT Press: Cambridge, MA, USA, 2012. [Google Scholar]
  57. Tomiyama, T. Service engineering to intensify service contents in product life cycles. In Proceedings of the Second International Symposium on Environmentally Conscious Design and Inverse Manufacturing, Tokyo, Japan, 11–15 December 2001; pp. 613–618. [Google Scholar]
  58. Cheng, Y.; Zhao, D.; Tao, F.; Zhang, L.; Liu, Y. Complex networks based manufacturing service and task management in cloud environment. In Proceedings of the 2015 IEEE 10th Conference on Industrial Electronics and Applications, Auckland, New Zealand, 15–17 June 2015; pp. 242–247. [Google Scholar]
  59. Chen, D.; Lü, L.; Shang, M.-S.; Zhang, Y.-C.; Zhou, T. Identifying influential nodes in complex networks. Phys. A 2012, 391, 1777–1787. [Google Scholar] [CrossRef] [Green Version]
  60. Lü, L.; Zhang, Y.-C.; Yeung, C.H.; Zhou, T. Leaders in social networks, the delicious case. PLoS ONE 2011, 6, e21202. [Google Scholar] [CrossRef] [Green Version]
  61. Gu, Y.; Zhu, Z. Node Ranking in Complex Networks based on LeaderRank and Modes Similarity. J. Univ. Electron. Sci. Technol. China 2017, 46, 441–448. [Google Scholar]
  62. Brin, S.; Page, L. The anatomy of a large-scale hypertextual web search engine. Comput. Netw. ISDN Syst. 1998, 30, 107–117. [Google Scholar] [CrossRef]
  63. Restrepo, J.G.; Ott, E.; Hunt, B.R. Characterizing the dynamical importance of network nodes and links. Phys. Rev. Lett. 2006, 97, 094102. [Google Scholar] [CrossRef] [Green Version]
  64. Li, Q.; Zhou, T.; Lü, L.; Chen, D. Identifying influential spreaders by weighted LeaderRank. Phys. A 2014, 404, 47–55. [Google Scholar] [CrossRef] [Green Version]
  65. Dudas, C.; Ng, A.H.; Pehrsson, L.; Boström, H. Integration of data mining and multi-objective optimisation for decision support in production systems development. Int. J. Comput. Integr. Manuf. 2014, 27, 824–839. [Google Scholar] [CrossRef]
  66. Manzini, E.; Vezzoli, C.; Clark, G. Product-service systems: Using an existing concept as a new approach to sustainability. J. Des. Res. 2001, 1, 27–40. [Google Scholar] [CrossRef]
Figure 1. Framework of the proposed method.
Figure 1. Framework of the proposed method.
Sustainability 12 00746 g001
Figure 2. Logical relationships among service activity nodes. (a) Serial relationship; (b) selective relationship; (c) parallel relationship.
Figure 2. Logical relationships among service activity nodes. (a) Serial relationship; (b) selective relationship; (c) parallel relationship.
Sustainability 12 00746 g002
Figure 3. Different connection situations of s a h . (a) Multiple inputs and multiple outputs; (b) Multiple inputs and single output; (c) Single input and multiple outputs.
Figure 3. Different connection situations of s a h . (a) Multiple inputs and multiple outputs; (b) Multiple inputs and single output; (c) Single input and multiple outputs.
Sustainability 12 00746 g003
Figure 4. The logical relationship between s a i and s a j . (a) It cannot be distinguished; (b) It is selective relationship; (c) It is parallel relationship.
Figure 4. The logical relationship between s a i and s a j . (a) It cannot be distinguished; (b) It is selective relationship; (c) It is parallel relationship.
Sustainability 12 00746 g004
Figure 5. Procedures of the improved Pareto-improved multiple-objective discrete cuckoo search algorithm (Pareto-IMODCS) algorithm.
Figure 5. Procedures of the improved Pareto-improved multiple-objective discrete cuckoo search algorithm (Pareto-IMODCS) algorithm.
Sustainability 12 00746 g005
Figure 6. Network status at time t.
Figure 6. Network status at time t.
Sustainability 12 00746 g006
Figure 7. The sparse matrixes of edges at time t . (a) Sparse matrixes of edges in N e t P ; (b) sparse matrixes of edges in N e t S ; (c) sparse matrixes of edges in N e t R ; (d) sparse matrixes of inter-layer edges from N e t S to N e t P ; (e) sparse matrixes of inter-layer edges from N e t S to N e t R ; (f) sparse matrixes of inter-layer edges from N e t R to N e t S .
Figure 7. The sparse matrixes of edges at time t . (a) Sparse matrixes of edges in N e t P ; (b) sparse matrixes of edges in N e t S ; (c) sparse matrixes of edges in N e t R ; (d) sparse matrixes of inter-layer edges from N e t S to N e t P ; (e) sparse matrixes of inter-layer edges from N e t S to N e t R ; (f) sparse matrixes of inter-layer edges from N e t R to N e t S .
Sustainability 12 00746 g007aSustainability 12 00746 g007b
Figure 8. Triggered service flows in the service layer at time t.
Figure 8. Triggered service flows in the service layer at time t.
Sustainability 12 00746 g008
Figure 9. Cost for all resources.
Figure 9. Cost for all resources.
Sustainability 12 00746 g009
Figure 10. e x r e 9 in different solutions in set A and set C.
Figure 10. e x r e 9 in different solutions in set A and set C.
Sustainability 12 00746 g010
Figure 11. Gantt chart of selected service activity nodes.
Figure 11. Gantt chart of selected service activity nodes.
Sustainability 12 00746 g011
Figure 12. w i j R S on edges that connect p r i to s a 8 , s a 9 , s a 19 , and s a 22 .
Figure 12. w i j R S on edges that connect p r i to s a 8 , s a 9 , s a 19 , and s a 22 .
Sustainability 12 00746 g012
Figure 13. w i j R S on edges that connect e r i to all s a .
Figure 13. w i j R S on edges that connect e r i to all s a .
Sustainability 12 00746 g013
Table 1. The executing time, cost, and reliability of related service activities.
Table 1. The executing time, cost, and reliability of related service activities.
s a i T i   ( h ) C i   ( 10 3   yuan ) R E i s a i T i   ( h ) C i   ( 10 3   yuan ) R E i
10 0 1.0000 132.6 3.10 0.9897
24.5 1.93 0.9721 141.3 1.61 0.9816
31.9 4.24 0.9984 152.4 2.17 0.9949
43.1 2.57 0.9811 161.6 2.31 0.9793
50 0 1.0000 172.0 1.98 0.9996
62.8 3.69 0.9846 183.1 3.73 0.9912
73.5 2.21 0.9930 191.8 4.19 0.9838
82.4 0.98 0.9840 202.8 1.91 0.9704
92.3 2.47 0.9949 211.3 3.56 0.9880
103.1 0.99 0.9974 222.3 2.87 0.9779
111.4 2.31 0.9859
122.9 1.81 0.9881
Table 2. Solutions in set A.
Table 2. Solutions in set A.
No.1231314
s a i
( r j )
11111
2(6,10,52,73,102)3(2,10,51,75,105)3(2,9,52,75,105)3(5,8,51,72,102)3(1,10,51,73,102)
4(18,25,42,85,90)4(19,25,42,85,90)4(21,25,42,85,90)4(21,25,42,85,90)4(19,29,42,83,90)
55555
6(4,15,35,72,93)6(9,15,39,73,94)6(4,16,39,73,95)6(10,16,35,74,95)6(2,14,35,75,94)
7(7,39,51,94,103)7(7,36,52,95,102)7(6,36,51,93,102)7(7,36,52,93,103)7(7,36,52,95,103)
8(16,28,38,43,77,88)8(14,31,36,47,80,87)8(14,31,37,46,79,86)8(16,29,38,44,79,86)8(16,31,34,44,80,88)
9(28,47,90,98)9(31,47,90,98)9(28,46,90,98)9(28,47,90,98)9(31,47,90,98)
10(2,38,74,95)10(1,38,72,92)10(7,38,72,94)12(13,29,38,94)12(13,25,37,91)
13(2,15,78)13(2,16,78)13(7,16,77)18(2,18,19,37,39,80,100)18(2,18,21,38,39,77,100)
14(12,37,79)14(15,38,77)14(13,38,78)21(31,52,89,102)21(32,52,88,103)
19(26,29,44,89)19(28,29,40,88)19(26,29,44,87)22(29,42,45,99)22(29,41,44,100)
22(29,43,44,99)22(27,44,46,100)22(29,44,47,100)
T S 29.727.127.125.625.6
C S 63.08 61.55 61.77 63.13 62.72
A V R E S 0.9760 0.9724 0.9746 0.9719 0.9699
Table 3. Solutions in set B.
Table 3. Solutions in set B.
No.1256
s a i
( r j )
1111
3(1,10,51,73,102)3(5,8,51,72,102)2(6,7,52,73,102)2(6,7,52,75,103)
4(19,29,42,83,90)4(21,25,42,85,90)4(19,25,42,83,90)4(21,25,42,83,90)
5555
6(2,14,35,75,94)6(10,16,35,74,95)6(8,16,35,75,95)6(10,16,35,72,95)
7(7,36,52,95,103)7(7,36,52,93,103)7(10,39,51,92,103)7(4,37,51,93,105)
8(16,31,34,44,80,88)8(16,29,38,44,79,86)8(16,28,39,43,78,86)8(16,28,39,44,78,86)
9(31,47,90,98)9(28,47,90,98)9(31,47,90,98)9(28,46,90,98)
12(13,25,37,91)12(13,29,38,94)12(12,28,36,93)12(15,29,38,91)
18(2,18,21,38,39,77,100)18(2,18,19,37,39,80,100)18(4,18,23,37,38,80,98)18(9,18,23,36,38,80,98)
21(32,52,88,103)21(31,52,89,102)21(32,52,89,103)21(29,51,88,102)
22(29,41,44,100)22(29,42,45,99)22(29,43,45,100)22(32,45,47,100)
T S 25.625.628.228.2
C S 62.72 63.13 65.64 64.95
A V R E S 0.9699 0.9719 0.9732 0.9721
Table 4. Resource quantity demands of s a
Table 4. Resource quantity demands of s a
k r 122 k r 124 k r 125 k r 126 k r 132 k r 134 k r 136
s a 9 0101011
s a 8 1111110
s a 19 0201010
s a 22 0102001
Table 5. r i that are allocated to s a 8 , s a 9 , s a 19 , and s a 22 .
Table 5. r i that are allocated to s a 8 , s a 9 , s a 19 , and s a 22 .
s a 8 s a 9 s a 19 s a 22
k r 124 k r 126 k r 134 k r 124 k r 126 k r 134 k r 136 k r 124 k r 126 k r 134 k r 124 k r 126 k r 136
A3 p r 31 p r 46 e r 86 p r 28 p r 46 e r 90 e r 98 p r 26 p r 29 p r 44 e r 87 p r 29 p r 44 p r 47 e r 100
C11 p r 30 p r 43 e r 90 p r 31 p r 47 e r 90 e r 97 p r 27 p r 32 p r 44 e r 86 p r 25 p r 44 p r 46 e r 98

Share and Cite

MDPI and ACS Style

Zhang, Z.; Xu, D.; Ostrosi, E.; Cheng, H. Optimization of the Product–Service System Configuration Based on a Multilayer Network. Sustainability 2020, 12, 746. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020746

AMA Style

Zhang Z, Xu D, Ostrosi E, Cheng H. Optimization of the Product–Service System Configuration Based on a Multilayer Network. Sustainability. 2020; 12(2):746. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020746

Chicago/Turabian Style

Zhang, Zaifang, Darao Xu, Egon Ostrosi, and Hui Cheng. 2020. "Optimization of the Product–Service System Configuration Based on a Multilayer Network" Sustainability 12, no. 2: 746. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020746

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