Yin Yong Wang Shuxin

(1School of Information Engineering, Wuhan University of Technology, Wuhan 430070, China)(2Shenzhen Research Institute, Wuhan University of Technology, Shenzhen 518000, China)(3School of Intelligent Manufacturing and Electronic Engineering, Wenzhou University of Technology, Wenzhou 325035, China)

Abstract：To compensate for the limitations of previous studies, a complex network-based method is developed for determining importance measures, which combines the functional roles of the components of a mechatronic system and their topological positions.First, the dependencies among the components are well-represented and well-calculated.Second, a mechatronic system is modeled as a weighted and directional functional dependency network(FDN), in which the node weights are determined by the functional roles of components in the system and their topological positions in the complex network whereas the edge weights are represented by dependency strengths.Third, given that the PageRank algorithm cannot calculate the dependency strengths among components, an improved PageRank importance measure(IPIM)algorithm is proposed, which combines the node weights and edge weights of complex networks.IPIM also considers the importance of neighboring components.Finally, a case study is conducted to investigate the accuracy of the proposed method.Results show that the method can effectively determine the importance measures of components.

Key words：importance measure; mechatronic system; dependency; complex network theory

Studies are ongoing for determining the importance measures of components in mechatronic systems.A well-known study is the theory of binary coherent systems[1-3].Zio et al.[4]improved the Birnbaum importance measure in binary and multistate systems by using the Markov or semi-Markov process.Subsequently, several other studies on multistate systems have been reported.Recently, methods for determining component importance have been extended to various variants considering different scenarios, such as cost-based importance measures considering the cost of maintaining each component[5], importance measures for systems with degrading components[6], importance measures for components with reconfigured systems[7], and many others[8-10].

Although studies on component importance measures have received considerable interest, some challenges still remain.For example, a serious drawback of the Birnbaum importance measure is that it ranks only individual components.The differential importance measure overcomes this challenge, but it does not consider component aging.For the state-space model, the expenditure to analyze a mechatronic system with a large state space is extremely high.Moreover, existing studies assume that the system components are independent, and none of these studies considers the influence of functional dependency.

Recently, studies on complex networks have gradually received research attention.Li et al.[11]systematically reviewed the literature to investigate the applications of complex networks in the manufacturing field.They found that the following two issues are their main negative factors: 1)Most previous studies on importance measures did not consider the coupling strength among components[12].Moreover, past research only considered models with unweighted networks.Alternatively, each edge was equally weighted.2)According to previous studies, the roles of neighboring nodes may affect a node’s degree of importance[13-14].That is, the importance of a node can be greatly affected by its neighbors, but it was not considered in previous studies.

To address the limitations in past research, this study investigates the dependency-based importance measures of components in mechatronic systems with complex network theory.Its main contributions are to:

1)Address the issue of characterizing functional dependencies among components.Although considerable efforts have been made to model functional dependencies among components, their adequate representations and calculations are still needed.

2)Model a mechatronic system as a functional dependency network(FDN), which is weighted and directional; the components of a mechatronic system are represented with nodes by considering their weights, and the dependencies among components are represented by edges.Furthermore, edge weights are represented by dependency strengths, and the dependencies are directional.

3)Propose a method to identify the important nodes in FDN, considering the functional dependencies among components with an improved PageRank importance measure(IPIM)algorithm.IPIM considers the contributions of the neighbor’s importance.

The paper illustrates the component functional dependency analysis referring to the example from Ref.[15]at first.Next, the mechatronic system is modeled as an FDN and the functional dependencies among components are represented as edge weights.The node weights by combining the functional roles of components in the system and the topological positions of the corresponding nodes are analyzed.Afterwards, an approach called IPIM to find the most important node in FDN is proposed.Finally, the validity of the methodology is given.

1 Component Functional Dependency Analysis

To clarify the dependencies of components in a mechatronic system, we present an example adapted from a previous study as follows[16]:

Suppose supported componentCjis a widget production machine, and supporting componentCisupplies coolant and lubricant fluids toCj.PiandPjare the actual operability levels of componentsCiandCj, respectively, on the condition of dependency linkage fromCitoCj.Suppose thatCjis completely operable when it produces 120 widgets per hour, thenPj(Pi=120)=100 utils.In the case of no supply of fluids fromCi,Cjcan only produce 80 widgets per hour.Suppose the production rate of 80 widgets per hour is worth 55 utils, thenPj(Pi=80)=55 utils when no support is obtained fromCi.This situation implies that the baseline operability level ofCjis 55 utils.Suppose the fluids fromCiare ideal for lowering the operating temperature of the engine and increasing the output ofCj.Without these fluids, the engine’s temperature would rise, resulting in part wearing.Cjwould also decline from its baseline operability level of 55 utils and eventually become completely inoperable(0 utils).

