Identifying Core Parts in Complex Mechanical Product for Change Management and Sustainable Design

Changes in customer needs are unavoidable during the design process of complex mechanical products, and may bring severely negative impacts on product design, such as extra costs and delays. One of the effective ways to prevent and reduce these negative impacts is to evaluate and manage the core parts of the product. Therefore, in this paper, a modified Dempster-Shafer (D-S) evidential approach is proposed for identifying the core parts. Firstly, an undirected weighted network model is constructed to systematically describe the product structure. Secondly, a modified D-S evidential approach is proposed to systematically and scientifically evaluate the core parts, which takes into account the degree of the nodes, the weights of the nodes, the positions of the nodes, and the global information of the network. Finally, the evaluation of the core parts of a wind turbine is carried out to illustrate the effectiveness of the proposed method in the paper. The results show that the modified D-S evidential approach achieves better performance regarding the evaluation of core parts than the node degree centrality measure, node betweenness centrality measure, and node closeness centrality measure.


Introduction
With the slowdown in global economic growth, different strategies have been proposed in recent years, such as "Industry 4.0", "Made in China 2025", etc.These strategies point out that the key to affecting the level of manufacturing is the competitiveness of the equipment manufacturing industry [1,2].There are many ways to improve the competitiveness of the equipment manufacturing industry.On the one hand, optimizing material flow and applying lean products in the production process can reduce the company's production costs and help obtain competitive prevalence [3,4].On the other hand, improving the design capability of the equipment manufacturing industry is the key issue for improving the company's production level and obtaining competitive advantages [5].This paper focuses on how to improve the design capability of the equipment manufacturing industry.In fact, the vast majority of products in the equipment manufacturing industry are complex mechanical products, such as large ships, weapons, and high-speed trains [6].Under this background, trying to improve the ability to design complex mechanical products is not only an issue that companies must face when meeting market challenges, but also a key approach to improving national core competitiveness [7].
Changes in customer needs often occur in the development of complex mechanical products, and are difficult to avoid [8,9].On the one hand, changes in customer needs are regarded as an important source of product innovation and sustainable design [10][11][12].On the other hand, changes in customer needs may bring severely negative impacts on product design, such as extra costs and delays, even an avalanche [13,14].Generally, a need change is firstly mapped to a change of one part or some parts.In addition, due to the large number of parts and the complicated relationship between parts, a change occurring in one part may cause a series of changes in other parts that are connected to it, or even spread farther, eventually resulting in an avalanche effect.So, effectively dealing with changes in customer needs, accelerating the speed of the change response, and improving the agility of product development have become the main ways to improve the competitiveness of enterprises in this area [15].
The existing research studies on changes in customer needs and engineering change management mainly focus on change propagation [16][17][18], change impact [19,20], and change control [21,22].There is no study of engineering change prevention.However, it is of great significance to carry out engineering change prevention for reducing the change impact, the market response time, and the product design costs.In the research studies on engineering change prevention, it is important to identify the exact core parts that control the propagation direction and the impact scope of the changes.Identifying core parts and focusing attention on them are important contents regarding change management, which is of great significance for reducing the impact of changes and preventing the avalanche effect.In this study, a core part is referred to as a part that plays a decisive role in change propagation speed and scope when changes in customer needs occur.Based on those analyses, which are oriented to the changes in customer needs, it is necessary to design an effective method that can be used to evaluate the core parts of complex mechanical products when changes in customer needs occur.
In recent years, studies have shown that systems with hierarchical structures and complex internal relations can be modeled by complex network theory [23][24][25].For example, complex networks were used to deal with intertwined problems among the logistics network, information network, and financial network.Combining the product characteristics, a structure network model is constructed based on the complex network theory [26][27][28].On the basis of the model, the evaluation of the core parts of a complex mechanical product can be regarded as the evaluation of core nodes in the network model.In general, node centrality is used to evaluate the influence degree of a node on other nodes, in which the bigger the node centrality, the more critical the node.
The node centrality assessment methods can be divided into four categories in the existing research studies.
(1) The node centrality assessment methods that consider the neighbors number of the node.Chen et al. proposed a semi-local centrality measure to evaluate the performance by using the spreading rate and the number of infected nodes [29].Chen etc. proposed a local ranking algorithm named 'cluster rank', which takes into account not only the number of neighbors and the neighbors' influences, but also the clustering coefficient [30].Kitsak et al. proposed a K-shell decomposition method, which provides a route for an optimal design of efficient dissemination strategies [31].Among them, the semi-local centrality and cluster rank algorithm are local assessment methods of node centrality, which cannot make a global assessment of node centrality, and easily lead to inaccurate assessment results.The K-shell decomposition method cannot effectively evaluate the node centrality for a complex network with a large number of nodes.
(2) The node centrality assessment methods depend on the path.Dolev et al. presented algorithms for computing the Routing Betweenness Centrality (RBC) of all the individual vertices in the network, and algorithms for computing the RBC of a given group of vertices [32].Freeman developed three methods for each concept: one absolute and one relative measure of the centrality of a position in a network, and one relenting the degree of centralization of the entire network [33].Brandes proposed a betweenness centrality method, which substantially increased the range of networks for which centrality analysis is feasible [34].Most of these methods are global optimization algorithms.The main disadvantage of these methods is their high time complexity.
(3) The node centrality assessment methods based on eigenvectors.Bonacich et al. proved that there is a necessary addition to the concepts of the degree, closeness, and betweenness centrality, as distinguished by Freeman [35].Brin and Page addressed building a practical large-scale system that could exploit the additional information present in hypertext [36].Lempel and Moran proved that the stochastic approach for link-structure analysis (SALSA) is equivalent to a weighted in-degree analysis of the link-structure of World Wide Web sub graphs, making it computationally more efficient than the mutual reinforcement approach [37].Chen et al. put forward the understanding model of complex product systems (CoPS) from the viewpoint of the width and depth of technology and the technology development model of CoPS [38].These methods require multiple iterative convergences, and most of are better for a directed network.
(4) The node centrality assessment methods based on node removal and shrinkage.Dangalchev proposed a new characteristic (residual closeness), which can measure the network resistance that is presented [39].Li et al. used the sum of the reciprocal of all the geodesics between a pair of disconnected actors caused by deleting an actor (set) to characterize the extent to which the network has been destroyed [40].Tan et al. proposed a node contraction method of evaluation of node importance in complex networks based on a new evaluation criterion [41].These methods are mostly used in power networks, internet, transportation networks, etc.
Based on the above analyses, it is found that the factors affecting the node centrality include three aspects: the number of the node neighbors, the strength of the node, and the position of nodes in the network.However, in addition to the above deficiencies, the main shortage of the existing methods is that most of the methods evaluate node centrality from a single point of view, which resulting in an inability to achieve a systematic and accurate assessment of node centrality.In order to solve these issues, a modified D-S evidential approach is proposed in this study.The effectiveness of the proposed method is verified by comparing it with other three methods of node centrality evaluation.Moreover, the case study indicates that the proposed method can be used to accurately evaluate the core parts.
The rest part of this paper is organized as follows.An undirected weighted network model is constructed in Section 2. Section 3 proposes a modified D-S evidential approach for core nodes evaluation.Section 4 is a case study; in this section, the evaluation of core parts of a wind turbine is used to verify the proposed approach.Finally, conclusions are represented.

