Next Article in Journal
Selected Applications of the Theory of Connections: A Technique for Analytical Constraint Embedding
Next Article in Special Issue
A Novel FMEA Model Based on Rough BWM and Rough TOPSIS-AL for Risk Assessment
Previous Article in Journal
Bounded Perturbation Resilience and Superiorization of Proximal Scaled Gradient Algorithm with Multi-Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment for Failure Mode and Effects Analysis Using the Bonferroni Mean and TODIM Method

1
School of Transportation and Logistics, Southwest Jiaotong University, Chengdu 611756, China
2
Department of System Engineering and Engineering Management, City University of Hong Kong, Hong Kong, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(6), 536; https://doi.org/10.3390/math7060536
Submission received: 20 May 2019 / Revised: 9 June 2019 / Accepted: 10 June 2019 / Published: 12 June 2019
(This article belongs to the Special Issue Recent Advances in Modeling for Reliability Analysis)

Abstract

:
As a safety and reliability analysis technique, failure mode and effects analysis (FMEA) has been used extensively in several industries for the identification and elimination of known and potential failures. However, some shortcomings associated with the FMEA method have limited its applicability. This study aims at presenting a comprehensive FMEA model that could efficiently handle the preference interdependence and psychological behavior of experts in the process of failure modes ranking. In this model, a linguistic variable expressed by the interval-valued Pythagorean fuzzy number (IVPFN) is utilized by experts to provide preference information with regard to failure modes’ evaluation and risk factors’ weight. Then, to depict the interdependent relationships between experts’ preferences, the Bonferroni mean operator is extended to IVPFN to aggregate the experts’ preference. Subsequently, an extended TODIM approach in which the dominance degree of failure modes is calculated by grey relational analysis is utilized to determine the risk priority of failure modes. Finally, a practical example concerning the risk assessment of a nuclear reheat valve system is provided to demonstrate the effectiveness and feasibility of the presented method. In addition, a sensitivity analysis and comparison analysis are conducted, and the results show that the preference interdependence and psychological behavior of experts have an important effect on the risk priority of failure modes.

1. Introduction

Failure mode and effects analysis (FMEA), which is a powerful engineering technique for risk management, has proven to be an effective methodology for accident prevention and risk analysis of a system/process to identify potential failure modes and assess the effects of different failure modes [1]. The formal FMEA method was originally proposed by the aerospace industry in the 1960s to satisfy the apparent requirements of reliability and safety [2]. The aim of implementing FMEA is to help identify failure modes that affect a system’s functioning and to correct these failures by giving some suggestions. Currently, FMEA has been widely applied in various fields, such as semiconductor [3], manufacturing [4], and nuclear power [5].
In conventional FMEA, the risk priority number (RPN), which is calculated by the product of severity (S), occurrence (O), and detection (D), is used to obtain the risk priority of failure modes. Although the RPN method is simple, in practical applications, it also suffers some drawbacks, some of which are listed as follows [6,7,8]: (1) Different evaluations of risk factors may produce the same RPN values, but their risk implications are entirely dissimilar; (2) direct and indirect relationships among various failure modes and effects are not considered; (3) obtaining precise and complete evaluations for qualitative risk factors are difficult or even impossible; and (4) uncertainty, subjectivity, and vagueness of the judgments of FMEA team members cannot be managed appropriately. To overcome the above weakness of the RPN method, improved approaches of FMEA have been numerously presented in the past decade and applied to deal with practical problems [6,8,9,10,11,12]. More information with respect to improved FMEA methods is provided in recent surveys [13].
In these improved FMEA methods, FMEA can be considered as multi-criteria decision making (MCDM) problems [14] as it takes into account many risk factors. In order to avoid an individual bias, FMEA activity is usually performed based on a team and should be regarded as a group decision behavior [10]. Therefore, FMEA activity can be considered as multi-criteria group decision making (MCGDM) problems [15]. Generally, the implementation of group-based FMEA consists of three stages, which includes failure modes’ evaluation, determination of the weights of risk factors, and failure modes ranking [16]. This study focuses on the following two stages: Failure modes evaluation and failure modes ranking.
Failure modes evaluation involves expression and aggregation expert preferences. To express precisely and completely experts’ preferences, many alternative methods, which incorporate fuzzy set theory and its extended forms [7,15,17,18,19], have been introduced into FMEA to deal with the uncertainty and ambiguity of experts’ judgment. However, no or little attention has been paid to the way of aggregating experts’ preferences, which usually impact on the final ranking of failure modes. Generally, before the ranking of alternatives, experts’ preferences are aggregated to collective preferences to prevent the loss of information [20], which indicates that the way of aggregation plays an important role in the aggregation preferences of MCGDM problems [21]. Therefore, it is a meaningful work to explore whether the method of aggregating preferences affects the final ranking results.
On the other hand, failure modes’ ranking aims at using an MCDM technique to obtain the risk priority of failure modes. In reality, FMEA is essentially a multi-criteria group decision making problem [15], and many uncertainties exist in an FMEA team expert’s subjective judgments and qualitative assessments due to different expertise and backgrounds. Meanwhile, many existing studies involving behavior experiments have indicated that the decision maker is not completely rational in a practical decision making process [22,23,24], and the psychological behavior of the decision maker, which plays an important role in decision analysis, should be taken into account in the decision making process [25]. However, most of the existing improved FMEA approaches based on the MCDM technique neglect the decision maker’s psychological behavior in the process of failure modes ranking. Consequently, how to determine the risk priority of failure modes based on the decision maker’s psychological behavior is a valuable research topic.
Based on the above analysis, the ranking results of failure modes will inevitably be influenced by the way preferences are aggregated and the decision maker’s psychological behavior. This study presents an integrated method by combing the Bonferroni mean and TODIM (an acronym in Portuguese interactive and multi-criteria decision making) to deal with the ranking of failure modes. In addition, interval-valued Pythagorean fuzzy numbers are used to characterize the evaluation information. The contribution of this paper is to utilize the extended Bonferroni mean operator to handle expert preferences that their preferences are interdependent, and apply an improved TODIM method to deal with the failure modes ranking, which considers the expert’s psychological behavior.
The rest of this paper is organized as follows. The relevant literature is reviewed in Section 2. The basic concepts and operational laws of an interval-valued Pythagorean fuzzy set are drawn in Section 3. The presented integrated FMEA risk assessment model is described in Section 4. A practical example is presented in Section 5 to illustrate the effectiveness and feasibility of the proposed method. Finally, Section 6 presents some conclusions and directions for future work.

2. Literature Review

2.1. Aggregation Expert Preferences

To form a collective assessment in FMEA, various information fusion operators have been utilized to aggregate expert preferences. For example, an arithmetic averaging operator, which assumes that the importance of each FMEA team expert is equal, is employed by many researchers to aggregate expert preferences [26,27,28]. In practice, FMEA team experts have different importance as they come from different departments and possess different backgrounds and expertise. Although the arithmetic averaging operator is extremely simple, it neglects the difference of expert weighting.
In order to overcome this problem, a weighted averaging (WA) operator and its extended form has been used by a number of researchers to construct a comprehensive assessment matrix in FMEA risk assessment [10,11,17,19]. On the other hand, Wang et al. [18] developed a hybrid MCDM model for evaluating the risk of failure modes, which applies an order weighted averaging (OWA) operator to fuse individual preferences under the interval-valued intuitionistic fuzzy context.
As is known, the main advantages of the WA and OWA operators consider the importance of the input argument itself and the location, respectively. To reflect the characteristics of the WA and OWA operators at the same time, Liu et al. [16] employed an interval 2-tuple hybrid averaging operator to construct a collective evaluation matrix in FMEA risk assessment. However, the above aggregation operators are only suitable to situations that input arguments that are independent, and they fail to reflect the interdependent relationships between expert preferences.

2.2. Bonferroni Mean

