Product service system module configuration optimization method under participation of multiple benefit related parties

Abstract

The invention discloses a product service system module configuration optimization method under the participation of multiple benefit related parties, and belongs to the technical field of product design, and the method comprises the steps: constructing a product service system module configuration optimization model as an upper layer model; constructing a cost optimization model of benefit-related parties as a lower-layer model, wherein the benefit-related parties comprise suppliers, manufacturers and distributors; and constructing a double-layer model by using the constructed upper-layer model and the constructed lower-layer model, and solving the double-layer model based on a chaos adaptive nested genetic algorithm. According to different optimization objectives and constraint conditions of PSS design and benefit-related parties, a one-master multi-slave bilevel programming model for the benefit-related parties to participate in PSS design is established, and a PSS module configuration optimization method under the participation of multiple benefit-related parties is provided.

Description

Technical Field

[0001] The present invention relates to the technical field of product design, and in particular to configuration optimization of a product service system with the participation of multiple stakeholders, and specifically to a product service system module configuration optimization method with the participation of multiple stakeholders. Background Art

[0002] The design of Product Service System (PSS) should not only meet the diverse and personalized needs of stakeholders, but also consider the industrial efficiency of large-scale production, which requires the module configuration optimization of PSS. Adopting appropriate management and manufacturing methods for different categories of modules and optimizing configuration in combination with module constraints can effectively resolve the conflict between personalization and scale benefits and achieve large-scale customization. PSS module configuration is to build a series of configuration schemes by selecting different module candidates based on a set of constraints according to the needs of different stakeholders and the importance of technical characteristics. The configuration scheme is crucial to meet the needs of various stakeholders and is a basic indicator for measuring market satisfaction.

[0003] PSS integrates product modules and service modules into a system. Product modules and service modules are composed of several composite modules respectively, and composite modules are further assembled from basic modules. Basic modules are the basic units that constitute PSS. They include mandatory modules and optional modules. Mandatory modules are necessary to constitute products, while optional modules provide additional functions or features and can be selected as needed. Each basic module has several module candidates, which provide different design, performance or cost options, bringing differentiation and flexibility to the product. In the configuration process of PSS, the selection of module candidates, the assembly method of composite modules, and the selection of manufacturers and distributors are all affected by stakeholders. Different suppliers may provide module candidates with different quality, price and delivery time, while different manufacturers may have different manufacturing capabilities and technical advantages. These factors will affect the cost, quality and delivery time of PSS. Therefore, when configuring PSS, it is necessary to fully consider the opinions and needs of stakeholders to ensure that the final product service system can meet the expectations of the market and customers. In addition, due to production capacity limitations, not all possible module combinations will be considered for production. This requires trade-offs and optimization during the module configuration process to select the combination with the greatest economic benefits and market potential.

[0004] In summary, considering the close interaction between physical products and services, services and services in the system, and the complexity of PSS, it is necessary to establish a systematic module configuration optimization plan, combining the importance of technical characteristics, as well as module relevance, mutual exclusivity and other constraints, to optimize the module configuration of services and physical products under the premise of ensuring compliance with the module division principle, and finally obtain a PSS configuration plan with the participation of multiple stakeholders. In order to simultaneously ensure the configuration constraints of PSS itself and the effective participation of stakeholders, improve production and service efficiency, and reduce resource consumption. Summary of the invention

[0005] To solve the above problems, the main purpose of the present invention is to provide a product service system module configuration optimization method with the participation of multiple stakeholders.

[0006] The technical solution of the present invention is as follows:

[0007] A method for optimizing product service system module configuration with multi-stakeholder participation includes the following steps:

[0008] S1. Construct a product service system module configuration optimization model as the upper model;

[0009] The upper-level model makes decisions on composite modules and the composition of modules;

[0010] The product configuration scheme is as follows:

[0011] From the module candidate set Select different module candidates to form different basic modules m j , all basic modules constitute the basic module set {m1,…,m j ,…,m J}; Select different basic modules from the basic module set to form different composite modules M k , all composite modules form a composite module set {M1,…,M k ,…,M K}; Select different composite modules from the composite module set to form different product solutions PSS i , all product solutions constitute the product solution set PSS;

[0012] Basic modules are divided into two categories, optional modules and mandatory modules. The set of all optional modules is the optional module set May, and the set of all mandatory modules is the mandatory module set Must.

[0013] Represent the product solution in vector form, X i ={x ikjl}, X i represents the i-th product configuration scheme (PSS i)’s decision vector; x ikjl represents the selection variable of the kth composite module, jth basic module, and lth module candidate in the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, x ikjl = 0, when the product configuration includes the selected variable, x ikjl =1;

[0014] The upper-level objective function is the basis for product configuration decisions. It makes decisions on the number and composition of configurable composite modules, and selects known required and optional modules and their corresponding module candidates to obtain the module structure of the configuration solution. In the configuration process, multiple factors need to be weighed to make the best decision. The upper-level objective function is used to integrate these factors to form a clear decision-making guidance standard. Maximizing customer and corporate interests is the most basic criterion in PSS configuration decisions, which is characterized by the utility of PSS as perceived by customers.

[0015]

[0016] Where U i Represents the i-th product configuration scheme PSS i The utility of jl is the utility of the lth module candidate of the jth basic module as considered by the customer; ω ij The technical characteristics corresponding to the jth basic module in the product PSS i The weight in; K represents the total number of composite modules; J represents the total number of basic modules; L j Represents the total number of candidate modules for the jth basic module;

[0017] The interests of other stakeholders are also a multi-dimensional comprehensive indicator. For suppliers, manufacturers and distributors, profit is the core embodiment of their interests, and it is also the most important and objective reflection, which directly reflects the economic benefits and operational efficiency of the enterprise. Therefore, the cost of stakeholders is used as a representation of their respective interests.

[0018] The cost function of stakeholders is:

[0019]

[0020] In the two-level optimization model, the minimum total cost of the lower-level suppliers, manufacturers, and distributors (respectively expressed as ) is the key decision variable. The profit cost rate of the upper-level architecture for the lower-level architecture is recorded as α, α∈(0,1). The lower level of the model is the configuration plan to optimize the selection of suppliers, manufacturers, and distributors in the supply, manufacturing, and distribution links, pursuing the overall cost optimization of the three.

[0021] The maximum value of the ratio of customer utility to stakeholder cost is taken as the objective function, that is, the utility per unit cost is maximized, as shown in the following formula:

[0022]

[0023] The constraints are:

[0024]

[0025] x ikjl ∈{0,1}

[0026] α∈(0,1)

[0027] i=1,…,I; k=1,…,K; j=1,…,J; l=1,…,L j

[0028] I,K∈N +

[0029] In the formula, m j represents the jth basic module; Must represents the set of mandatory modules; May represents the set of optional modules; I represents the total number of configuration schemes; h represents the total number of mandatory modules; a represents the total number of optional modules; N + represents a positive integer;

[0030] S2, construct the stakeholder cost optimization model as the lower model;

[0031] The lower-level functions include manufacturers, suppliers, and distributors, stakeholders in the supply chain. The model costs compete with each other in the optimization process, and the maximum value of the utility of the upper-level customers is finally achieved to output the product service system configuration plan.

[0032] S21. The supplier cost function should comprehensively consider the raw material cost, ordering cost, inventory cost and transportation cost corresponding to each module candidate, and its corresponding objective function is as follows:

[0033]

[0034] The corresponding constraints are as follows:

[0035]

[0036]

[0037] In the formula, represents the material cost of the lth module candidate of the jth basic module in the sth supplier; d jl PSS iThe demand for the lth module candidate of the jth basic module; represents the unit ordering cost of the lth module candidate of the jth basic module from the sth supplier; represents the unit inventory cost of the lth module candidate of the jth basic module in the sth supplier; represents the unit transportation cost of the lth module candidate of the jth basic module in the sth supplier; represents the allocation function of the selected suppliers on the supply amount; Q is the total number of suppliers who can provide the lth module candidate of the jth basic module; j PSS i The demand for the jth basic module;

[0038] For the upper model, PSS i Select the corresponding basic module m j ,Right now Then any supplier s in the supplier set S corresponding to the module needs to make a decision for this situation, deciding whether to supply this basic module candidate. The decision variable is recorded as When the sth supplier provides the jth basic module and the lth module candidate, then otherwise

[0039] In the actual design and production process of PSS, there may be compatibility and incompatibility constraints between modules. For incompatibility constraints, when two different optional modules cannot appear in the same product at the same time, it constitutes module incompatibility. If two modules cannot be selected at the same time, the constraints are as follows:

[0040]

[0041] In the formula, The selection variable representing the lth module candidate of the j1th basic module of the kth composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, The selection variable representing the lth module candidate of the j2th basic module of the k'th composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable,

[0042] As for compatibility constraints, especially when two optional modules must appear at the same time, this usually means that the two modules are interdependent, or they together constitute a basic functional unit of the product. If two modules need to be selected at the same time, the module compatibility constraints are as follows:

[0043]

[0044] S22. The manufacturer's cost function should comprehensively consider the fixed costs corresponding to each composite module, the variable costs such as mechanical equipment and labor input. The corresponding objective function is as follows:

[0045]

[0046] The corresponding constraints are as follows:

[0047]

[0048] In the formula, represents the operating fixed cost of the rth manufacturer of the kth composite module; R k represents the total number of manufacturers; represents the unit operating variable cost of the rth manufacturer to assemble and manufacture the lth module candidate of the jth basic module; d k represents the demand for the kth composite module; It represents the process variation cost of the rth manufacturer to assemble and manufacture the required module j as the module candidate l;

[0049] For the upper model, PSS i Select the corresponding basic module m j ,Right now Then any manufacturer r in the manufacturer set R corresponding to the module needs to make a decision for this situation, deciding whether to manufacture the required or optional module candidate. The decision variable is recorded as When the rth manufacturer provides the kth composite module, then otherwise

[0050] represents the proportion of the demand for the kth composite module produced by the rth manufacturer; Select the total number of manufacturers for the kth composite module; is the production capacity range of the kth composite module by the rth manufacturer;

[0051] S23. The distributor's cost function should comprehensively consider the ordering cost, inventory cost, and shipping cost of the product. The corresponding objective function is as follows:

[0052]

[0053] The corresponding constraints are as follows:

[0054]

[0055] y b ∈{0,1}

[0056] In the formula, for the upper model, the manufacturer selects the corresponding mandatory module m i or optional module m j To form a composite module and further assemble it into a product, any distributor b in the distributor set B corresponding to the product needs to make a decision for this situation, deciding whether to distribute the product assembled from this composite module. The decision variable is denoted as y b , when the bth distributor decides to distribute, then y b =1, otherwise y b =0;

[0057] N B Select the number of distributors for each product;

[0058] RC b represents the unit ordering cost of the b-th distributor;

[0059] HC b represents the unit inventory cost of the product in the bth distributor;

[0060] TCP b represents the unit transportation cost of the product in the bth distributor;

[0061] f(y b ) represents the proportion function of the b-th distributor in the distribution of products;

[0062] d represents the total demand for the product;

[0063] W b It represents the distribution capacity range of the product by the b-th distributor;

[0064] The upper model constructed by S3 and S1 and the lower model constructed by S2 constitute a double-layer model, which is solved based on the chaotic adaptive nested genetic algorithm;

[0065] The specific steps of the chaotic adaptive nested genetic algorithm process are as follows:

[0066] S31. Parameter initialization: Determine the size of the initial population M, the genetic generation G, the chaotic state control parameter μ, the elite retention number m, the upper and lower bounds p of the adaptive crossover and mutation operators c0 、p c1 、p m0 and pm1 ;

[0067] S32, chaos initialization population: performing chaos optimization operation on the initial populations of the upper and lower layer models respectively;

[0068] The calculation formula for determining the chaotic variables in the upper model is as follows:

[0069] x k+1 =μx k (1-x k )

[0070] k=0,1,…0≤x k ≤1

[0071] In the formula, x k is the result of the chaotic variable x after k iterations, μ is the chaotic state control parameter, μ∈[0,4], when μ=4, and When , the system is in a completely chaotic state, the chaotic variable traverses and searches in the range of (0,1), and the generated initial sequence value will have completely chaotic characteristics.

[0072] The calculation formula for determining the chaotic variables in the lower model is as follows:

[0073] t ik =Int(nμx k-1 (1-x k-1 ))i=0,1,…,m-1

[0074] t ik is the value of the i-th optimization variable after the k-th iteration in the mathematical model. Its value is the number of candidate items. If there are n+1 candidate items, the value range is: ; n is the number of module candidates corresponding to each basic module;

[0075] S33. Assume that the population size in the algorithm is M, sort all individuals in the population according to their fitness values ​​from large to small, select the first m individuals with the best fitness values ​​and directly pass them on to the next generation; for the remaining Mm individuals, use a mixed ranking strategy to select individuals participating in crossover and mutation from these individuals.

[0076] The calculation formula for the fitness value of the i-th individual is as follows:

[0077]

[0078] S34, combining the selection mechanism of elite retention and mixed ranking to screen out the parent individuals for subsequent crossover and mutation operations;