The Construction of a Network Model
The expression of product structure is the basis of identifying core parts.Due to complex mechanical products containing a large number of parts and there being complex correlations between different parts, complex network theory was used to model the product structure.The constructed network model can systematically and quantitatively express the product structure [28].In the network model, the parts are defined as nodes, while the relationship between parts are defined as edges, and the strength of the edges is used to express the network weight.Therefore, the network model G s is defined as follows: V = {v i , i = 1, 2, . . ., n} is the set of nodes, where v i is the i-th part of a complex mechanical product; E = e ij , i, j = 1, 2, 3, . . ., n, i = j is the set of edges, where e ij denotes the relationship between v i and v j ; and W = w ij , i, j = 1, 2, 3, . . ., n, i = j is the strength of edge weights, where w ij denotes the strength of e ij .
The determination of weights is a key step for constructing a network model.In this study, the weights of the network model are determined by considering functional and structural aspects.The comprehensive weight can be calculated as follows: where α and β correspond to the weights of w f ij and w s ij , respectively.w f ij is the total functional load value that the v j can bear relative to v i .ws ij denotes the structural correlation strength between v i and v j .The weight w f ij is determined by the functional load, and the value is calculated as follows: where f h is the functional load that v j can bear for the realization of function h, and F h is the total functional capacity.The weight w s ij is determined by the structural constraints, such as size constraints, shape constraints, etc., and it can be calculated as follows: where w s ij−l denotes the correlation strength of the l-th kind of structural correlations.L represents the total number of types of structural correlations.
It should be noted that the greater the correlation strength between two nodes, the greater the weight of the node.

A Modified D-S Evidential Approach for Core Nodes Evaluation
A core part is refers to a part that plays a decisive role in change propagation speed and scope when changes occur.Determining whether a part is a core part is based on its influence on other parts.Correspondingly, the identification of core parts refers to identifying those nodes that have an important influence for the network model.Due to the node centrality indicating the influence degree of the node on other nodes, the identification of the influence nodes can be seen as the calculation of the node centrality.The greater the node centrality, the greater the influence degree of this node on the other nodes in network.
The classical D-S evidential method is a node centrality calculation method in the complex network theory, which is often used to deal with uncertain problems [4].The classical D-S evidential method has considered the number of the node neighbors and the strength of the node for node centrality calculation; however, this method ignores the position of the node and the global structure of the network.Therefore, in order to systematically and scientifically evaluate the node centrality, a modified D-S evidential approach is proposed in this study.
The modified approach is embodied in two aspects.Firstly, the position information of the node is considered by introducing the semi-local centrality measure.According to the definition and calculation formula of the semi-local centrality method [29], the semi-local centrality method reflects the position information of the node in the network by considering both the nearest and the next nearest neighbors.In this study, the value of node centrality is calculated by the D-S evidential approach, and is further processed by the idea of the semi-local centrality method.Secondly, the global information of the network is considered by adding the calculation of the shortest path.In this study, the basic probability assignments (BPAs) are calculated not only based on the node degree, but also based on the shortest path.In short, the calculation results of the modified D-S evidential approach not only reflect the degree and weight of the node, they also reflect the position information and global information of the network, which makes the results more systematic and accurate.
This section has two parts: Section 3.1 includes a definition of all the parameters involved in the proposed method, and Section 3.2 is the detailed calculation process of the modified D-S evidential approach.

Parameters Definition
All of the parameters that are involved in the modified D-S evidential are defined in this section, as shown in Table 1.The table contains the notations of all of the parameters, as well as a specific description of the parameters.The evaluation indices for the degree of a node as "a core node" and "a non-core node" m ki (high)/m ki (low) The probabilities of "high" and "low" influence for the degree of node v i α i The corrected parameter of node degree calculation k i The degree of node v i k max /k min The maximum and minimum values of a node degree in the network The probabilities of "high" and "low" influence for the shortest path between v i and other nodes The shortest path between v i and v j d imax /d imin The maximum and minimum values of the shortest path between v i and other nodes Γ(i) The set of nodes that connected (nearest neighbors) with v i p(k i ) The degree distribution H(j) The set of nodes that is a degree is lower than The BPAs of the node v i with respect to the degree and shortest path M(i) Integration of the BPAs of v i m i (high)/m i (low) The integration probabilities of "high" and "low" influence for v i m i (θ) The probability of "high" or "low" The final probabilities of "high" and "low" of v i σ/δ Constants ENC(i) The evidential node centrality of v i ENC(i) min The minimum value of evidential node centrality NC(i) The node centrality of v i Γ(i) The set of nodes that is next nearest neighbor of v i

The Calculation Process of Node Centrality
A frame of discernment η is defined as follows: The basic probability assignment (BPA) for each node is given as follows.
w ij (10) Suppose node v i with degree k i follows a degree distribution p(k i ).Thus, α i is defined as follows: The BPAs of the v i are calculated based on the degree and shortest path, as follows: Integration of the BPAs of v i is calculated as follows: Generally, in the Equation ( 13), let m i (θ) assign to m i (high) and m i (low) averagely, then: Therefore, the node centrality of v i is shown as follows: In order to ensure that ENC(i) is a positive number, the numerical treatment and normalization are represented as follows: According to the computing idea of a semi-local centrality measure, the value of ENC (i) is further optimized, and the value of NC(i) is obtained.The value of NC(i) is equal to the sum of ENC (i), the nearest and the next nearest neighbors' ENC' of node v i .
NC(i) is calculated as follows: The higher the value of NC(i), the bigger the node centrality of v i .In other words, the parts with bigger values of NC(i) are regarded as core parts.

Case Study
The evaluation of core parts of a wind turbine is used to verify the proposed method.Specifically, the purpose of this study is to verify that the modified proposed has a positive effect on reducing the impact of changes in customer needs on wind turbine design, preventing the change avalanche, decreasing the design delay, and reducing extra costs.
In this section, firstly, the structure of the wind turbine is described and transformed into a network structure based on complex network theory.Secondly, the node centrality is calculated with the proposed method; this means that the core parts are evaluated, too.Finally, the proposed method and calculation results are discussed by comparing them with the existing methods.

The Construction of the Network Model
According to the research in Section 2, the structure of the wind turbine is described and transformed as a network model by using complex network theory.In the network model, the parts of the wind turbine are defined as nodes, the relationships between the parts are defined as the edges, and the strength of the edges constitute the network weight.
By considering the representation of this study, the main parts of the wind turbine (as listed in Table 2) are selected.According to the research in Section 2 and combining with the characteristics of wind turbines, the weights of the structure of a wind turbine are determined by considering both the functional and structural elements.Due to the influence of different factors, the weights of the factors evaluated by experts and calculated based on the entropy method are enumerated in Table 2.
For the calculation of the weights, firstly, 10 designers and experts took these two factors (functional and structural) as the evaluation index of the relationship between parts, and then judged the correlation strength between the parts using linguistic variables.The relationship strength of the linguistic variables are divided into six levels, including very strong, strong, middle, weak, very weak, and no.The value of these six levels are 0.9 (very strong), 0.7 (strong), 0.5 (middle), 0.3 (weak), 0.1 (very weak), and 0 (no), respectively.Secondly, on this basis, the edge weights are calculated according to Equations ( 2)-( 4).The calculation results are shown in Table 3.Based on this, the network model is represented by UCINET 6, which systematically and intuitively depicts the structure of the wind turbine (as shown in Figure 1).

The Evaluation of Core Parts
Based on the constructed network model, the

( )
NC i values of the nodes are calculated by using the modified D-S evidential approach according to equations ( 4)- (18).The calculated results are shown in Figure 2.

The Evaluation of Core Parts
Based on the constructed network model, the NC(i) values of the nodes are calculated by using the modified D-S evidential approach according to Equations ( 4)- (18).The calculated results are shown in Figure 2.

The Evaluation of Core Parts
Based on the constructed network model, the ( ) NC i values of the nodes are calculated by using the modified D-S evidential approach according to equations ( 4)- (18).The calculated results are shown in Figure 2.  As shown in Figure 2, node v 2 (Rotor hub) has the largest value of node centrality in the network model.From this, it can be seen that the node v 2 is a core node in the structure network model of the wind turbine.When node v 2 changes, it will cause other nodes to change in the structure network model, which will cause the greatest impact on the model.In addition, other nodes with greater centrality include node v 18 (Nacelle support), v 12 (Deflector), v 13 (Rotor lock insert), v 20 (Nacelle frame), and v 14 (Blade bearing), where the corresponding values of NC are 151.954,147.603, 136.178, 134.664, and 127.803.These nodes can be seen as core nodes in the network model.Correspondingly, rotor hub, nacelle support, deflector, rotor lock insert, nacelle frame, and blade bearing are core parts of the wind turbine.It can be seen that the larger the node centrality of the node, the higher influence of the node.
In order to reduce the impact of changes, as well as decrease the market response time and design costs, the change buffer of core parts should be added during the wind turbine design process.Specifically, in the design process of the wind turbine, the core parts of the wind turbine should have a certain design margin as a change buffer.The change is amplified when the change of a part that is caused by changes exceeds its change buffer; otherwise, the change is absorbed.

Discussion
In order to verify the effectiveness and the advantages of the proposed method, the analysis and verification of the calculation results are carried out from two aspects.On the one hand, the validity of the proposed method is verified by analyzing and discussing the results obtained in this paper.On the other hand, the results obtained in this paper are compared with the results obtained by other methods to verify the advantages of the proposed method.Due to differences in the calculations, the results cannot be directly compared and discussed.Therefore, in order to verify the advantages and effectiveness of the proposed method, in this paper, the SI model of a complex network was introduced.
The SI model is used to describe the change propagation process in the network.It can also be understood as the description of the nodes' state (unchanged node and changed node) in the network at each moment of the change propagation process.The SI model is described as follows.
Each node in the SI model is allowed to have two states, namely the unchanged node and the changed node.The change propagation process refers to the process in which a changed node causes a change in an unchanged node.The change propagation capability is the probability that the changed node causes a change in the other unchanged nodes; that is, the change propagation rate can be used to reflect the change propagation capability of a node.The greater the change propagation rate of a node, the more changes it will cause.
The change propagation rate between two nodes is defined: where λ ij denotes the propagation rate from the changed node v j to the unchanged node v i .w M denotes the maximum weight in the structure network.
The change propagation rate of node v i can be calculated as follows: where V t denotes the set of changed nodes in the neighbor of the node v i .
In addition, the change propagation threshold λ c of structure network is defined as follows: where Λ m represents the feature vector value of the parts' relationship matrix.In this paper, the SI model was used to verify the proposed method through the change propagation process in the network.
Firstly, the calculation results of the modified D-S evidential are analyzed and discussed.
The top 20% nodes of node centrality value in the network are selected for analyzing; that is, v 2 , v 18 , v 12 , v 13 , v 20 , and v 14 are selected for analyzing.These nodes are selected as the initially changed nodes and input into the SI model.In addition, the change propagation rates at each time under different initially changed nodes are calculated, and then, the number of changed nodes at each time can be determined.Based on this, the change propagation process is simulated.The propagation and impact of the initial change nodes on the network are analyzed based on the simulated results.The comparison of the change propagation influence of different initially changed nodes (v 2 , v 18 , v 12 , v 13 , v 20 , and v 14 ) is shown in Figure 3.As shown in Figure 3, the abscissa represent the times of change in the propagation process.The ordinate represents the number of changed nodes, which is used to describe the change propagation capability of the initially changed node in the structure network of a wind turbine.The change propagation process of these six nodes as initially changed nodes are shown in Figure 3.The six curves represents the number of changed nodes at different times of different initially changed nodes.t 0 is the initial time at which a change occurs.The change propagation is completed at t 8 .
As shown in Figure 3, node v 2 , as the initially changed node in the structure network of a wind turbine, has the fastest propagation speed compared to others nodes as initially changed nodes.Moreover, node v 2 causes the most changes to the nodes in the network (23 nodes are changed in the network).Therefore, node v 2 is considered to have the strongest change propagation capability.The change propagation speed of node v 18 as the initially changed node is slower than the change propagation speed of node v 2 ; meanwhile, node v 18 also causes the most changes to the nodes in the network (the same as node v 2 ).The change propagation speed of node v 12 as the initially changed node is slower than that of node v 2 and node v 18 , and is faster than node v 13 , node v 20 , and node v 14 .In addition, the number of node changes causing by node v 12 is less than that of node v 2 and node v 18 , and is more than node v 13 , node v 20 , and node v 14 .Finally, the change propagation speed and changed nodes are the same as those of node v 13 , node v 20 , and node v 14 as initially changed nodes, eventually causing 16 nodes to be changed in the network.
The compared results are in agreement with the results of the proposed method.From the analyses results, it can be seen that the nodes identified by the modified D-S evidential approach have a greater impact on the network.In other words, the proposed method is effective.
Secondly, in order to verify the advantages of the proposed method, the proposed method was compared with other node centrality evaluation methods.
In this section, the degree centrality (DC), the betweenness centrality (BC), and the closeness centrality (CC) [42] are compared with the proposed method in this study.Firstly, we calculate the value of the DC, BC, and CC of nodes in the network model of the wind turbine.Secondly, the top 20% nodes regarding node centrality value, as calculated by these methods, are selected for analysis.Some nodes are selected as the initially changed nodes and input into the SI model.In addition, the change propagation rates at each time under different initially changed nodes are calculated, and then, the number of changed nodes at each time can be determined.Based on this, the proposed method that is verified has advantages over the other methods.
The DC, BC, and CC of wind turbine parts are calculated, respectively.The calculation results are shown in Table 4.As shown in Table 4, the degree centrality, betweenness centrality, and closeness centrality of nodes in the network model are calculated.It can be seen that the node with the largest degree centrality is node v 13 , which is followed by nodes v 12 , v 2 , v 18 , v 20 , and v 8 .Meanwhile, the node with the largest betweenness centrality is node v 17 , followed by nodes v 12 , v 18 , v 3 , v 20 , and v 2 .In additional, the node with largest closeness centrality is node v 18 , followed by node v 17 , v 11 , v 19 , v 7 , and v 4 .Based on the calculation results, the ranking of core nodes in a network with the modified D-S evidential approach, DC, BC, and CC are obtained.The top 20% of nodes are selected for sorting, which are shown in Table 5.
Table 5.The ranking of core nodes in the network model by different methods.

Rank
The Proposed Method DC BC CC Obviously, the calculation results regarding node centrality with these four methods are different; however, the advantages of the four methods cannot be directly compared and discussed.For this reason, the SI model is introduced.
According to the ranking results of the modified D-S evidential approach and DC (as shown in Table 5), it can be seen that there are seven nodes (v 2 , v 18 , v 12 , v 13 , v 20 , v 14 , v 8 ) in the two rankings.Due to the effectiveness of the modified D-S evidential approach being proven (as shown in Figure 3), the rank of the change propagation capability of these nodes Therefore, the node v 14 (the weakest change propagation capability in ranking of the modified D-S evidential approach) is selected to compare with node v 8 (node v 8 has not been analyzed).According to the ranking results of the modified D-S evidential approach and BC, it can be seen that there are eight nodes (v 2 , v 18 , v 12 , v 13 , v 20 , v 14 , v 8 , v 17 ) in the two rankings.Based on the existing analysis, similarly, the node v 14 is selected to compare with node v 17 (node v 17 has not been analyzed).According to the ranking results of the modified D-S evidential approach and BC, it can be seen that there are 11 nodes (v 2 , v 18 , v 12 , v 13 , v 20 , v 14 , v 17 , v 11 , v 19 , v 7 , v 4 ) in the two rankings.Based on the existing analysis, similarly, node v 14 is selected to compare with node v 11 , node v 19 , node v 7 , and v 4 (these nodes have not been analyzed).
Therefore, the propagation influence of six groups of nodes are compared (as shown in Figure 4).These nodes are selected as the initially changed nodes, and are input into the SI model.After that, the change propagation process is simulated.The propagation and impact of the initial change nodes on the network are analyzed based on the simulated results.
The simulation results are shown in Figure 4.As shown in Figure 4, the abscissa represents the times of change propagation.The ordinate represents the number of changed nodes, which is used to describe the change impact of the initially changed node in the wind turbine structure network.
In Figure 4b,c, it can be seen that the node v 14 propagates more quickly than node v 7 and v 11 .The propagation scopes corresponding to node v 14 , node v 7 and node v 11 are the same.In Figure 4a,d-f, it can be seen that node v 14 propagates more quickly and infects more largely than v 4 , v 19 v 8 , and v 17 .However, considering the change propagation speed and change impact scope, it is easy to see and know whether the comprehensive change propagation capability of node v 14 is greater than that of nodes v 4 , v 7 , v 8 , v 11 , v 17 , and v 19 .
Moreover, the two nodes with the largest node centrality calculated by the four methods are selected as the initially changed nodes and are input into the SI model.The propagation and impact of the initially changed nodes are analyzed based on the simulated results, as shown in Figure 5.In Figure 4b,c  In Figure 5, the four curves represents the number of changed nodes of the four methods at different times.t 0 is the initial time when the change occurs.The change propagation is completed at t 7 .It should be noted that when the change propagation is completed, the propagation scopes corresponding to the proposed method and betweenness centrality are same.However, considering the change propagation speed and change impact scope, it can be noticed that the proposed method is better than the other three methods.
Based on the above analysis, it can be seen that the nodes identified by the modified D-S evidential approach have a greater impact on the network.In other words, the proposed method has advantages over other methods.
In summary, the change of these core nodes is the root cause of avalanche transmission.Controlling and managing core nodes are the core of change management.It is an important measure to reduce the impact of customer's needs changes, and thus decrease the market response time and design costs.Therefore, identifying core parts and giving large buffers of core parts are both effective measures to reduce the change influence.
Moreover, the two nodes with the largest node centrality calculated by the four methods are selected as the initially changed nodes and are input into the SI model.The propagation and impact of the initially changed nodes are analyzed based on the simulated results, as shown in Figure 5. t .It should be noted that when the change propagation is completed, the propagation scopes corresponding to the proposed method and betweenness centrality are same.However, considering the change propagation speed and change impact scope, it can be noticed that the proposed method is better than the other three methods.
Based on the above analysis, it can be seen that the nodes identified by the modified D-S evidential approach have a greater impact on the network.In other words, the proposed method has advantages over other methods.
In summary, the change of these core nodes is the root cause of avalanche transmission.Controlling and managing core nodes are the core of change management.It is an important measure to reduce the impact of customer's needs changes, and thus decrease the market response time and design costs.Therefore, identifying core parts and giving large buffers of core parts are both effective measures to reduce the change influence.

Concluding Remarks
Meeting customer needs is an integral part of the sustainable design of complex mechanical products.In order to reduce the impact of changing customer needs and prevent an avalanche, it is of great importance to identify the core parts that have great impact in product design.Complex network theory is used for network modeling and the evaluation of core parts in the study.In evaluation, the core parts are transformed to calculate the centrality of the nodes.An undirected weighted network model is established that is used to systematically describe the structure of the products.Based on the constructed network, the modified D-S evidential approach is proposed, which is used to calculate the node centrality.By introducing the semi-local centrality measure and adding the calculation of the shortest path, the calculation results of the modified D-S evidential approach not only reflect the degree and weight of the node, they also reflect the position

Concluding Remarks
Meeting customer needs is an integral part of the sustainable design of complex mechanical products.In order to reduce the impact of changing customer needs and prevent an avalanche, it is of great importance to identify the core parts that have great impact in product design.Complex network theory is used for network modeling and the evaluation of core parts in the study.In evaluation, the core parts are transformed to calculate the centrality of the nodes.An undirected weighted network model is established that is used to systematically describe the structure of the products.Based on the constructed network, the modified D-S evidential approach is proposed, which is used to calculate the node centrality.By introducing the semi-local centrality measure and adding the calculation of the shortest path, the calculation results of the modified D-S evidential approach not only reflect the degree and weight of the node, they also reflect the position information and global information of the network, which makes the results more systematic and accurate.Moreover, it has been proven that the modified D-S evidential is effective.Furthermore, the study has proved that the nodes identified by the modified D-S evidential approach have a greater change propagation capability and a greater change impact on the network.In other words, the proposed method has advantages over other methods.
In addition, the redundant design of core parts is the research direction for the further exploration oriented toward changes in customer needs and engineering.

Figure 1 .
Figure 1.The structure network model of the wind turbine.

Figure
FigureThe structure network model of the wind turbine.

Figure 1 .
Figure 1.The structure network model of the wind turbine.

Figure 2 .
Figure 2. The NC(i) of nodes in network model.

Figure 3 .
Figure 3.The comparison of propagation influence of different initially changed nodes.

Figure 4 .
Figure 4. Comparison of the propagation influence of different initially changed nodes.

Figure 4 .
Figure 4. Comparison of the propagation influence of different initially changed nodes.

Figure 5 .
Figure 5.The influence of customer's needs changes with the four methods.

Figure 5 .
Figure 5.The influence of customer's needs changes with the four methods.

Table 1 .
The description of parameters.BPAs: basic probability assignments.

Table 2 .
The weights of the influencing factors.

Table 3 .
The weights of the network model.

Table 4 .
The degree centrality (DC), betweenness centrality (BC), and closeness centrality (CC) of nodes in a structure network model.
, it can be seen that the node 14 v propagates more quickly than node 7 v and 11 v 17 , and17v .However, considering the change propagation speed and change impact scope, it is easy to see and know whether the comprehensive change propagation capability of node 14