The Bonferroni mean (BM), initially developed by Bonferroni [29], is an effective tool for the aggregation of information of correlative arguments. The main advantage of the BM operator is its ability to capture the interrelationships among input arguments [30]. In practice, team experts may come from different departments or industries. Their subjective preferences usually indicate some interactive characteristic because these preferences are often influenced by their social status, power, knowledge, and other factors. Therefore, many studies in recent years have focused on the Bonferroni mean operator due to its ability to deal with interdependent input arguments. For example, Liu et al. [31] presented some BM operators of multi-valued neutrosophic numbers. Liu et al. [32] defined some new operational laws for an interval-valued 2-tuple linguistic and developed some new interval-valued 2-tuple linguistic BM operators based on these laws. He et al. [33] proposed two interval-valued hesitant fuzzy power BM operators by combining the BM with the power average operator. Fan et al. [34] presented linguistic neutrosophic numbers (LNNs) and a normalized weighted BM operator by merging the LNN and BM operator. Besides, the BM operator has also been extended to other fuzzy environment to aggregate various fuzzy information [35,36], such as hesitant fuzzy sets and linguistic intuitionistic fuzzy numbers.

2.3. Failure Mode Ranking

The risk priority of failure modes is a typical MCDM problem, so it is necessary to adopt an MCDM method to determine the ranking of failure modes [37]. In order to obtain a reasonable result, various MCDM methods have been applied in FMEA risk evaluation. These methods can be divided into two categories, including failure modes independence and failure modes interdependence.
For the first category, Franceschini and Galetto [14] developed an MCDM technique to determine the risk priority of failure modes in FMEA without requiring arbitrary and artificial numerical conversion. Since then, various MCDM methods have been introduced in traditional FMEA to improve its capability. For instance, Braglia et al. [9] proposed an MCDM method based on the TOPSIS (Technique for Order Preference by Similarity to an Ideal Solution) method, which obtains the prioritization of failure modes according to the descending order of the close coefficient. Chin et al. [10] presented a new FMEA model using the evidential reasoning method, which can capture the diversity opinions of FMEA team experts under different types of uncertainties well. Liu et al. [38] extended the VIKOR (VIsekriterijumska optimizacija i KOm-promisno Resenje) method, which employs the concepts of “maximum group utility” and “minimum individual regret” to determine the ranking of failure modes, thereby obtaining a compromise solution that can be accepted by most experts. Liu et al. [16] introduced the elimination and choice expressing reality (ELECTRE) to improve the ranking accuracy of failure modes in FMEA. The main characteristics of this method are that it is not compensative and it can deal with the uncertainty and fuzziness that affects the performance evaluation. Wang et al. [18] presented a new risk priority model by combing complex prorportational assessment (COPRAS) and an analytic network process to obtain the ranking of failure modes in FMEA. Huang et al. [11] developed a new FMEA model by integrating linguistic distribution assessments and a modified TODIM (an acronym in Portuguese for the interactive and MCDM) method.
On the other hand, many methods have been proposed to reflect the interdependent relationships among failure modes in the risk priority of failure modes. For instance, Xu et al. [39] developed a fuzzy-logic-based FMEA technique to solve interdependencies among various failure modes. Seyed-Hosseini et al. [40] first introduced the decision-making trial and evaluation laboratory (DEMATEL) method, which can deal with the interdependent relationships among failure modes well to rank the risk priority of failure modes. Since then, some new methods based on the DEMATEL, which is integrated with other methods, have been developed to rank the risk of failure modes. For example, Liu et al. [41] proposed a new method for the risk priority of failure modes by combining the fuzzy weighted averaging with the DEMATEL method. Motivated by the advantages of the TOPSIS method, Chang et al. [6] presented a new method that integrates TOPSIS and DEMATEL to determine the risk priority of failure modes. Liu et al. [42] combined the VIKOR, DEMATEL, and AHP method to analyze the risk of failure modes, thereby improving the assessment capability of the conventional FMEA.

2.4. TODIM Method

The TODIM method was originally proposed by Gomes and Lima [43], which is the first MCDM method based on prospect theory [24]. The prominent characteristic of the TODIM method is that it can determine the gain and loss from the reference point and also reflect the fact that the decision maker is more sensitive to loss compared to gain. Therefore, it was demonstrated to be an effective and powerful tool for handing the MCDM problems and considers the expert’s psychological behavior. Recently, the TODIM method has been widely applied to obtain the optimal ranking of alternatives in various MCDM problems. For instance, Qin et al. [44] developed an extended TODIM method based on interval type-2 fuzzy sets and applied it to deal with the green supplier selection problem. Zhang and Xu [45] presented a hesitant fuzzy TODIM method, which can capture the influences of experts’ psychological factors on evaluation results, to evaluate the sustainable water management efficiency. Furthermore, the TODIM method has also been extended to other fuzzy environments to solve MCDM problems [46,47,48], such as the multi-valued neutrosophic set, multi granulatrity linguistic terms set, and Pythagorean fuzzy set.
The literature review above indicate that (1) some aggregation operators are used to fuse the expert preferences, but no or little attention has been paid to the interdependent relationship of preferences, and (2) various MCDM methods are applied to obtain the risk priority of failure modes, but only Huang et al. [11] focused on the expert’s psychological behavior. However, Huang et al. [11] failed to consider the interdependent relationship between expert preferences and neglected the correlation between failure modes. Therefore, this study presents a risk priority model based on the extended Bonferroni mean and improved TODIM method, where the risk priority of failure modes can be determined. In this model, on the one hand, we used the extended Bonferroni mean operator to aggregate the expert preferences; on the other hand, we used an improved TODIM method to obtain the ranking of failure modes. Besides, a linguistic variable expressed by interval-valued Pythagorean fuzzy numbers is employed to represent the uncertainty and fuzziness of the evaluation information of experts.

3. Preliminaries