[0079] S35, adaptive crossover and mutation operation: calculate the corresponding crossover and mutation probabilities based on the fitness value of the individual, and complete the crossover and mutation operations;

[0080] The sire is determined by roulette selection, and the probability of being selected as the sire is determined according to the fitness value of the individual; the probability of the i-th individual being selected can be calculated as:

[0081]

[0082] Adaptive crossover probability p c , adaptive mutation probability p m The calculation formulas are:

[0083]

[0084] In the formula, p c0 、p c1 are the initial adjustment factors for the lower and upper bounds of the crossover probability, which together determine the fluctuation range of the crossover probability; p m0 、p m1 are the initial adjustment factors for the lower and upper bounds of the mutation probability, respectively, which limit the upper and lower limits of the mutation probability; f is defined as the individual fitness value; f max is the optimal fitness value of individuals in the current population; f min is the worst fitness value of individuals in the population; f avg is the average fitness value of individuals in the population;

[0085] S36. If the number of iterations reaches the upper limit of the genetic algebra G preset in step S31, the algorithm ends and outputs the optimal module configuration solution; otherwise, go to step S33.

[0086] Beneficial effects:

[0087] Aiming at the master-slave association optimization problem of PSS module configuration, the present invention establishes a one-master-multiple-slave double-layer configuration optimization model in which stakeholders participate in PSS design.

[0088] In order to effectively solve the model, a chaotic adaptive nested genetic algorithm was developed. The population was initialized by chaos theory, and the advantages of chaos algorithm and genetic algorithm were integrated to form complementarity and improve the algorithm's optimization ability, so as to obtain the optimal solution for PSS module configuration with the participation of multiple stakeholders, helping enterprises to design products with core competitiveness that meet the needs of stakeholders.

[0089] The main advantages of this method include the following:

[0090] (1) Aiming at the different optimization objectives and constraints of PSS design and stakeholders, the present invention establishes a one-master-multiple-slave two-level planning model for stakeholders to participate in PSS design, and proposes a PSS module configuration optimization method with the participation of multiple stakeholders.

[0091] (2) The present invention proposes a chaotic adaptive nested genetic algorithm to effectively solve the one-master-multiple-slave two-level programming model in which stakeholders participate in the design of PSS. The algorithm integrates the advantages of the chaotic algorithm and the genetic algorithm to form a complementary relationship, improves the algorithm's optimization ability, can avoid individuals from concentrating in a certain local area, and is more likely to contain individuals close to the global optimal solution, which helps to narrow the scope of subsequent searches; it can dynamically adjust its parameters according to the characteristics of the problem and the search process, and conduct a detailed search near the optimal solution, further improving the algorithm's convergence speed and global optimization ability; it can maintain the diversity of the group without increasing the population size, which helps to reduce randomness and avoid premature convergence, and improve the probability of the algorithm finding the global optimal solution. BRIEF DESCRIPTION OF THE DRAWINGS

[0092] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below.

[0093] Figure 1 1 is a fitness value comparison diagram of three composite module division methods in an embodiment of the present invention;

[0094] Figure 2 It is the fitness value, customer utility and stakeholder cost change curve of the optimal composite module partitioning method in the case of the present invention;

[0095] Figure 3 It is a schematic diagram of a one-master-multiple-slave two-layer planning model in which stakeholders of the present invention participate in the design of PSS;

[0096] Figure 4 It is a schematic diagram of the chaotic adaptive nested genetic algorithm of the present invention. DETAILED DESCRIPTION

[0097] The present invention will be further described in detail below in conjunction with the embodiments and drawings.

[0098] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. The described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0099] Example

[0100] A method for optimizing product service system module configuration with multi-stakeholder participation includes the following steps:

[0101] S1. Construct a product service system module configuration optimization model as the upper model;

[0102] The product configuration scheme in the upper model is as follows: different module candidates are selected from the module candidate set to form different basic modules, and all basic modules form a basic module set; different basic modules are selected from the basic module set to form different composite modules, and all composite modules form a composite module set; different composite modules are selected from the composite module set to form different product solutions, and all product solutions form a product solution set; basic modules are divided into two categories, namely optional modules and mandatory modules, the set composed of all optional modules is the optional module set, and the set composed of all mandatory modules is the mandatory module set;

[0103] The objective function of the upper model is:

[0104]

[0105] in,

[0106] In the formula, F(x ikjl ) represents the cost-utility function; U i Represents the i-th product configuration scheme PSS i The utility of TC i represents the cost of stakeholders in the i-th product configuration scheme; K represents the total number of composite modules; J represents the total number of basic modules; L j represents the total number of module candidates for the jth basic module; ω ij The technical characteristics corresponding to the jth basic module in the product PSS i The weight in μ jl x is the utility of the lth module candidate of the jth basic module that the customer considers; ikjl represents the selection variable of the lth module candidate of the jth basic module of the kth composite module of the ith product configuration scheme. When the product configuration scheme does not contain this selection variable, x ikjl =0, when the product configuration includes the selected variable, x ikjl =1; α is the profit cost rate of the upper structure based on the benefit analysis of the lower structure; denote the minimum total cost of the lower-level suppliers, manufacturers, and distributors, respectively;

[0107] The constraints of the upper model include:

[0108]

[0109] x ikjl ∈{0,1}

[0110] α∈(0,1)

[0111] i=1,…,I; k=1,…,K; j=1,…,J; l=1,…,L j

[0112] I,K∈N +

[0113] In the formula, m j represents the jth basic module; Must represents the set of mandatory modules; May represents the set of optional modules; h represents the total number of mandatory modules; a represents the total number of optional modules; I represents the total number of configuration schemes; N + represents a positive integer;

[0114] S2. Constructing a stakeholder cost optimization model as a lower-level model; the stakeholders include suppliers, manufacturers, and distributors;

[0115] The supplier's objective function in the lower model is:

[0116]

[0117] The constraints are:

[0118]

[0119] In the formula, S jl is the total number of suppliers who provide candidates, is a decision variable. When the sth supplier provides the lth module candidate of the jth basic module, otherwise represents the material cost of the lth module candidate of the jth basic module in the sth supplier; d jl PSS i The demand for the lth module candidate of the jth basic module; represents the unit ordering cost of the lth module candidate of the jth basic module from the sth supplier; represents the unit inventory cost of the lth module candidate of the jth basic module in the sth supplier; represents the unit transportation cost of the lth module candidate of the jth basic module in the sth supplier; represents the allocation function of the selected suppliers on the supply amount; is the total number of suppliers who can provide the lth module candidate of the jth basic module; Q j PSS i The demand for the jth basic module;

[0120] In addition, compatibility and incompatibility constraints are included;

[0121] The constraints for the two modules to be incompatible (cannot be selected at the same time) are:

[0122]

[0123] In the formula, The selection variable representing the lth module candidate of the j1th basic module of the kth composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, The selection variable representing the lth module candidate of the j2th basic module of the k'th composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable,

[0124] The constraints for two modules to be compatible (selected at the same time) are:

[0125]

[0126] In the formula, The selection variable representing the lth module candidate of the j3rd basic module of the kth composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, The selection variable representing the lth module candidate of the j4th basic module of the k'th composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable,

[0127] The objective function of the manufacturer in the lower model is:

[0128]

[0129] The constraints are:

[0130]

[0131] In the formula, R k represents the total number of manufacturers; represents the operating fixed cost of the rth manufacturer of the kth composite module; represents the unit operating variable cost of the rth manufacturer to assemble and manufacture the lth module candidate of the jth basic module; d k represents the demand for the kth composite module; It represents the process variation cost of the rth manufacturer to assemble and manufacture the required module j as the module candidate l; is a decision variable. When the rth manufacturer provides the kth composite module, otherwise represents the proportion of the demand for the kth composite module produced by the rth manufacturer; Select the total number of manufacturers for the kth composite module; is the production capacity range of the kth composite module by the rth manufacturer;

[0132] The objective function of the distributor in the lower model is:

[0133]

[0134] The constraints are:

[0135]

[0136] y b ∈{0,1}

[0137] Where B is the total number of distributors; N B Select the number of distributors for each product;

[0138] RC b represents the unit ordering cost of the b-th distributor;

[0139] HC b represents the unit inventory cost of the product in the bth distributor;

[0140] TCD b represents the unit transportation cost of the product in the bth distributor;

[0141] f(y b ) represents the proportion function of the b-th distributor in the distribution of products;

[0142] d represents the total demand for the product;

[0143] W b It represents the distribution capacity range of the product by the b-th distributor;

[0144] The upper model constructed by S3 and S1 and the lower model constructed by S2 constitute a double-layer model, which is solved based on the chaotic adaptive nested genetic algorithm;

[0145] The specific steps of the chaotic adaptive nested genetic algorithm process are as follows:

[0146] S31. Parameter initialization: Determine the size of the initial population M, the genetic generation G, the chaotic state control parameter μ, the elite retention number m, the upper and lower bounds p of the adaptive crossover and mutation operators c0 、p c1 、pm0 and p m1 ;

[0147] S32, chaos initialization population: performing chaos optimization operation on the initial populations of the upper and lower layer models respectively;

[0148] The calculation formula for determining the chaotic variables in the upper model is as follows:

[0149] x k+1 =μx k (1-x k )

[0150] k=0,1,…0≤x k ≤1

[0151] In the formula, x k is the result of the chaotic variable x after k iterations, μ is the chaotic state control parameter, μ∈[0,4], when μ=4, and When , the system is in a completely chaotic state, the chaotic variable traverses and searches in the range of (0,1), and the generated initial sequence value will have completely chaotic characteristics.

[0152] The calculation formula for determining the chaotic variables in the lower model is as follows:

[0153] t ik =Int(nμx k-1 (1-x k-1 ))i=0,1,…,m-1

[0154] t ik The value of the i-th optimization variable after the k-th iteration in the mathematical model is the number of module candidates. If there are n+1 module candidates, the value range is: 0, 1, ..., n; n is the number of module candidates corresponding to each basic module;

[0155] S33. The population size in the algorithm is M. All individuals in the population are sorted from large to small according to their fitness values, and the first m individuals with the best fitness values ​​are selected to be directly inherited to the next generation; for the remaining Mm individuals, a mixed ranking strategy is used to select individuals participating in crossover and mutation from these individuals.

[0156] The calculation formula for the fitness value of the i-th individual is as follows:

[0157]

[0158] S34, combining the selection mechanism of elite retention and mixed ranking to screen out the parent individuals for subsequent crossover and mutation operations;

[0159] S35, adaptive crossover and mutation operation: calculate the corresponding crossover and mutation probabilities based on the fitness value of the individual, and complete the crossover and mutation operations;

[0160] The sire is determined by roulette selection, and the probability of being selected as the sire is determined according to the fitness value of the individual; the probability of the i-th individual being selected can be calculated as:

[0161]

[0162] Adaptive crossover probability p c , adaptive mutation probability p m The calculation formulas are:

[0163]

[0164] In the formula, p c0 、p c1 are the initial adjustment factors for the lower and upper bounds of the crossover probability, which together determine the fluctuation range of the crossover probability; p m0 、p m1 are the initial adjustment factors for the lower and upper bounds of the mutation probability, respectively, which limit the upper and lower limits of the mutation probability; f is defined as the individual fitness value; f max is the optimal fitness value of individuals in the current population; f min is the worst fitness value of individuals in the population; f avg is the average fitness value of individuals in the population;

[0165] S36. If the number of iterations reaches the upper limit of the genetic algebra G preset in step S31, the algorithm ends and outputs the optimal module configuration solution; otherwise, go to step S33.

[0166] Beneficial effects:

[0167] Aiming at the master-slave association optimization problem of PSS module configuration, the present invention establishes a one-master-multiple-slave double-layer configuration optimization model in which stakeholders participate in PSS design.

[0168] In order to effectively solve the model, a chaotic adaptive nested genetic algorithm was developed. The population was initialized by chaos theory, and the advantages of chaos algorithm and genetic algorithm were integrated to form complementarity and improve the algorithm's optimization ability, so as to obtain the optimal solution for PSS module configuration with the participation of multiple stakeholders, helping enterprises to design products with core competitiveness that meet the needs of stakeholders.

[0169] In order to illustrate the practical application of the proposed method, the present invention takes the configuration optimization of the smart speaker product service system as an example, establishes a one-master-multiple-slave two-layer configuration optimization model in which stakeholders participate in the design of the product service system, solves the model through the designed improved CANGA algorithm, and outputs the optimal smart speaker product service system module configuration scheme.

[0170] The detailed steps are as follows:

[0171] Table 1 Classification of product service system modules

[0172]

[0173]

[0174] Three reasonable composite module division schemes are selected to facilitate subsequent analysis, as shown in Tables 2, 3, and 4:

[0175] Table 2 Division of the first composite module

[0176] Composite Modules Basic Modules Security Module <![CDATA[m4、m5、m7、m 10 ]]> Personalization module <![CDATA[m6、m8、m9]]> Power Module <![CDATA[m1、m2、m3]]>

[0177] Table 3 Division of the second composite module

[0178] Composite Modules Basic Modules Structural Modules <![CDATA[m5、m9]]> Communication sensor module <![CDATA[m4、m6、m7、m8、m 10 ]]> Power Module <![CDATA[m1、m2、m3]]>

[0179] Table 4 Division of the third composite module

[0180] Composite Modules Basic Modules Physical Product Module <![CDATA[m4、m5、m6、m7]]> Service Module <![CDATA[m8、m9、m 10 ]]> Power Module <![CDATA[m1、m2、m3]]>

[0181] By using SPSS software to perform joint analysis and calculation, the utility values ​​of customers for module candidates were obtained, as shown in Table 5.

[0182] Table 5 Utility values ​​of basic module candidates based on conjoint analysis

[0183]

[0184]

[0185] Table 6 Module Candidate-Supplier Information

[0186]

[0187]

[0188] Table 7 Composite Modules - Manufacturer Information

[0189]

[0190]

[0191] Table 8 Product-Distributor Information

[0192]

[0193] Based on the chaotic adaptive nested genetic algorithm, the fitness values ​​of the three composite module division methods are calculated respectively, and the output images of the calculation results of the three division methods are compared. Figure 1 As shown in the figure, the optimal composite module division method is selected as the second scheme with the largest fitness value. After the composite module division method is selected, the corresponding basic module and its module candidate configuration scheme are calculated.

[0194] The upper utility function, lower cost function and fitness value change trend of the second optimal composite module partitioning scheme are output as an image, such as Figure 2 It can be seen that the upper-layer customer utility function shows an upward trend, the lower-layer cost function shows a downward trend, and the fitness value shows an upward trend under the interaction of the two until the optimal value of the PSS module configuration scheme is found, as shown in Table 9.

[0195] Table 9 Optimal configuration scheme

[0196]

[0197] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed in the embodiments of the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for optimizing product service system module configuration with the participation of multiple stakeholders, characterized in that: The following steps are involved: S1. Construct a product service system module configuration optimization model as the upper model; The product configuration scheme in the upper model is as follows: different module candidates are selected from the module candidate set to form different basic modules, and all basic modules form a basic module set; different basic modules are selected from the basic module set to form different composite modules, and all composite modules form a composite module set; different composite modules are selected from the composite module set to form different product solutions, and all product solutions form a product solution set; basic modules are divided into two categories, namely optional modules and mandatory modules, the set composed of all optional modules is the optional module set, and the set composed of all mandatory modules is the mandatory module set; The objective function of the upper model is: in, In the formula, F(x ikjl ) represents the cost-utility function; U i Represents the i-th product configuration scheme PSS i The utility of TC i represents the cost of stakeholders in the i-th product configuration scheme; K represents the total number of composite modules; J represents the total number of basic modules; L j represents the total number of module candidates for the jth basic module; The technical characteristics corresponding to the jth basic module in the product PSS i The weight in μ jl x is the utility of the lth module candidate of the jth basic module that the customer considers; ikjl represents the selection variable of the lth module candidate of the jth basic module of the kth composite module of the ith product configuration scheme. When the product configuration scheme does not contain this selection variable, x ikjl = 0, when the product configuration includes the selected variable, x ikjl =1; α is the profit cost rate of the upper structure to the lower structure benefit analysis; TC i S ,TC i R ,TC i B denote the minimum total cost of the lower-level suppliers, manufacturers, and distributors, respectively; S2. constructing a stakeholder cost optimization model as a lower-level model, wherein the stakeholders include suppliers, manufacturers, and distributors; S3, construct a double-layer model with the upper model constructed by S1 and the lower model constructed by S2, and solve the double-layer model based on the chaotic adaptive nested genetic algorithm.

2. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 1, characterized in that: The supplier's objective function in the lower model is: The constraints are: In the formula, S jl The total number of suppliers who provide candidates; is a decision variable. When the sth supplier provides the lth module candidate of the jth basic module, otherwise represents the material cost of the lth module candidate of the jth basic module in the sth supplier; d jl PSS i The demand for the lth module candidate of the jth basic module; represents the unit ordering cost of the lth module candidate of the jth basic module from the sth supplier; represents the unit inventory cost of the lth module candidate of the jth basic module in the sth supplier; represents the unit transportation cost of the lth module candidate of the jth basic module in the sth supplier; represents the allocation function of the selected suppliers on the supply amount; is the total number of suppliers who can provide the lth module candidate of the jth basic module; Q j PSS i The demand for the jth basic module.

3. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 2, characterized in that: The constraints on suppliers in the lower model also include: The constraint that two modules cannot be selected at the same time is calculated as follows: In the formula, The selection variable representing the lth module candidate of the j1th basic module of the kth composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, Represents the selection variable of the lth module candidate of the j2th basic module of the k'th composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, The constraint condition for selecting two modules at the same time is calculated as follows: In the formula, The selection variable representing the lth module candidate of the j3rd basic module of the kth composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, The selection variable representing the lth module candidate of the j4th basic module of the k'th composite module of the i-th product configuration scheme. When the product configuration scheme does not contain this selection variable, When the product configuration contains this selection variable, May represents a collection of optional modules.

4. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 1, characterized in that: The objective function of the manufacturer in the lower model is: The constraints are: In the formula, R k represents the total number of manufacturers; represents the operating fixed cost of the rth manufacturer of the kth composite module; represents the unit operating variable cost of the rth manufacturer to assemble and manufacture the lth module candidate of the jth basic module; d k represents the demand for the kth composite module; It represents the process variation cost of the rth manufacturer to assemble and manufacture the required module j as the module candidate l; is a decision variable. When the rth manufacturer provides the kth composite module, otherwise represents the proportion of the demand for the kth composite module produced by the rth manufacturer; Select the total number of manufacturers for the kth composite module; is the production capacity range of the rth manufacturer for the kth composite module; Must represents the set of mandatory modules.

5. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 1, characterized in that: The objective function of the distributor in the lower model is: The objective function of the distributor in the lower model is: The constraints are: y b ∈{0,1} Where B is the total number of distributors for each product; N B Select the number of distributors for each product; RC b represents the unit ordering cost of the b-th distributor; HC b represents the unit inventory cost of the product in the bth distributor; TCD b represents the unit transportation cost of the product in the b-th distributor; f(y b ) represents the proportion function of the b-th distributor in distributing the product; d represents the total demand for the product; W b represents the distribution capacity of the bth distributor for the product; y b is a decision variable. When the bth distributor decides to distribute, y b =1, otherwise y b =0.

6. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 1, characterized in that: Step S3 includes: S31, parameter initialization: determine the size of the initial population M, the genetic generation G, the chaotic state control parameter μ, the elite retention number m, and the upper and lower bounds of the adaptive crossover and mutation operators; S32, chaos initialization population: performing chaos optimization operation on the initial populations of the upper and lower layer models respectively; S33, sort all individuals in the population according to their fitness values ​​from large to small, select the first m individuals with the best fitness values ​​to be directly inherited to the next generation; for the remaining Mm individuals, use a mixed ranking strategy to select individuals participating in crossover and mutation from these individuals; The calculation formula for the fitness value of the i-th individual is as follows: S34, combining the selection mechanism of elite retention and mixed ranking to screen out the parent individuals for subsequent crossover and mutation operations; S35, adaptive crossover and mutation operation: calculate the corresponding crossover and mutation probabilities based on the fitness value of the individual, and complete the crossover and mutation operations; The sire is determined by roulette selection, and the probability of being selected as the sire is determined according to the fitness value of the individual; the probability of the i-th individual being selected can be calculated as: Adaptive crossover probability p c , adaptive mutation probability p m The calculation formulas are: In the formula, p c0 、p c1 are the initial adjustment factors for the lower and upper bounds of the crossover probability, which together determine the fluctuation range of the crossover probability; p m0 、p m1 are the initial adjustment factors for the lower and upper bounds of the mutation probability, respectively, which limit the upper and lower limits of the mutation probability; f is defined as the individual fitness value; f max is the optimal fitness value of individuals in the current population; f min is the worst fitness value of individuals in the population; f avg is the average fitness value of individuals in the population; S36. If the number of iterations reaches the upper limit of the genetic algebra G preset in step S31, the algorithm ends and outputs the optimal module configuration solution; otherwise, go to step S33.

7. The method for optimizing product service system module configuration with multi-stakeholder participation according to claim 5, characterized in that: The calculation formula for determining the chaotic variables in the upper model in step S32 is as follows: x k+1 =μx k (1-x k ) In the formula, x k is the result of chaotic variable x after iterating k times; x k+1 is the result of the chaotic variable x after iterating k+1 times; The calculation formula for determining the chaotic variables in the lower model is as follows: t ik =Int(nμx k-1 (1-x k-1 )) t ik is the value of the i-th optimization variable after the k-th iteration in the mathematical model; x k-1 is the result of chaotic variable x after iterating k-1 times; n is the number of module candidates corresponding to each basic module.

Images