In the example cited above, the authors introduce a conceptual word “util” to express the value of functioning performance for a component, which is a dimensionless number similar to the von Neumann-Morgenstern utility.Dependency strength parameterαijis defined to represent the reliance degree of componentCjonCi, andpjis the baseline operability level of componentCj.Thus,

(1)

wherePiandPjare the actual operability levels of componentsCiandCj, respectively.Eq.(1)reveals that if the dependency strengthαijfromCitoCjis 0, thenCjcan only function at its baseline levelPj.Moreover, ifαijequals 1, thenCjcan function with 100 utils.

However, the support from componentCitoCjmay be intermittent; that is, the interactive relationship between componentsCiandCjis not always continuously effective.Under this condition, two additional indicators, contact frequency and duration, are used to depict dependency strengths accurately.According to a previous study[15], another time-related parameterβijis defined as follows:

(2)

wheredijrepresents the empirical contact duration of the type of functional dependency between componentsCiandCj, andfijdenotes their contact frequency.s(Ci)is the number of times the operation state ofCichanges, whereass(Ci|Cj)is the number of times the operation state ofCjarising fromCichanges.Thus, the normalized dependency strength between componentsCiandCjcan be expressed as

(3)

2 FDN Modeling and Node Weight Analysis

2.1 FDN modeling and edge weight

As discussed in the introduction section, if the components of a mechatronic system are treated as nodes and the dependencies among components as edges, the mechatronic system can be represented as an FDN.To facilitate FDN modeling, we first transfer the bidirectional dependence among components into a one-way dependence, as illustrated in Fig.1.

(a)

In Fig.1, a bidirectional link exists between componentsviandvjwith the dependencies ofeijandeji, respectively.Virtual nodevkis inserted into the network to help transfer the bidirectional link into a one-way one.In Fig.1(a),eijis the dependency from nodesvitovjindicated by a right arrow, whereasejiis that from nodesvjtovipointed out by a left arrow.To avoid the bidirectional link, virtual nodevkis inserted into the network to help transfer the bidirectional link into a one-way one, as displayed in Fig.1(b)whereeijis the dependency from nodesvitovjindicated by a right arrow.ejkis set to be equal toeji, whereasekiis set as 1.Thus, the dependency from nodesvjtovican be calculated asejk*1=eji.We define an FDN as a networked organization of components connected by functional dependencies.Functional dependencies among components are represented by directional edges.Considering the importance of components and the strengths of functional dependencies, an analytical model is established as follows:

FDN=(V,E,W)

(4)

and

(5)

whereVis the set of nodes representing components;nis the total number of nodes;Eis the set of edges;eijrepresents the functional dependency strength(edge weight)between nodesviandvj;Wis the set of node weights forV, andW(vi)represents the weight for nodevi.The edge weight between nodesviandvjcan be calculated using Eq.(3)withαijandβijintroduced in Section 1.

2.2 Node weight analysis

The node weight of a component in a mechatronic system is determined by the functional role of the component in the system and the topological position of the corresponding node in the FDN, as shown in Fig.2.

Fig.2 Factors influencing the node weight

The indexes to evaluate the functional role of a component include the expense, lifetime, mean time between failure(MTBF), and failure probability(FB).These indexes can be obtained through a statistical analysis of failure maintenance data or basic performance parameters in the datasheet for the component.

The topological positions of nodes in an FDN can be determined using the evaluation metrics of complexity theory.For example, when evaluating the importance of nodeviin an FDN, the number of the tail ends(not head ends)linked with nodeviis considered, which is called the outdegree ofvi, as shown in Fig.3.

Fig.3 Outdegree of vi

The outdegree of nodevican be defined as

(6)

whereVout(vi)represents the set of all nodes directed by nodevi.Thekoutiofviin Fig.3 is 3.

Other metrics for evaluating the topological importance of nodes from different aspects include clustering coefficient, closeness centrality, and eigenvector centrality.Specifically, these metrics can be used to evaluate the topological importance of nodes on the basis of the actual application requirements.

Herein, we combine the functional roles of components in a system and the topological positions of the corresponding nodes in an FDN by using an aggregation operator(AO)[17].An AO is a generalized operator that aggregates multiple generalized importance measure indexes into a single comprehensive one to identify the key components in a system.Assuming that setA={a1,a2,…,ai,…,an} exists and the basic measure set of each elementaiinAis {αi1,αi2,…,αim}, the comprehensive value ofaican be expressed as follows:

Ti=AO(αi1,αi2,…,αim)

(7)

Herein, we use the Choquet integral(CI)as the AO because it considers the importance of the critical factors of components and the interrelated and restrictive influences of the factors.On this basis and by combining each importance measure metric shown in Fig.2, CI is used to evaluate the weight of nodevias follows:

i=1,2,…,N

(8)

whereTi={αi1,αi2,…,αim} is the importance index set for nodevi;Nis the number of nodes in the FDN;jrepresents thej-th importance index(j=1, 2, …,m, wheremis the number of importance indexes for nodevi);αi(j)-αi(j-1)is the vector transformation forαi(j), such thatαi(1)<αi(2)< … <αi(m); andμ(Ai(j))is the weight of each importance index for nodevi.

3 IPIM Method

Herein, we propose an approach called IPIM to determine the most important node in an FDN.IPIM is used for evaluating the importance of a component by considering functional dependencyeijdescribed usingαijandβijand node weightW(vi), which are introduced in Section 1 and Section 2.

The PageRank algorithm is used by Google Search and has recorded huge success in sorting websites on the Internet.It is expressed as follows[18-19]:

(9)

where PR(pi)is the PageRank value of the evaluated webpagepi; PR(pj)is the PageRank value of the neighboring webpagepjlinked with webpagepi;p1,p2, …,pNare the considered pages;Nis the number of pages;M(pi)represents the set of pages;Cout(pj)shows the number of outbound links on pagepj; anddis the damping factor.Notably, when calculating the Page-Rank value, pages with no outbound links are assumed to be linked to all other pages in the collection.Therefore, the PageRank values are divided evenly among all other pages.Alternatively, to be fair to pages that have no links, random transitions are added to all nodes on the website, usually with adthat is equal to 0.85, estimated from the frequency with which users continuously follow a link[20].

The PageRank algorithm can be applied to any collection of entities with reciprocal quotations or links.However, when used to evaluate important nodes in an FDN, the algorithm cannot analyze or measure the coupling strength among components.Herein, we introduce IPIM for importance measures to assess node importance in an FDN, which is a directed-weighted complex network.IPIM considers the node weights and edge weights and has three improvements compared with that in Eq.(9).

(10)

wherew(vi)is the normalized value ofW(vi), which is the node weights calculated using Eq.(8).The expression ofw(vi)is given as follows:

(11)

(12)

whereejiis the strength of dependence of the evaluated webpagepion webpagepj, and ∑ejkis the sum of all the dependency strengths between webpagepjand its out-linking webpagepk.

3)In Eq.(9),dis the damping factor indicating the probability that the user reaches the evaluated webpagepiand continues to browse the neighboring webpagepjlinked withpi.It is also a constant value.However, the influence of a component on others in the FDN owing to functional dependencies may vary and spread in a large scope, blocking or even paralyzing the whole network operation.We define the influence propagation probability of nodevias follows[15]:

(13)

where FR(vi)is the functional influence of nodevi, and it can be obtained by adding the normalized values of expense, lifetime, MTBF, and FB.Vi+is the set of neighboring nodes connected to nodevi.Thus, the influence propagation damping of nodevican be calculated as

d(vi)=1-Fp(vi)

(14)

Moreover, IPIM synthetically considers the functional dependency, the node weight of each component, the component functional influence, and the topological structure of the FDN.To calculate the importance of nodevi, the IPIM algorithm is expressed as follows:

(15)

where IPIM(vi)is the IPIM ranking value ofvi;d(vi)is the influence propagation damping;w(vi)is the normalized value ofW(vi)calculated using Eq.(11).

4 Case Study

4.1 System description and network modeling

A case study on an elevator system is conducted to investigate the proposed models and methods.The core units of an elevator include traction subsystem, guide subsystem, car and door subsystem, electrical control subsystem, and safety protection subsystem, as shown in Fig.4.All subsystems have various functional components.Herein, we consider only the traction(a typical example of mechanical subsystems)and safety protection subsystems(a typical example of electrical subsystems)to demonstrate component dependency modeling and its network models.

Fig.4 Core units of an elevator

When all subsystems are devoted to constituting an elevator system and the dependencies among the components in different subsystems are considered, the FDN of the elevator system is obtained, as depicted in Fig.5.

Fig.5 FDN for an elevator system

4.2 Component importance evaluation

In the component dependency network of the elevator system, the strength of the functional dependency of each component is represented aseij.eijis determined using parametersαandβ, whereαrepresents the reliance degree among components in their operation states, andβrepresents the interactive frequency and contact duration of functional dependencies among components.Herein,αcan be collected from the design or manufacturing datasheets or the repair and maintenance databases of the elevator system, andβcan be obtained from experience.As reported in a previous study[15], linguistic terms are used to describe the fuzzy probabilities regarding the interactive frequency and contact duration of functional dependencies by triangular fuzzy numbers.Reliance degree parameterαand the fuzzy probabilities ofβare listed in Tab.1.Here,αandβare derived from the original data collected by Elevator World, Inc., Educational Division.The data must be selected and cleared before use.

Tab.1 Parameters for determining functional dependency strength

The factors used for assessing the functional roles of components include their expense, lifetime, MTBF, and FB, as shown in Fig.2.Such information is collected from the historical databases of the elevator system.The data is studied by eliminating irrelevant items, and the result is presented in Tab.2.The topological positions of components in the elevator system can be evaluated using metrics such as the outdegree and betweenness centrality of the FDN.With the FDN of the elevator system, the metric values for each component can be calculated using Pajek toolkits(a set of tools for complex network parameter calculations)[21], and the outcome is presented in Tab.3.Herein, an iteration termination condition is set as follows:

Tab.2 Historical data assessing the functional roles of components

Tab.3 FDN metric values for components

|IPIM(vn)-IPIM(vn-1)|≤ε

(16)

whereεis a small value to terminate the iterative arithmetic.Whenε= 0.01, the calculation is performed in 15 seconds on a MATLAB platform in a Windows 10 operating system.

When IPIM is used to calculate the importance of a node in an FDN(considering the functional dependency of each node, the node weight of each component, component functional influence, and the topological structure of the FDN), the number of iterations must be determined because the computing complexity increases drastically as the number of components increases.

We compare the importance ranking results obtained using IPIM with those using the Birnbaum binary importance measure reported in Ref.[3], the Birnbaum measure for multistate components in Ref.[4], and the ranking results obtained from the maintenance records of Shanghai Aoan Elevator Co., Ltd.The comparison is shown in Fig.6.

Fig.6 Comparison of ranking results

Thex-axis represents the component number, standing for each component(15 components), whereas they-axis represents the ranking order(1 to 15).The smaller the ranking order, the more important the component is.

The curve presenting the importance ranking obtained from IPIM coincides with that of the maintenance records of Shanghai Aoan Elevator Co.Ltd., indicating a high evaluation accuracy of the IPIM method.Moreover, all the methods identify the traction motor as the most important component, which is consistent with the result obtained from the maintenance records.

However, considerable differences are observed among the ranking results for the reducer, traction rope, and speed governor(Nos.6, 8, and 11, respectively)obtained using different methods.Taking the traction rope(No.8)as an example, IPIM ranks it as important with a ranking of 10, the Birnbaum binary importance measure ranks it as important with a ranking of 2, and the Birnbaum measure for multistate components ranks it as important with a ranking of 6.This discrepancy can be attributed to the Birnbaum binary importance measure method considering only two states for each component, including perfect and completely invalid states.As the traction rope drives the traction wheel to lift the elevator car, the binary importance measure method ranks it high.The Birnbaum measure for multistate components evaluates the traction rope through a stochastic process considering the set of all possible states, the so-called state space.Thus, it provides the traction rope with a fair ranking.However, when we investigate the functional strength of the traction rope with other components, itseijis only 0.66.Furthermore, the expense, lifetime, MTBF, and FB are 150 yuan, 7×104h, 5×104h, and 15%, whereas the outdegree and betweenness centrality are 2 and 0.11, respectively.This results in a smallW(vi)value.All these factors make the IPIM method rank the traction rope as a relatively unimportant component.The close agreement with the ranking results of the maintenance records of Shanghai Aoan Elevator Co., Ltd.validates the correctness of IPIM.

5 Conclusions

1)A complex network-based method is proposed for determining importance by combining the functional role of each component of a system and their topological position in the mechatronic system based on an FDN.

2)In the FDN, the components of the mechatronic system are treated as nodes by considering their weights and the dependencies among components; the dependencies are represented by edges.An IPIM algorithm is established to identify important components by considering the importance of the neighboring components.

3)A case study is conducted to investigate the accuracy of the proposed method, and the results indicate that the method can effectively determine the importance measures of components.

4)Although the results of the proposed method are satisfactory, further studies must obtain:(i)a more accurate model, such as networks of networks, for complex mechatronics and(ii)more parameters for models considering the node weights of components in a system.