In this section, some basic concepts related to interval-valued Pythagorean fuzzy sets (IVPFS) are introduced briefly as follows.
Definition 1.
[49] Let A = { a 1 , a 2 , , a n } be an ordinary finite nonempty set and an IVPFS in A is an expression, P, shown as follows:
P = { < a , μ p ( a ) , ν p ( a ) > , a A } ,
where μ p ( a ) = [ μ _ p ( a ) , μ ¯ p ( a ) ] , μ p ( a ) [ 0 , 1 ] and ν p ( a ) = [ ν _ p ( a ) , ν ¯ p ( a ) ] , ν p ( a ) [ 0 , 1 ] are interval values, and μ _ p ( a ) [ 0 , 1 ] , μ ¯ p ( a ) [ 0 , 1 ] , ν _ p ( a ) [ 0 , 1 ] and ν ¯ p ( a ) [ 0 , 1 ] , satisfying μ _ p 2 ( a ) + ν ¯ p 2 ( a ) 1 .
For every a A , π _ p ( a ) = 1 μ ¯ p 2 ( a ) ν ¯ p 2 ( a ) and π ¯ p ( a ) = 1 μ _ p 2 ( a ) ν _ p 2 ( a ) . For convenience, Peng and Yang [50] designated P = ( [ μ _ p , μ ¯ p ] , [ ν _ p , ν ¯ p ] ) as an interval-valued Pythagorean fuzzy number (IVPFN), where μ _ p 2 + ν ¯ p 2 1 . Notably, the IVPFS reduces into PFS when μ _ p = μ ¯ p and ν _ p = ν ¯ p .
Definition 2.
[49] Let P = ( [ μ _ P , μ ¯ P ] , [ ν _ P , ν ¯ P ] ) , P 1 = ( [ μ _ P 1 , μ ¯ P 1 ] , [ ν _ P 1 , ν ¯ P 1 ] ) , and P 2 = ( [ μ _ P 2 , μ ¯ P 2 ] , [ ν _ P 2 , ν ¯ P 2 ] ) be three IVPFNs and λ > 0 . Then, the basic operational rules can be expressed as follows:
(1)
P 1 P 2 = ( [ μ _ p 1 2 + μ _ p 2 2 μ _ p 1 2 μ _ p 2 2 , μ ¯ p 1 2 + μ ¯ p 2 2 μ ¯ p 1 2 μ ¯ p 2 2 ] , [ ν _ p 1 ν _ p 2 , ν ¯ p 1 ν ¯ p 2 ] ) ;
(2)
P 1 P 2 = ( [ ν _ p 1 ν _ p 2 , ν ¯ p 1 ν ¯ p 2 ] , [ μ _ p 1 2 + μ _ p 2 2 μ _ p 1 2 μ _ p 2 2 , μ ¯ p 1 2 + μ ¯ p 2 2 μ ¯ p 1 2 μ ¯ p 2 2 ] ) ;
(3)
λ P = ( [ 1 ( 1 μ _ p 2 ) λ , 1 ( 1 μ ¯ p 2 ) λ ] , [ ( ν _ p ) λ , ( ν ¯ p ) λ ] ) ;
(4)
P λ = ( [ ( ν _ p ) λ , ( ν ¯ p ) λ ] , [ 1 ( 1 μ _ p 2 ) λ , 1 ( 1 μ ¯ p 2 ) λ ] ) .
Definition 3.
[49] Let P 1 = ( [ μ _ P 1 , μ ¯ P 1 ] , [ ν _ P 1 , ν ¯ P 1 ] ) and P 2 = ( [ μ _ P 2 , μ ¯ P 2 ] , [ ν _ P 2 , ν ¯ P 2 ] ) , which are two IVPFNs; then, the distance between them is defined as follows:
d ( P 1 , P 2 ) = 1 4 ( | μ _ p 1 2 μ _ p 2 2 | + | μ ¯ p 1 2 μ ¯ p 2 2 | + | ν _ p 1 2 ν _ p 2 2 | + | ν ¯ p 1 2 ν ¯ p 2 2 | + | π _ p 1 2 π _ p 2 2 | + | π ¯ p 1 2 π ¯ p 2 2 | ) .
Definition 4.
[50] Let P = ( [ μ _ P , μ ¯ P ] , [ ν _ P , ν ¯ P ] ) be an IVPFN. Then, the score and accuracy functions of P is defined respectively as follows:
s ( P ) = 1 2 ( μ _ P 2 + μ ¯ P 2 ν _ P 2 ν ¯ P 2 ) ,
a ( P ) = 1 2 ( μ _ P 2 + μ ¯ P 2 + ν _ P 2 + ν ¯ P 2 ) .
According to the above two equations, a ranking method for two IVPFNs can be introduced as follows:
(1)
If s ( P 1 ) > s ( P 2 ) , then P 1 P 2 ;
(2)
If s ( P 1 ) = s ( P 2 ) , then:
(a)
If a ( P 1 ) > a ( P 2 ) , then P 1 P 2 ;
(b)
If a ( P 1 ) = a ( P 2 ) , then P 1 = P 2 .

4. The Proposed FMEA Method

Generally, we assume that an FMEA team consists of L cross-functional experts, E k ( k = 1 , 2 , , L ) , that come from different departments and domains, and the FMEA team is responsible for the assessment of a finite set of m potential failure modes, F M i ( i = 1 , 2 , , m ) (Failure modes), with respect to a finite set of n risk factors, R F j ( j = 1 , 2 , , n ) (Risk factors). In practical risk analysis, the relative weights of FMEA team experts may not be equal because of their differing domain knowledge and experiences. Therefore, let λ = ( λ 1 , λ 2 , , λ L ) T be the weight vector of the team experts, where λ k [ 0 , 1 ] and k = 1 L λ k = 1 . Furthermore, suppose that L k = ( l i j k ) m × n and Z k = ( z j k ) 1 × n express the linguistic assessment matrix of the kth team expert, where l i j k and z j k are the linguistic ratings obtained from E k and express the assessment of F M i with respect to R F j and the relative importance of the risk factors, respectively.
Based on the assumptions and notations above, we developed an FMEA risk assessment model under the IVPFS environment. A flowchart of the proposed model is shown in Figure 1.
The main steps of the proposed method are explained as follows.
Step 1. Evaluate the failure modes’ and risk factors’ relative weights using linguistic terms.
In reality, FMEA team experts generally come from different departments and may have different practical experiences, knowledge structures, and evaluation levels [17]; the evaluation information provided by them may be certain or uncertain, complete or incomplete, and precise or imprecise [10]. Interval-valued Pythagorean fuzzy sets are a good tool to handle the imprecise and ambiguous information, and manage complex uncertainty in real world applications [49,51]. It is highly suitable for descripting the uncertainty and vagueness of experts’ judgments in FMEA. Therefore, the ratings of failure modes and the relative importance of risk factors were evaluated by linguistic variables expressed in IVPFNs in Table 1 and Table 2, respectively. Then, the linguistic evaluation matrices, L k and Z k , can be transformed into the assessment matrices, P k = ( p i j k ) m × n and W k = ( w j k ) 1 × n , respectively, according to Table 1 and Table 2, where p i j k and w j k are IVPFNs.
Step 2. Aggregate the preferences of the FMEA team experts into a group decision evaluation.
To consider the interdependent relationships between expert preferences, it is necessary to adopt a suitable aggregation operator to fuse expert preferences. The Bonferroni mean is a well-known operator because it is capable of handling interdependent input arguments. However, the conventional Bonferroni mean does not deal with IVPFN information. Hence, we extended the Bonferroni mean to the IVPFS environment to define a weighted interval-valued Pythagorean fuzzy BM operator, which is similar to the weighted interval-valued intuitionistic fuzzy BM operator [52].
Definition 5.
Let x , y 0 , where x and y do not simultaneously take the value of 0. Let p i = ( [ μ _ P i , μ ¯ P i ] , [ ν _ P i , ν ¯ P i ] ) ( i = 1 , 2 , , n ) be a collection of IVPFNs, and w = ( w 1 , w 2 , , w n ) T be the weight vector of p i , where w i indicates the importance degree of p i and satisfies w i [ 0 , 1 ] and i = 1 n w i = 1 . If:
WIVPFBM w x , y ( p 1 , p 2 , , p n ) = i , j = 1 , i j n ( 1 n ( n 1 ) × ( ( w i p i ) x ( w j p j ) y ) ) 1 x + y ,
then WIVPFBM w x , y is regarded as a WIVPFBM (Weighted interval-value Pythagorean fuzzy Bonferroni mean) operator.
Theorem 1.
Let x , y 0 , where x and y do not simultaneously take the value of 0. Let p i = ( [ μ _ P i , μ ¯ P i ] , [ ν _ P i , ν ¯ P i ] ) ( i = 1 , 2 , , n ) be a collection of IVPFNs, and w = ( w 1 , w 2 , , w n ) T be the weight vector of p i , where w i indicates the importance degree of p i and satisfies w i [ 0 , 1 ] and i = 1 n w i = 1 . Then, the aggregated value using the WIVPFBM operator is also an IVPFN, and:
WIVPFBM w x , y ( p 1 , p 2 , , p n ) = ( [ ( 1 i , j = 1 i j n ( 1 ( 1 ( 1 μ _ P i 2 ) w i ) x ( 1 ( 1 μ _ P j 2 ) w j ) y ) 1 n ( n 1 ) ) 1 x + y , ( 1 i , j = 1 i j n ( 1 ( 1 ( 1 μ ¯ P i 2 ) w i ) x ( 1 ( 1 μ ¯ P j 2 ) w j ) y ) 1 n ( n 1 ) ) 1 x + y ] , [ 1 ( 1 i , j = 1 i j n ( 1 ( 1 ν _ P i 2 w i ) x ( 1 ν _ P j 2 w j ) y ) 1 n ( n 1 ) ) 1 x + y , 1 ( 1 i , j = 1 i j n ( 1 ( 1 ν ¯ P i 2 w i ) x ( 1 ν ¯ P j 2 w j ) y ) 1 n ( n 1 ) ) 1 x + y ] ) .
According to the matrix, P k = ( p i j k ) m × n , all individual assessments can be aggregated into a group assessment matrix, P = ( p i j ) m × n , by the WIVPFBM operator, which can depict the interdependent relationships of expert preferences.
P = WIVPFBM λ x , y ( p i j 1 , p i j 2 , , p i j L ) .
Similarly, a group assessment matrix, W = ( w j ) 1 × n , for the weights of risk factors can be calculated as follows:
W = WIVPFBM λ x , y ( w i j 1 , w i j 2 , , w i j L ) .
Step 3. Determine the combination weights of risk factors.
To sufficiently reflect experts’ expertise and intrinsic information, this paper adopts subjective and objective weighting methods to determine the weights of risk factors. Numerous researchers have focused on the entropy method to obtain the objective weights of criteria [53] as it can handle the situation where conflict exists between the weight values, and eliminate the subjectivity of experts’ judgments in the valuation of weights. Therefore, we defined an interval-valued Pythagorean fuzzy entropy, which is parallel to interval-valued intuitionistic fuzzy entropy [54], to calculate the objective weight of risk factors.
Definition 6.
Let P be an IVPFS defined in the universe of discourse A. An interval-valued Pythagorean fuzzy entropy (IVPFE) is given as follows:
E ( P ) = 1 n i = 1 n 2 | μ _ P i 2 ν _ P i 2 | | μ ¯ P i 2 ν ¯ P i 2 | + π _ P i 2 + π ¯ P i 2 2 + | μ _ P i 2 ν _ P i 2 | + | μ ¯ P i 2 ν ¯ P i 2 | + π _ P i 2 + π ¯ P i 2 .
Step 3.1. Obtain the subjective weights of risk factors.
Based on the aggregated weights, W = ( w j ) 1 × n ( j = 1 , 2 , , n ) , of risk factors, the score values of weights of each risk factor can be calculated as follows:
s ( w j ) = 1 2 ( μ _ P i 2 + μ ¯ P i 2 ν _ P i 2 ν ¯ P i 2 ) .
The normalized subjective weight of each risk factor can be obtained as:
w j s = s ( w j ) j = 1 n s ( w j ) .
Step 3.2. Obtain the objective weights of risk factors using the entropy method.
The objective weights, w j o ( j = 1 , 2 , , n ) , of risk factors can be determined through the entropy method. The entropy value of each failure mode with respect to the risk factor is calculated as follows:
E i j = 2 | μ _ P i 2 ν _ P i 2 | | μ ¯ P i 2 ν ¯ P i 2 | + π _ P i 2 + π ¯ P i 2 2 + | μ _ P i 2 ν _ P i 2 | + | μ ¯ P i 2 ν ¯ P i 2 | + π _ P i 2 + π ¯ P i 2 .
Thus, the normalized objective weight of each risk factor can be obtained by the following equation:
w j o = i = 1 m ( 1 E i j ) j = 1 n i = 1 m ( 1 E i j ) .
Step 3.3. Calculate the comprehensive weights of risk factors.
It is significant to combine the subjective and objective weights to calculate the comprehensive weight of each risk factor, which can consider simultaneously the importance of risk factors and the experiences of the experts [28]. Therefore, based on the w j s and w j o , the comprehensive weight is obtained as follows:
w ^ j = ζ w j s + ( 1 ζ ) w j o , ( j = 1 , 2 , , n ) ,
where parameter ζ [ 0 , 1 ] is the relative importance coefficient for the subjective weight. In this study, two types of weights were assumed to be equally important, that is, ζ = 0.5 .
Step 4. Obtain the ranking of failure modes by using the improved TODIM method.
TODIM, which is an effective behavior MCDM method based on prospect theory, was used by Huang et al. [11] to rank the risk priority of failure modes, but it cannot consider the interrelationships between failure modes when calculating the dominance degree of failure modes. In the real world, interdependent relationships may exist between failure modes because of the complexity of the failure system. Grey relational analysis is part of a grey system, which is suitable for handling problems with complicated interrelationships between multiple factors and variables [55]. Therefore, grey relational analysis is introduced in the TODIM method to compute the dominance degree of failure modes.
Step 4.1. Calculate the dominance degree of failure modes.
The dominance degree of F M i over F M s regarding the R F j can be obtained as follows:
Φ j ( F M i , F M s ) = { ( ξ ( p i j , p 0 ) ξ ( p s j , p 0 ) ) w ^ i j j = 1 n w ^ i j , i f   s ( p i j ) > s ( p s j ) 0 , i f   s ( p i j ) = s ( p s j ) 1 θ ( ξ ( p i j , p 0 ) ξ ( p s j , p 0 ) ) j = 1 n w ^ i j w ^ i j , i f   s ( p i j ) < s ( p s j ) ,
where Φ j ( F M i , F M s ) represents the contribution of the risk factor, R F j , to the function, Ω j ( F M i , F M s ) , when comparing F M s with F M i . If s ( p i j ) > s ( p s j ) , then a gain occurs. If s ( p i j ) < s ( p s j ) , then a loss is incurred. Finally, if s ( p i j ) = s ( p s j ) , then the result is zero. Moreover, w ^ i j = w ^ j / w ^ ( j = 1 , 2 , , n ) denotes the relative weight of R F j to the reference, R F l , and w ^ = max { w ^ j | j = 1 , 2 , , n } . Parameter θ represents the attenuation coefficient of the loss, which can be tuned according to the practical problem. In this study, the value, θ , was assumed to be 0.45. The term, ξ ( p i j , p 0 ) , denotes the grey relation coefficient between the group evaluation value and the standard value, which is calculated as follows:
ξ ( p i j , p 0 ) = min min | s ( p i j ) s ( p 0 ) | + δ max max | s ( p i j ) s ( p 0 ) | | s ( p i j ) s ( p 0 ) | + δ max max | s ( p i j ) s ( p 0 ) | .
Step 4.2. Obtain the overall dominance degree of failure modes.
The overall dominance degree of F M i over F M s can be calculated according to the following expression:
Ω ( F M i , F M s ) = j = 1 n Φ j ( F M i , F M s ) , ( i , s = 1 , 2 , , m ) .
Step 4.3. Calculate the global dominance degree of failure modes.
The global dominance degree of F M i can be obtained as:
ψ ( F M i ) = s = 1 m Ω ( F M i , F M s ) min i m { s = 1 m Ω ( F M i , F M s ) } max i m { s = 1 m Ω ( F M i , F M s ) } min i m { s = 1 m Ω ( F M i , F M s ) } .
Step 4.4. Rank the failure modes according to ψ ( F M i ) .

5. Case Illustration

A nuclear reheat valve system in a nuclear steam turbine [28,56] was selected as the case study to demonstrate the applicability of the proposed method. This valve system is crucial to the operation of nuclear power stations, and its failure would decrease the reliability of the entire power station. Therefore, the reheat valve system must be quickly closed during abnormal operating conditions to cut off the steam entering the low-pressure cylinder and ensure system safety and reliability.
An FMEA team with four experts was formed to identify the most significant failure modes. The four team experts were given the following relative weights as 0.15, 0.30, 0.35, and 0.20 because of their different background knowledge and experiences. The team identified eight failure modes through brainstorming and selected three risk factors, namely, S, O, and D. The eight failure modes were a long closing time or no valve action, valve nor closing tightly, large leak around valve shaft, valve fluctuations, valve jam during operation, valve shaft fracture, fault of valve shaft support bearing, and abnormal noise from valve system.

5.1. Implementation of the Proposed Method

In this section, the proposed method was applied to solve the valve system risk assessment problem, and the implementation procedure is presented as follows.
Step 1. Evaluate the failure modes and risk factors’ relative weights using linguistic terms.
The evaluations of the eight failure modes and the risk factors’ relative weights provided by experts are shown Table 3 and Table 4, respectively. Subsequently, the linguistic ratings were transformed into the corresponding IVPFNs.
Step 2. Establish a group decision matrix by utilizing the WIVPFBM operator.
The aggregated fuzzy ratings of each failure mode were calculated by the WIVPFBM operator to establish the fuzzy collective evaluation matrix, P = ( p i j ) 8 × 3 , as shown in Table 5. Subsequently, the aggregated fuzzy weights vector of risk factors was calculated by Equation (8), as shown in Table 5.
Step 3. Determine the comprehensive weights of risk factors.
According to Equations (10)–(13), the subjective and objective weights of risk factors were obtained as w s = ( 0.3935 , 0.2767 , 0.3298 ) T and w o = ( 0.2563 , 0.3635 , 0.3802 ) T , respectively. Then, the comprehensive weights of risk factors were determined as w = ( 0.3249 , 0.3201 , 0.3550 ) T .
Step 4. Use the improved TODIM method to rank the risk priority of failure modes.
Step 4.1. Calculate the dominance degree of the failure modes.
Applying Equations (15) (i.e., θ = 0.45 ) and (16) (i.e., δ = 0.5 ), the dominance degree matrices concerning the risk factor, S, O, and D, were obtained, respectively.
Φ S = ( 0 3.0274 2.5009 0.2140 0.8967 0.6298 1.2359 1.9135 0.4426 0 0.2494 0.3131 0.4228 0.4329 0.4041 0.3430 0.3656 1.7061 0 0.1892 0.3413 0.3539 0.3179 0.2354 0.3129 2.1414 1.2942 0 0.2841 0.2990 0.2554 0.1401 0.1311 2.8916 2.3346 1.9431 0 0.0933 0.8506 1.6904 0.0921 2.9612 2.4203 2.0452 0.6382 0 1.0634 1.8069 0.1807 2.7637 2.1742 1.7470 0.1244 0.1555 0 1.4608 0.2798 2.3460 1.6103 0.9582 0.2471 0.2642 0.2136 0 ) ;
Φ O = ( 0 0.9524 0.2775 0.2016 1.2758 0.2250 0.2558 0.2815 0.1372 0 0.3095 0.2439 0.8488 0.2635 0.2903 0.3131 1.9262 2.1488 0 1.3232 2.3104 1.1275 0.7463 0.0473 1.3997 1.6930 0.1906 0 1.8939 0.0998 0.1574 0.1964 0.1838 0.1223 0.3328 0.2728 0 0.2905 0.3150 0.3361 1.5617 1.8292 0.1642 0.6927 2.0166 0 0.1217 0.1692 1.7758 2.0150 0.1075 1.0927 2.1865 0.8451 0 0.1175 1.9540 2.1738 0.3286 1.3634 2.3336 1.1744 0.8154 0 ) ;
Φ D = ( 0 2.4641 2.1292 2.0989 1.8226 2.2266 1.8351 2.4234 0.3936 0 0.1981 0.2062 0.2649 0.1686 0.2627 0.0712 0.3401 1.2404 0 0.0571 0.1758 0.6514 0.1725 1.1574 0.3353 1.2909 0.3577 0 0.1663 0.7431 0.1627 1.2114 0.2912 1.6584 1.1007 1.0410 0 1.2790 0.2143 1.5972 0.3557 1.0556 0.1041 0.1187 0.2043 0 0.2014 0.9567 0.2932 1.6445 1.0797 1.0187 0.0342 1.2609 0 1.5828 0.3871 0.4460 0.1849 0.1935 0.2552 0.1528 0.2529 0 ) .
Step 4.2. Compute the overall dominance degree of the failure modes.
The overall dominance degree matrix can be obtained based on Equation (17).
Ω = ( 0 6.4440 4.3527 4.0373 3.9950 2.6315 2.8153 4.0555 0.9735 0 0.7571 0.7632 0.1611 0.8651 0.9570 0.7274 1.2204 5.0952 0 1.0769 1.7932 1.4250 0.2559 0.8746 0.7515 5.1254 1.4613 0 1.4435 0.3443 0.5756 0.8749 0.6060 4.4277 3.1026 2.7113 0 0.8952 0.7499 2.9515 1.1140 5.8460 2.1538 2.6192 2.4505 0 0.7402 2.5944 1.3019 6.4232 3.1464 3.8584 2.0279 1.9506 0 2.9262 1.2871 0.9658 1.7540 2.1281 1.8313 0.7574 0.3490 0 ) .
Step 4.3. Determine the global dominance degree of each failure mode.
The global dominance of the failure mode, F M i , was computed by Equation (18), namely, ψ ( F M 1 ) = 0.9795 , ψ ( F M 2 ) = 0 , ψ ( F M 3 ) = 0.6613 , ψ ( F M 4 ) = 0.6483 , ψ ( F M 5 ) = 0.7046 , ψ ( F M 6 ) = 0.8924 , ψ ( F M 7 ) = 1 , and ψ ( F M 8 ) = 0.7089 .
Step 4.4. Rank the failure modes.
Based on the global dominance degree of failure modes, the risk priority of failure modes is F M 7 F M 1 F M 6 F M 8 F M 5 F M 3 F M 4 F M 2 . According to the results, F M 7 is the most serious failure mode and should be given the top risk priority, followed by F M 1 , F M 6 , F M 8 , F M 5 , F M 3 , F M 4 , and F M 2 .

5.2. Sensitivity Analysis

A sensitivity analysis on the expert preferences’ interdependence and failure modes’ interdependence was conducted to illustrate their impact on the risk priority of failure modes. In the sensitivity analysis, we considered four methods as follows: (1) The proposed method in this paper (case 0); (2) GRA (Grey relational analysis)is not employed to compute the dominance degree of failure modes (case 1); (3) an interval-valued Pythagorean fuzzy weighted averaging (IVPFWA) operator is utilized to aggregate the expert preferences in Step 2 (case 2), and (4) the IVPFWA operator is used to aggregate the expert preferences and GRA is not applied to calculate the dominance degree of failure modes (case 3). The ranking results of the eight failure modes under the four cases are shown in Figure 2.
As we can see in Figure 2, the risk priority of failure modes in case 0 is very different from the ranking results in the other three cases. This divergence can be demonstrated by the Spearman rank-correlation coefficients between the rankings of failure modes in case 0 and case 1, case 2, and case 3, and the correlation coefficients were 0.7143, 0.5714, and 0.6429, respectively. The ranking of failure modes in case 0 is different from those in case 1 except for F M 3 . The largest difference between the two ranking results happens at F M 7 , and F M 7 has a difference of three ranking places. The main reason for this difference is that the failure modes’ interdependence was not considered in the calculation process of the dominance degree of failure modes in case 1. There is some difference between the ranking results in case 0 and case 2. In these two ranking results of the failure modes, F M 2 and F M 6 have a difference of four ranking places and three ranking places, respectively. The reason for this difference is that the IVPFWA operator was applied in case 2, rather than the WIVPFBM operator, to aggregate expert preferences, which neglects the interdependent relationship between the expert preferences. In addition, it can be clearly seen from Figure 1 that the eight failure modes have different risk priorities in case 0 and case 3. In case 3, the IVPFWA operator was utilized to aggregate expert preferences and GRA was not employed to compute the dominance degree of failure modes, which may lead to this difference. The sensitivity analysis shows that the expert preferences’ interdependence and the failure modes’ interdependence have a great influence on the final risk priority of failure modes.

5.3. Comparisons and Discussion

To further demonstrate the effectiveness of the proposed method, comparisons of the results with the fuzzy TOPSIS [28], fuzzy VIKOR [12], and fuzzy digraph and matrix (FDM) [56] were made. Figure 3 shows the ranking results of all eight failure modes obtained by the four methods.
The rankings of failure modes obtained by the proposed method and fuzzy TOPSIS method show a notable difference. The ranking of F M 7 differs significantly for the fuzzy TOPSIS and the proposed method. F M 7 only ranked in fourth place in the former method, while it ranked first place in the latter method. Rankings determined by the fuzzy TOPSIS method may be irrational because each expert was assigned the same weight. In reality, experts usually come from several organizations and have different backgrounds, knowledge structures, and practical experience. Hence, they should be given different weights during risk assessments. In addition, the experts’ psychological behavior, expert preferences’ interdependence, and failure modes’ interdependence were not considered in the fuzzy TOPSIS.
The ranking determined by the fuzzy VIKOR significantly differed from the ranking obtained by the proposed method. The Spearman rank-correlation coefficient between the two rankings of failure modes was 0.4286, which further validated this point. For example, F M 5 has a difference of four ranking places, and F M 2 , F M 6 , and F M 7 have a divergence of three ranking places. The reasons for this difference between the two ranking results are that the expert preferences’ interdependence, experts’ psychological behavior, and failure modes’ interdependence were ignored in the fuzzy VIKOR methods. Furthermore, F M 8 ranks behind F M 5 with the fuzzy VIKOR method. However, the former is more important in reality; therefore, F M 8 merits a higher priority in comparison with F M 5 .
The risk priority of failure modes determined by the FDM method are slightly different from the results yielded by the proposed method, except for F M 2 , F M 6 , and F M 4 . The inconsistent ranking results may be explained by the following: The FDM method considered the interrelations between risk factors by a risk factors fuzzy digraph, while it neglected the expert preferences’ interdependence and the experts’ psychological behavior. Furthermore, the FDM method only considered the subjective weights of risk factors, which may also cause this difference between the two ranking results of failure modes.
Based on the comparative analysis above, compared with the fuzzy TOPSIS, fuzzy VIKOR, and FDM method, the advantages of the proposed method are summarized as follows: (1) The proposed method adopts the WIVPFBM operator to fuse the evaluations of experts into a group assessment, which reflects the interdependent relationships among expert preferences; (2) the risk priority of failure modes is obtained by applying the improved TODIM method, which considers the experts’ psychological behavior and reflects the interdependence between failure modes.

6. Conclusions

In this study, we proposed a new risk assessment method of FMEA by combing the WIVPFBM operator and an improved TODIM method to evaluate and determine the risk priority of failure modes under the IVPFS environment. In this method, all evaluation information of FMEA team experts was initially provided in the form of linguistic terms represented by the IVPFNs. Then, the WIVPFBM operator was applied to aggregate the individual evaluation matrix into a group decision matrix. Subsequently, the weights of risk factors were obtained by the subjective and objective weighting method. Finally, the risk ranking of failure modes was obtained by the improved TODIM method. In addition, a sensitivity analysis indicated that the expert preferences’ interdependence and failure modes’ interdependence had a great influence on the final risk priority of failure modes.
The flexibility and effectiveness of the presented method were demonstrated by the application of the nuclear reheat valve system. The results indicate that the proposed method, which considered the team experts’ psychological behaviors and the interrelations among expert preferences, is effective and flexible for real applications compared with other improved FMEA methods. As a direction for future research, extending the proposed method by considering the consensus reaching process in the risk assessment process is recommended.

Author Contributions

Conceptualization, J.Z. and B.S.; methodology, J.Z. and B.S.; validation, R.W. and K.-S.C.; formal analysis, R.W. and K.-S.C.; data curation, J.Z. and B.S; writing—original draft preparation, J.Z.; writing—review and editing, J.Z. and K.-S.C.; visualization, R.W.; supervision, B.S.; project administration, B.S.

Funding

This research was funded by National Natural Science Foundation of China, grant number (71173177) and the Doctoral Innovation Fund Program of Southwest Jiaotong University (D-CX201729).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stamatis, D.H. Failure Mode and Effect Analysis: FMEA from Theory to Execution; ASQ Quality Press: Milwaukee, WI, USA, 2003. [Google Scholar]
  2. Bowles, J.B.; Pelaez, C.E. Fuzzy logic prioritization of failures in a system failure mode, effects and criticality analysis. Reliab. Eng. Syst. Saf. 1995, 50, 203–213. [Google Scholar] [CrossRef]
  3. Yeh, T.M.; Chen, L.Y. Fuzzy-based risk priority number in FMEA for semiconductor wafer processes. Int. J. Prod. Res. 2014, 52, 539–549. [Google Scholar] [CrossRef]
  4. Karim, M.A.; Smith, A.J.R.; Halgamuge, S. Empirical relationships between some manufacturing practices and performance. Int. J. Prod. Res. 2008, 46, 3583–3613. [Google Scholar] [CrossRef]
  5. Braaksma, A.J.J.; Meesters, A.J.; Klingenberg, W.; Hicks, C. A quantitative method for failure mode and effects analysis. Int. J. Prod. Res. 2012, 50, 6904–6917. [Google Scholar] [CrossRef]
  6. Chang, K.H.; Chang, Y.C.; Lee, Y.T. Integrating TPOSIS and dematel methods to rank the risk of failure of FMEA. Int. J. Inf. Technol. Decis. Mak. 2014, 13, 1229–1257. [Google Scholar] [CrossRef]
  7. Bozdag, E.; Asan, U.; Soyer, A.; Serdarasan, S. Risk prioritization in failure mode and effects analysis using interval type-2 fuzzy sets. Expert Syst. Appl. 2015, 42, 4000–4015. [Google Scholar] [CrossRef]
  8. Certa, A.; Hopps, F.; Inghilleri, R.; Fata, C.M.L. A Dempster-shafer theory-based approach to the failure mode, effects and criticality analysis (FMECA) under epistemic uncertainty: Application to the propulsion system of a fishing vessel. Reliab. Eng. Syst. Saf. 2017, 159, 69–79. [Google Scholar] [CrossRef]
  9. Braglia, M.; Frosolini, M.; Montanari, R. Fuzzy TOPSIS approach for failure mode, effects and criticality analysis. Qual. Reliab. Eng. Int. 2010, 19, 425–443. [Google Scholar] [CrossRef]
  10. Chin, K.S.; Wang, Y.M.; Poon, G.K.K.; Yang, J.B. Failure mode and effects analysis using a group-based evidential reasoning approach. Comput. Oper. Res. 2009, 36, 1768–1779. [Google Scholar] [CrossRef]
  11. Huang, J.; Li, Z.; Liu, H.C. New approach for failure mode and effect analysis using linguistic distribution assessments and TODIM method. Reliab. Eng. Syst. Saf. 2017, 167, 302–309. [Google Scholar] [CrossRef]
  12. Safari, H.; Faraji, Z.; Majidian, S. Identifying and evaluating enterprise architecture risks using FMEA and fuzzy vikor. J. Intell. Manuf. 2016, 27, 475–486. [Google Scholar] [CrossRef]
  13. Liu, H.C.; Liu, L.; Liu, N. Risk evaluation approaches in failure mode and effects analysis: A literature review. Expert Syst. Appl. 2013, 40, 828–838. [Google Scholar] [CrossRef]
  14. Franceschini, F.; Galetto, M. A new approach for evaluation of risk priorities of failure modes in FMEA. Int. J. Prod. Res. 2001, 39, 2991–3002. [Google Scholar] [CrossRef] [Green Version]
  15. Guo, J. A risk assessment approach for failure mode and effects analysis based on intuitionistic fuzzy sets and evidence theory. J. Intell. Fuzzy Syst. 2016, 30, 869–881. [Google Scholar] [CrossRef]
  16. Liu, H.C.; You, J.X.; Chen, S.; Chen, Y.Z. An integrated failure mode and effect analysis approach for accurate risk assessment under uncertainty. Iie Trans. 2016, 48, 1027–1042. [Google Scholar] [CrossRef]
  17. Liu, H.C.; You, J.X.; Duan, C.Y. An integrated approach for failure mode and effect analysis under interval-valued intuitionistic fuzzy environment. Int. J. Prod. Econ. 2017, 207, 163–172. [Google Scholar] [CrossRef]
  18. Wang, L.E.; Liu, H.C.; Quan, M.Y. Evaluating the risk of failure modes with a hybrid MCDM model under interval-valued intuitionistic fuzzy environment. Comput. Ind. Eng. 2016, 102, 175–185. [Google Scholar] [CrossRef]
  19. Zhao, H.; You, J.X.; Liu, H.C. Failure mode and effect analysis using multimoora method with continuous weighted entropy under interval-valued intuitionistic fuzzy environment. Soft Comput. 2017, 21, 5355–5367. [Google Scholar] [CrossRef]
  20. Kabak, Ö.; Ervural, B. Multiple attribute group decision making: A generic conceptual framework and a classification scheme. Knowl. Based Syst. 2017, 123, 13–30. [Google Scholar] [CrossRef]
  21. Ju, D.; Ju, Y.; Wang, A. Multi-attribute group decision making based on power generalized Heronian mean operator under hesitant fuzzy linguistic environment. Soft Comput. 2018, 23, 3823–3842. [Google Scholar] [CrossRef]
  22. Camerer, C. Bounded rationality in individual decision making. Exp. Econ. 1998, 1, 163–183. [Google Scholar] [CrossRef]
  23. Charness, G.; Rabin, M. Expressed preferences and behavior in experimental games. Games Econ. Behav. 2005, 53, 151–169. [Google Scholar] [CrossRef] [Green Version]
  24. Kahneman, D.; Tversky, A. Prospect theory: An analysis of decision under risk. Econometrica 1979, 47, 263–291. [Google Scholar] [CrossRef]
  25. Zhang, X.; Xu, Z. The TODIM analysis approach based on novel measured functions under hesitant fuzzy environment. Knowl. Based Syst. 2014, 61, 48–58. [Google Scholar] [CrossRef]
  26. Mohsen, O.; Fereshteh, N. An extended VIKORmethod based on entropy measure for the failure modes risk assessment – a case study of the geothermal power plant (GPP). Saf. Sci. 2017, 92, 160–172. [Google Scholar] [CrossRef]
  27. Selim, H.; Yunusoglu, M.G.; Balaman, Ş.Y. A dynamic maintenance planning framework based on fuzzy TOPSIS and FMEA: Application in an international food company. Qual. Reliab. Eng. Int. 2016, 32, 795–804. [Google Scholar] [CrossRef]
  28. Song, W.; Ming, X.; Wu, Z.; Zhu, B. Failure modes and effects analysis using integrated weight-based fuzzy TOPSIS. Int. J. Comput. Integr. Manuf. 2013, 26, 1172–1186. [Google Scholar] [CrossRef]
  29. Bonferroni, C. Sulle medie multiple di potenze. Boll. Dell’Unione Mat. Ital. 1950, 5, 267–270. [Google Scholar]
  30. Liang, D.; Xu, Z.; Darko, A.P. Projection model for fusing the information of Pythagorean fuzzy multicriteria group decision making based on geometric Bonferroni mean. Int. J. Intell. Syst. 2017, 32, 966–987. [Google Scholar] [CrossRef]
  31. Liu, P.; Zhang, L.; Liu, X.; Wang, P. Multi-valued neutrosophic number Bonferroni mean operators with their applications in multiple attribute group decision making. Int. J. Inf. Technol. Decis. Mak. 2016, 15, 1181–1210. [Google Scholar] [CrossRef]
  32. Liu, X.; Tao, Z.; Chen, H.; Zhou, L. A new interval-valued 2-tuple linguistic Bonferroni mean operator and its application to multiattribute group decision making. Int. J. Fuzzy Syst. 2017, 19, 86–108. [Google Scholar] [CrossRef]
  33. He, Y.; Zhen, H.; Shi, L.; Meng, S. Multiple attribute group decision making based on ivhfpbms and a new ranking method for interval-valued hesitant fuzzy information. Comput. Ind. Eng. 2016, 99, 63–77. [Google Scholar] [CrossRef]
  34. Fan, C.; Ye, J.; Hu, K.; Fan, E. Bonferroni mean operators of linguistic neutrosophic numbers and their multiple attribute group decision-making methods. Information 2017, 8, 107. [Google Scholar] [CrossRef]
  35. He, Y.; Zhen, H.; Wang, G.; Chen, H. Hesitant fuzzy power Bonferroni means and their application to multiple attribute decision making. IEEE Trans. Fuzzy Syst. 2015, 23, 1655–1668. [Google Scholar] [CrossRef]
  36. Liu, P.; Xi, L. Multiattribute group decision making methods based on linguistic intuitionistic fuzzy power Bonferroni mean operators. Complexity 2017, 2017, 1–15. [Google Scholar] [CrossRef]
  37. Lolli, F.; Ishizaka, A.; Gamberini, R.; Rimini, B.; Messori, M. Flowsort-gdss—A novel group multi-criteria decision support system for sorting problems with application to FMEA. Expert Syst. Appl. Int. J. 2015, 42, 6342–6349. [Google Scholar] [CrossRef]
  38. Liu, H.C.; Liu, L.; Liu, N.; Mao, L.X. Risk evaluation in failure mode and effects analysis with extended VIKOR method under fuzzy environment. Expert Syst. Appl. 2012, 39, 12926–12934. [Google Scholar] [CrossRef]
  39. Xu, K.; Tang, L.C.; Xie, M.; Ho, S.L.; Zhu, M.L. Fuzzy assessment of FMEA for engine systems. Reliab. Eng. Syst. Saf. 2002, 75, 17–29. [Google Scholar] [CrossRef]
  40. Seyed-Hosseini, S.M.; Safaei, N.; Asgharpour, M.J. Reprioritization of failures in a system failure mode and effects analysis by decision making trial and evaluation laboratory technique. Reliab. Eng. Syst. Saf. 2006, 91, 872–881. [Google Scholar] [CrossRef]
  41. Liu, H.C.; You, J.X.; Lin, Q.L.; Li, H. Risk assessment in system FMEA combining fuzzy weighted average with fuzzy decision-making trial and evaluation laboratory. Int. J. Comput. Integr. Manuf. 2015, 28, 701–714. [Google Scholar] [CrossRef]
  42. Liu, H.C.; You, J.X.; Ding, X.F.; Su, Q. Improving risk evaluation in FMEA with a hybrid multiple criteria decision making method. Int. J. Qual. Reliab. Manag. 2015, 32, 763–782. [Google Scholar] [CrossRef]
  43. Gomes, L.; Lima, M. TODIM: Basics and application to multicriteria ranking of projects with environmental impacts. Found. Comput. Decis. Sci. 1992, 16, 113–127. [Google Scholar]
  44. Qin, J.; Liu, X.; Pedrycz, W. An extended TODIM multi-criteria group decision making method for green supplier selection in interval type-2 fuzzy environment. Eur. J. Oper. Res. 2017, 258, 626–638. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Xu, Z. Efficiency evaluation of sustainable water management using the HF-TODIM method. Int. Trans. Oper. Res. 2019, 26, 747–764. [Google Scholar] [CrossRef]
  46. Ji, P.; Zhang, H.Y.; Wang, J.Q. A projection-based TODIM method under multi-valued neutrosophic environments and its application in personnel selection. Neural Comput. Appl. 2018, 29, 221–234. [Google Scholar] [CrossRef]
  47. Qin, Q.; Liang, F.; Li, L.; Wei, Y.M. Selection of energy performance contracting business models: A behavioral decision-making approach. Renew. Sustain. Energy Rev. 2017, 72, 422–433. [Google Scholar] [CrossRef]
  48. Ren, P.; Xu, Z.; Gou, X. Pythagorean fuzzy TODIM approach to multi-criteria decision making. Appl. Soft Comput. 2016, 42, 246–259. [Google Scholar] [CrossRef]
  49. Zhang, X. Multicriteria Pythagorean fuzzy decision analysis: A hierarchical QUALIFLEX approach with the closeness index-based ranking methods. Inf. Sci. 2016, 330, 104–124. [Google Scholar] [CrossRef]
  50. Peng, X.; Yang, Y. Fundamental properties of interval-valued Pythagorean fuzzy aggregation operators. Int. J. Intell. Syst. 2016, 31, 444–487. [Google Scholar] [CrossRef]
  51. Chen, T.Y. An interval-valued Pythagorean fuzzy outranking method with a closeness-based assignment model for multiple criteria decision making. Int. J. Intell. Syst. 2018, 33, 126–168. [Google Scholar] [CrossRef]
  52. Xu, Z.; Chen, Q. A multi-criteria decision making procedure based on interval-valuedintuitionistic fuzzy Bonferroni mean. J. Syst. Sci. Syst. Eng. 2011, 20, 217–228. [Google Scholar] [CrossRef]
  53. Bazzazi, A.A.; Osanloo, M.; Karimi, B. Deriving preference order of open pit mines equipment through MADM methods: Application of modified vikor method. Expert Syst. Appl. 2011, 38, 2550–2556. [Google Scholar] [CrossRef]
  54. Wei, C.P.; Wang, P.; Zhang, Y.Z. Entropy, similarity measure of interval-valued intuitionistic fuzzy sets and their applications. Inf. Sci. 2011, 181, 4273–4286. [Google Scholar] [CrossRef]
  55. Wei, G. Grey relational analysis model for dynamic hybrid multiple attribute decision making. Knowl. Based Syst. 2011, 24, 672–679. [Google Scholar] [CrossRef]
  56. Liu, H.C.; Chen, Y.Z.; You, J.X.; Li, H. Risk evaluation in failure mode and effects analysis using fuzzy digraph and matrix approach. J. Intell. Manuf. 2016, 27, 805–816. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed FMEA (Failure modes and effective analysis) model.
Figure 1. Flowchart of the proposed FMEA (Failure modes and effective analysis) model.
Mathematics 07 00536 g001
Figure 2. The ranking results of the sensitivity analysis.
Figure 2. The ranking results of the sensitivity analysis.
Mathematics 07 00536 g002
Figure 3. Comparative rankings of failure modes.
Figure 3. Comparative rankings of failure modes.
Mathematics 07 00536 g003
Table 1. Linguistic terms for rating the failure modes [50].
Table 1. Linguistic terms for rating the failure modes [50].
Linguistic VariablesAbbreviationIVPFNs
Very highVH([0.8000, 0.9000], [0.1000, 0.2000])
HighH([0.7000, 0.8000], [0.2000, 0.3000])
Medium highMH([0.6000, 0.7000], [0.3000, 0.4000])
MediumM([0.5000, 0.6000], [0.4000, 0.5000])
Medium lowML([0.3000, 0.4000], [0.6000, 0.7000])
LowL([0.2000, 0.3000], [0.7000, 0.8000])
Very lowVL([0.1000, 0.2000], [0.8000, 0.9000])
Table 2. Linguistic terms for rating the weights of risk factors [50].
Table 2. Linguistic terms for rating the weights of risk factors [50].
Linguistic VariablesAbbreviationIVPFNs
Very highVH([0.8000, 0.9000], [0.1000, 0.2000])
HighH([0.7000, 0.8000], [0.2000, 0.3000])
MediumM([0.5000, 0.6000], [0.4000, 0.5000])
LowL([0.3000, 0.4000], [0.6000, 0.7000])
Very lowVL([0.1000, 0.2000], [0.8000, 0.9000])
Table 3. Assessment information on the eight failure modes by FMEA team experts.
Table 3. Assessment information on the eight failure modes by FMEA team experts.
Risk FactorsSeverity (S)Occurrence (O)Detection (D)
Team ExpertsE1E2E3E4E1E2E3E4E1E2E3E4
FM1VHVHVHHMLMLMLHHMHM
FM2LVLVHMLMLMLMMLMLMMML
FM3MLMLLMLMHMHMMMLMLLM
FM4MHMHMMMLLMLMLMLMLML
FM5VHVHMHHMLMLMMMLMLLML
FM6VHVHMHVHLMLLMLLMLMLM
FM7MHMHVHHMMMHMMLLMLML
FM8MMHMHMMHMMHMLMMLML
Table 4. Linguistic evaluations of risk factor weights.
Table 4. Linguistic evaluations of risk factor weights.
Team ExpertsE1E2E3E4
OccurrenceLMMH
SeverityMHHVH
DetectionMHMH
Table 5. Fuzzy collective evaluation matrix and aggregated fuzzy weights of risk factors.
Table 5. Fuzzy collective evaluation matrix and aggregated fuzzy weights of risk factors.
Risk FactorsSOD
FM1([0.4517,0.5520],
[0.6053,0.6972])
([0.1743,0.2292],
[0.8708,0.9062])
([0.3416,0.4138],
[0.7297,0.7875])
FM2([0.1943, 0.2609],
[0.8627, 0.9045])
([0.1847, 0.2403],
[0.8620, 0.8982])
([0.2200, 0.2770],
[0.8384, 0.8772])
FM3([0.1357, 0.1884],
[0.8936, 0.9264])
([0.2898, 0.3540],
[0.7764, 0.8248])
([0.1619, 0.2163],
[0.8768, 0.9123])
FM4([0.2898, 0.3540],
[0.7764, 0.8248])
([0.1371, 0.1897],
[0.8931, 0.9258])
([0.1597, 0.2136],
[0.8785, 0.9136])
FM5([0.4077, 0.4978],
[0.6539, 0.7330])
([0.2120, 0.2689],
[0.8426, 0.8813])
([0.1357, 0.1884],
[0.8936, 0.9264])
FM6([0.4264, 0.5231],
[0.6286, 0.7165])
([0.1268, 0.1791],
[0.9000, 0.9321])
([0.1698, 0.2243],
[0.8722, 0.9077])
FM7([0.3818, 0.4630],
[0.6928, 0.7586])
([0.2796, 0.3422],
[0.7872, 0.8332])
([0.1371, 0.1897],
[0.8931, 0.9258])
FM8([0.3141, 0.3818],
[0.7586, 0.8099])
([0.2924, 0.3568],
[0.7745, 0.8231])
([0.1829, 0.2385],
[0.8625, 0.8988])
Weights([0.3736, 0.4573],
[0.6817, 0.7526])
([0.2677, 0.3366],
[0.7913, 0.8386])
([0.3177, 0.3873],
[0.7421, 0.7980])

Share and Cite

MDPI and ACS Style

Zhu, J.; Shuai, B.; Wang, R.; Chin, K.-S. Risk Assessment for Failure Mode and Effects Analysis Using the Bonferroni Mean and TODIM Method. Mathematics 2019, 7, 536. https://doi.org/10.3390/math7060536

AMA Style

Zhu J, Shuai B, Wang R, Chin K-S. Risk Assessment for Failure Mode and Effects Analysis Using the Bonferroni Mean and TODIM Method. Mathematics. 2019; 7(6):536. https://doi.org/10.3390/math7060536

Chicago/Turabian Style

Zhu, Jianghong, Bin Shuai, Rui Wang, and Kwai-Sang Chin. 2019. "Risk Assessment for Failure Mode and Effects Analysis Using the Bonferroni Mean and TODIM Method" Mathematics 7, no. 6: 536. https://doi.org/10.3390/math7060536

APA Style

Zhu, J., Shuai, B., Wang, R., & Chin, K. -S. (2019). Risk Assessment for Failure Mode and Effects Analysis Using the Bonferroni Mean and TODIM Method. Mathematics, 7(6), 536. https://doi.org/10.3390/math7060536

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop