Next Article in Journal
Optimal Effort on Self-Insurance-Cum-Protection: A New Analysis Using Yaari’s Dual Theory
Next Article in Special Issue
On the Search for a Measure to Compare Interval-Valued Fuzzy Sets
Previous Article in Journal
SEIR Mathematical Model of Convalescent Plasma Transfusion to Reduce COVID-19 Disease Transmission
Previous Article in Special Issue
Fuzzy Logical Algebra and Study of the Effectiveness of Medications for COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computing the Number of Failures for Fuzzy Weibull Hazard Function

by
Hennie Husniah
1,* and
Asep K. Supriatna
2
1
Department of Industrial Engineering, Langlangbuana University, Bandung 40261, Indonesia
2
Department of Mathematics, Padjadjaran University, Jatinangor 45363, Indonesia
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(22), 2858; https://doi.org/10.3390/math9222858
Submission received: 30 September 2021 / Revised: 30 October 2021 / Accepted: 5 November 2021 / Published: 10 November 2021
(This article belongs to the Special Issue Fuzzy Sets, Fuzzy Logic and Their Applications 2021)

Abstract

:
The number of failures plays an important factor in the study of maintenance strategy of a manufacturing system. In the real situation, this number is often affected by some uncertainties. Many of the uncertainties fall into the possibilistic uncertainty, which are different from the probabilistic uncertainty. This uncertainty is commonly modeled by applying the fuzzy theoretical framework. This paper aims to compute the number of failures for a system which has Weibull failure distribution with a fuzzy shape parameter. In this case two different approaches are used to calculate the number. In the first approach, the fuzziness membership of the shape parameter propagates to the number of failures so that they have exactly the same values of the membership. While in the second approach, the membership is computed through the α-cut or α-level of the shape parameter approach in the computation of the formula for the number of failures. Without loss of generality, we use the Triangular Fuzzy Number (TFN) for the Weibull shape parameter. We show that both methods have succeeded in computing the number of failures for the system under investigation. Both methods show that when we consider the function of the number of failures as a function of time then the uncertainty (the fuzziness) of the resulting number of failures becomes larger and larger as the time increases. By using the first method, the resulting number of failures has a TFN form. Meanwhile, the resulting number of failures from the second method does not necessarily have a TFN form, but a TFN-like form. Some comparisons between these two methods are presented using the Generalized Mean Value Defuzzification (GMVD) method. The results show that for certain weighting factor of the GMVD, the cores of these fuzzy numbers of failures are identical.

1. Introduction

Uncertainty is present in almost all decision problems, including in the field of reliability and maintenance. This is due to unknown future events and imprecision as well as human subjectivity in a decision process [1]. There are some important factors that significantly affect the decision-making in any field. In the field of reliability and maintenance, the number of failures plays important roles in the study of maintenance strategy of a manufacturing system. In the real situation, this number is often affected by some uncertainties. Many of the uncertainties fall into the possibilistic uncertainty, which is different from the probabilistic uncertainty. In many cases, at least one of the parameters or variables of the decision function has fuzzy value, instead of crisp value. This uncertainty is commonly modeled by applying fuzzy theoretical framework, e.g., the variable and parameter have fuzzy values and the calculation is done using extension principle approach [2].
As an important factor, the number of failures is essential to obtain, and subsequently is used as a base for further decision processes in reliability and maintenance analysis. As an example, this “number” is used in the calculation to design optimal maintenance strategies which are directed to minimize the number of failures while also minimizing the costs of operation [3,4,5]. For this reason, the knowledge on how to compute or predict the number of failures becomes vital. Considering the occurrence of uncertainty and imprecision—together with complexity of the system under investigation, failure data are often difficult to obtain. In this case, the theory of fuzzy sets has been widely used to provide a framework to deal with these uncertainty and imprecision [6]. Among the important questions needed to be addressed related to the number of failures of a system having a possibilistic uncertainty is, first, how to compute this number for a given possibility distribution with fuzzy parameters. Nowadays, some calculator for fuzzy numbers are readily available [2]. Second, it is also important to know how the degree of uncertainty of the parameters propagates to the resulting failure numbers. This is commonly known as the propagation of fuzziness, which is defined as “the way in which the amount of imprecision in the model’s inputs affects the changes in the model’s output” [7], (p. 163). Technically the propagation of uncertainty happens through mathematical operations involved in the model and in the computation. Knowing the method to calculate the number and its degree of uncertainty, will significantly improve the quality of the decision being sought (see also [8,9] for similar cases in other area). In general, fuzziness propagation in complex engineering systems may constitute a significant challenge [10].
The aims of the paper are two-fold, namely, to calculate the number of failures for a system which has Weibull failure distribution with a fuzzy shape parameter and to understand how the fuzziness of this shape parameter propagates to the resulting number of failures. These two objectives constitute the importance and contributions of the work presented in this paper. In addition, in this paper we look for the number of failures and two different approaches are used to calculate this number. In the first approach, the fuzziness membership of the shape parameter propagates to the number of failures so that they have exactly the same values of the membership. While in the second approach, the membership is computed through the α-cut or α-level of the shape parameter.

Literature Review

As it is explained earlier, the motivation of the paper is due the importance of finding the number of failures in the field of maintenance strategy. Some examples of such importance can be seen in [11,12,13,14,15,16] from various perspectives. It is often found that most of the problems in maintenance engineering are finding optimal strategies that minimize the cost of operation to manipulate the system as well as minimize the number of failures of the system (e.g., [17]). In many cases, the number of failures is represented in its distribution function. Several type of distribution functions are commonly used to model the failures of an industrial equipment, among others is the Weibull distribution function together with its hazard function [18,19]. This distribution function could appear either in two-parameters model or in three-parameters model [20,21]. Hence, the Weibull distribution function plays vital roles in areas of research related to maintenance strategy in which understanding a system, predicting the outcome of a system, and prescribing an optimal intervention to obtain the best performance of a system are being sought. In fact, the spectrum of the area applications of the Weibull distribution is quite broad from engineering, social sciences, to biological and health problems.
Apart from the abundant usage of the Weibull distribution in many areas of research, especially in maintenance strategy, most of the analysis only consider the crisp form of data, i.e., ignoring the presence of possibilistic uncertainty which might be often found in many real phenomena. For example, in maintenance engineering, most maintenance models in literature mainly consider certain or crisp condition, e.g., [22,23]. However, as mentioned earlier, these kinds of models do not seem to fit in the real condition. Readers may find a brief review of the importance of the possibilistic uncertainty in [24]. In reliability and maintenance problems, uncertainty may affect the models, the nonhomogeneous Poisson process (NHPP), and the Weibull generalized renewal process parameters [25,26], and the probability distribution parameters [8], and it is important to know how this uncertainty propagates through the models which likely affect the insight and prediction from the models [27].
There are several approaches to model the possibilistic uncertainty, one of them is by applying the fuzzy number theory, as suggested by Zadeh [28,29] and his followers. The popularity of fuzzy number theory in reliability and maintenance literatures is now getting higher [30] resulting in new modeling approach in many aspects, like maintenance risk-based inspection interval optimization with fuzzy failure interaction for two-component repairable system [31] and many others. The authors in [32] are among the few authors who consider fuzzy Weibull distribution in their works. They consider the Weibull distribution function as a fuzzy function and use it for analyzing the behavior of an industrial system stochastically by utilizing vague, imprecise, and uncertain data, which in turn result in the reliability indices (such as hazard function, maintainability, etc.) of time of time varying failure rate instead of the constant failure rate for the system.
In general, the author in [33] (pp. 152–157) shows several methods on how to implement a fuzzy function in addressing problems with possibilistic uncertainty. The methods can be classified into three different ways depending to which aspect of the crisp function the fuzzy concept was applied, namely (i) crisp function with fuzzy constraint, (ii) crisp function which propagates the fuzziness of independent variable to dependent variable, and (iii) function that is itself fuzzy. However, in this paper we will only look for the number of failures by using the first and the second approaches above. The fuzziness of the shape parameter is assumed to propagate to the number of failures with the same form of fuzzy number membership in the first approach, as found in [34,35]. While in the second approach, the concept of α-cut or α-level of the fuzziness of the shape parameter is used in the computation to calculate the number of failures, as found in [36]. An example of the methodology on how to compare fuzzy numbers, such as those resulting from different approaches of fuzzy function concepts above can be seen in [37].
In this paper we re-visit the model in [34,35] by giving some more detail analysis and results discussed in those papers. The authors in [34] discussed the Weibull hazard function by assuming a fuzzy shape parameter, which conceptually can be used to compute the number of failures without actually showing the resulting number of failures (either in crisp number form or fuzzy number form). They show how to compute the fuzzy number of failures of Weibull hazard function in [35] by assuming a fuzzy shape parameter in the Weibull hazard function via the second approach in [33], (p. 154), i.e., by considering the Weibull function as a crisp function which propagates the fuzziness of independent variable to dependent variable. In this paper we use different approaches by considering the fuzziness of the shape parameter in the computation of the number of failures directly, through the concept of α-cut or α-level [33] (p. 130) and [38], (pp. 7–16). Further we discuss the generalized mean value defuzzification (GMVD) and use it to compare the resulting fuzzy number of failures from different approaches of computation. The proposed defuzzification method (GMVD) is able to find a crisp number which is close to the core of the triangular fuzzy number (TFN).
We organize the presentation of the paper as follows. Section 2 presents briefly some basic methods that are utilized in the preceding sections, namely, the Weibull distribution function, fuzzy number and its membership function, α-cut of a fuzzy number, defuzzification process with Generalized Mean Value Defuzzification (GMVD), and the number of failures for Weibull hazard function with fuzzy parameter. Section 3 gives the main results together with numerical examples to show the visual illustration of the main results. This includes the comparisons from two different methods, i.e., the results from the method considering propagation of the fuzziness of independent variable to dependent variable and the results from the α-cut method. Section 4 presents the discussions of the results and it is finally followed by concluding remarks and further direction of research in Section 5.

2. Materials and Methods

The object being investigated in this paper is the Weibull distribution function as a mathematical model describing the deterioration of life cycle of an industrial system or an equipment. This deterioration or failure data are commonly modeled by the Weibull distribution function such as found in [39]. The reason of popularity of the Weibull function is its flexibility, so that it can be regarded as the generalization of exponential and Rayleigh distribution functions, which are also commonly used in reliability and maintenance studies [40]. The Weibull distribution is a continuous probability distribution function having the form:
f ( t ) = β θ t θ β 1 exp t θ β , x > 0 ,   θ > 0 ,   β > 0 , 0 , otherwise .
where θ is the scale parameter and β is the shape parameter. The first mathematician who described it in detail is Waloddi Weibull in 1951. The Weibull distribution has a flexibility to model various lifetime data by changing the value of the shape parameter, e.g., if β = 1 the Weibull distribution is reduced to an exponential distribution and if β = 2 the Weibull distribution is identical to Rayleigh distribution [40]. Throughout the paper we will assume the scale parameter θ = 1 for some reasons. For example, this choice is sufficient in our context of maintenance modeling if we assume that the average of first failure of the equipment/system under investigation happens within one unit of time—say one month or one year, because of its warranty and good quality control. The authors in [41] give an example with θ = 1 in their simulation.
Using this Weibull distribution function we can calculate some reliability indices, such as hazard function, number of failures, mean time between failures, preventive maintenance time, and replacement time. The standard methods on the calculation of these indices, both for standard and complex systems, can be found among others in [42,43]. Details theory and applications of the Weibull distribution function can be found in [44,45]. In the Section 2.1 we present some concepts of fuzzy theoretical framework which are used in the subsequent method and analysis, namely, fuzzy number and its membership function, α-cut of a fuzzy number, defuzzification process with Generalized Mean Value Defuzzification.

2.1. Fuzzy Number and Its Membership Function

As an introduction to the section that follows we define several concepts of fuzzy number theory that will be used later on. A fuzzy number can be regarded as an extension of a real number in the sense that it has a membership function other than binary to represent uncertainty. Binary membership gives a crips value for the membership, i.e., either a member or not a member. Fuzzy number gives a wide spectrum of membership from zero (definitely not a member of a set S) to one (definitely a member of a set S). Technically, a fuzzy number A ˜ refers to a connected set of possible values, where each possible value of A ˜ , say a, has its own membership value in the interval [0,1]. This value that measures the degree of possibility for a to be a member of A ˜ is called the membership function, usually written as μ : a A x [ 0 , 1 ] . This fuzzy number is commonly written with the symbol A ˜ = ( A , μ ( A ) ) or alternatively A ˜ = ( x , μ A ˜ ( x ) ) x X representing the underlying connected set A with the membership function μ ( A ) . In this regards, the fuzzy number is viewed as a pair of mathematical objects comprising of a set together with its grade or membership function. The fuzzy number is purportedly designed to represent the possibilistic uncertainty and to quantify the unclear and inaccuracies of the abundance of information.
The membership of a fuzzy number can be determined by several functional approaches, which can be classified into the linear and the non-linear functional forms. Among the most popular functional form of fuzzy number are the Triangular Fuzzy Number (TFN) which is often written as (a;b;c) and the Trapezoidal Fuzzy Number (TrFN) which is often written as (a;b;c;d). The functional forms or the membership functions of these fuzzy numbers are given in Equations (2) and (3), which are graphically shown in Figure 1 and Figure 2. Note that for the TrFN in Equation (3), the membership function within the intervals [a,b] and [c,d] are given by increasing and decreasing linear curves respectively. This concept is generalized by the LR-Flat Fuzzy Number which is then used as a new method for solving fuzzy transportation problems [2,46,47].
  • The membership function of a triangular fuzzy number (TFN):
μ A ˜ x = 0   ,   x a x a b a   ,   a x b x c b c   ,   b x c 0   ,   x c .
  • The membership function of a trapezoidal fuzzy number (TrFN):
μ A ˜ x = 0   ,   x a x a b a   ,   a x b 1   ,   b x c x d c d   ,   c x d 0   ,   x d .
In Equation (2), a, b, and c are real numbers satisfying a < b < c which constitute the TFN core and support components. In this case b is called the core of the fuzzy number and the sets [a,b) and (b,c] are called the support of the fuzzy number. Similarly, for TrFN, in Equation (3) the core of the fuzzy number is given by [b,c] and the support is given by the set [a,b) and (c,d]. Other forms of fuzzy numbers are piecewise quadratic fuzzy number [48], pentagonal fuzzy number [49], Bell shaped fuzzy number [50], parabolic trapezoidal fuzzy number [51], new bell shaped fuzzy number [52], and many others. A good reference on how some new methods and techniques are developed to advance fuzzy numbers concepts for modern analytics can be found in [46]. However, for simplicity, to emphasize the methodological aspect all examples in this paper assume the triangular fuzzy numbers (TFN). In the next section we briefly describe the α-cut of a triangular fuzzy numbers (α-cut, α-level cut, α-level set or sometimes simply is called a cut).

2.2. The α-Cut of a Fuzzy Number

Each fuzzy number is associated with its α-cut. This α-cut sometimes is also called the α-level set. It is technically defined as the set of objects in the associated fuzzy set which have the membership with the values which are at least α. This actually can be seen as a crisp set representation of a fuzzy number. Following this definition, it can be shown that the α-cut of the triangular fuzzy number (1) is given by:
A ˜ α = [ a 1 α , a 2 α ] = ( b - a ) α + a , ( b - c ) α + c
for all α [ 0 , 1 ] .

2.3. Generalized Mean Value Defuzzification

For some reasons, the information regarding the best representation of a crisp number for a fuzzy number is needed. In this case, defuzzification of the fuzzy number is done. It is a mathematical calculation which converts the fuzzy number into a single crisp value with respect to a fuzzy set. Some defuzzification formulas are available in literature, such as basic defuzzification distributions, center of area, center of gravity, fuzzy mean, last of maxima, weighted fuzzy mean, etc., [53,54,55]. In this paper we will use the generalized mean value defuzzification method (GMVD) which is defined as
N ( A ˜ ) = a + n b + c n + 2 ,
where A ˜ = ( a ; b ; c ) is a TFN and n can be regarded as the weight of the core of the fuzzy number. The larger the weight of the core, the closer the resulting crisp number from the GMVD to the core of the fuzzy number. The properties of this GMVD will be discussed later on and used in the comparation of the resulting number of failures.

3. Results

3.1. Number of Failures for Weibull Hazard Function with Fuzzy Parameter

As explained in the previous section we consider the one-parameter Weibull distribution function, since this choice is sufficient in our context of maintenance modeling if we assume that the average of first failure of the equipment/system under investigation happens within one unit of time—say one month or one year—because of its warranty and good quality control. By considering this assumption (θ = 1) and fuzzy parameter β ˜ the number of failure is computed using the first method, in which the calculation of the fuzzy number is done point-wise (will be defined later), and we only need the crisp function for the computation. From Equation (1) we have the following one-parameter Weibull cumulative distribution, g, and its hazard function, h:
g ( t ) = 1 e t β ,
and
h ( t ) = β t β 1 ,
so that the number of failures is given by
N ( t ) = t β .
The parameter β ˜ is the fuzzy number of the shape parameter of the Weibull function. We will treat the fuzziness of the shape parameter in two different approaches: (i) Crisp function which propagates the fuzziness of independent variable to dependent variable and, in which the computation is done point-wise; (ii) crisp function with fuzzy constraint through the level-set computation.
The First Method (Point-wise Method): Let β ˜ be a TFN which is identified by three crisp numbers a, b, and c, i.e., β ˜ = a ; b ; c satisfying Equation (2). We compute the number of failures point-wise, i.e., by substituting these crisp numbers one at a time to obtain the crisp output, say a’, b’, and c’. By assuming the same fuzzy measure propagates to the output, we will have µ(a’) = µ(a), µ(b’) = µ(b), and µ(c’) = µ(c), which give a TFN fuzzy output (a’; b’; c’) for the function g(t), h(t), and N(t) [34].
The Second Method (α-Cut Method): In the second approach, the fuzzy number β ˜ is identified as an α-cut satisfying Equation (4). As it is explained in [34], the fuzzy number of the shape parameter is approximated by a sequence of interval associated with the number α in [0,1]. This sequence consists of crisp numbers in the interval indicating the support of the fuzzy number for every α in [0, 1). If α is one then the supports converge to/become the core of the fuzzy number. The calculation to determine the number of failures is done at the end points of the interval. In this case, the stack of the end points of the intervals need not to be a TFN, which in many cases forms a TFN-like form (see numerical examples for the details).
To facilitate comparison between the results from the two methods, we use the GMVD defined in Equation (5). This GMVD has the properties as described in Theorem 1.
Theorem 1.
Let a TFN is given by (a;b;c), then the generalized mean value defuzzification (GMVD) defined by Equation (5) has the following properties:
  • For a symmetrical case, i.e., b a = c b = Δ then N ( A ˜ ) = b
  • For an asymmetrical case, i.e., b a = Δ a c b = Δ c then
    • N ( A ˜ ) > b if Δ a < Δ c
    • N ( A ˜ ) < b if Δ a > Δ c
  • If n then N ( A ˜ ) = b regardless the value of p and q .
Proof of Theorem 1:
  • Case 1: symmetrical TFN, i.e., b a = c b = Δ then
    N ( A ˜ ) = a + n b + c n + 2 = a + n b + ( a + 2 Δ ) n + 2 = 2 a + n b + 2 Δ n + 2 = 2 ( a + Δ ) + n b n + 2 = 2 b + n b n + 2 = ( 2 + n ) b n + 2 = b .
    Hence, N ( A ˜ ) = b .
  • Case 2: non-symmetrical TFN, i.e., b a = Δ a c b = Δ c then
    • if Δ a < Δ c then
      N ( A ˜ ) = a + n b + c n + 2 = a + n b + ( a + Δ a + Δ c ) n + 2 > a + n b + ( a + 2 Δ a ) n + 2 = 2 a + n b + 2 Δ a n + 2
      Hence, N ( A ˜ ) > b .
    • if Δ a > Δ c then
      Ν ( A ˜ ) = a + n b + c n + 2 = a + n b + ( a + Δ a + Δ c ) n + 2 < a + n b + ( a + 2 Δ a ) n + 2 = 2 a + n b + 2 Δ a n + 2
      Hence, N ( A ˜ ) < b .
  • If n then lim n N ( A ˜ ) = lim n a + n b + c n + 2 = b . □
As shown by the theorem, the GMVD above has a special characteristic, i.e., it is able to find a crisp number which is close to the core of the triangular fuzzy number (TFN). As examples, first consider the symmetrical TFN in Figure 2(left), i.e., β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) . It has GMVD = 1.55 for n = 1 and GMVD = 1.55 for n = 1000. Since it is symmetrical, the values of GMVD are the same as the core of the TFN for all n. However, for the non-symmetrical TFN, such as skewed left TFN β ˜ = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) in Figure 2(right), it has GMVD = 2.6833 for n = 1 and GMVD = 2.74980 for n = 1000. In this case, the larger is n the closer it is to the core of the TFN, i.e., 2.75. We will use this method of defuzzification for comparing the fuzzy output from two different methods in this paper.
Next we look at the fuzzy number of failures generated by the Weibull distribution via the α-cut method. Let us recall the α-cut of the triangular fuzzy number A ˜ = ( a ; b ; c ) is given by A ˜ α = [ a 1 α , a 2 α ] = ( b - a ) α + a , ( b - c ) α + c then the shape parameter, the Weibull cumulative distribution, the Weibul hazard function, and the number of failures are, respectively, given in the form of α-cut as follows. The shape parameter will have the form
β α = [ x 1 + x 3 α , x 2 x 3 α ] ,
for some x 1 ,   x 2 ,   x 3   R .By considering the α-cut in Equation (9) and substituting it into Equations (6) and (7) using the fuzzy arithmetic give rise to the cumulative distribution
g ( t ) α = [ 1 - exp ( - t y 1 + y 3 α ) , 1 - exp ( - t y 2 - y 3 α ) ] ,
for some y 1 ,   y 2 ,   y 3   R and the hazard function
h ( t ) α = [ ( z 1 + z 3 α ) t z 4 + z 6 α , ( z 2 z 3 α ) t z 5 z 6 α ] ,
for some z 1 ,   z 2 ,   z 3 ,   z 4 ,   z 5 ,   z 6   R . So that by integrating both sides of Equation (11) we end up with the number of failures, which is given by
N ( t ) α = [ t u 1 + u 3 α , t u 2 u 3 α ]
for some u 1 ,   u 2 ,   u 3   R .
The following theorem shows that as time goes, the GMVD of the number of failures increases and the support of the number of failures becomes wider. This means that the degree of uncertainty becomes larger.
Theorem 2.
For Δ t > 0 let N ( t ) α and N ( t + Δ t ) α be the fuzzy number of failures at time t and t + Δ t , respectively, then:
  • N ( t ) α = t p α ,   t s α and N ( t + Δ t ) α = t + Δ t p α ,   t + Δ t s α ,
  • G M V D N ( t + Δ t ) α G M V D N ( t ) α for all t R + ,
  • ( t + Δ t ) p α ( t + Δ t ) s α t p α t s α 0 for all t Z + .
Proof of Theorem 2:
  • It is clear.
  • It can be proved by using Theorem 1.
  • Note that for every α [ 0 , 1 ] , the interval in Equation (12) has the form t p α ,   t s α for some p α ,   s α     R . Without loss of generality, we will drop the index α, so that to prove the theorem we need ( t + Δ t ) p ( t + Δ t ) s t p t s 0 .
Consider the following binomial rule,
( x + Δ x ) n = k = 0 n n ! ( n k ) ! k ! Δ x n k x k . Then we have
( x + Δ x ) n = k = 0 n 1 ( n 1 ) ! ( ( n 1 ) k ) ! k ! Δ x ( n 1 ) k x k + n ! ( n k ) ! k ! Δ x n k x k = k = 0 n 1 ( n 1 ) ! ( ( n 1 ) k ) ! k ! Δ x ( n 1 ) k x k + x n .
Using this rule then for p , s Z + we have
( t + Δ t ) p = k = 0 p 1 ( p 1 ) ! ( ( p 1 ) k ) ! k ! Δ t ( p 1 ) k y k + t p
( t + Δ t ) s = k = 0 s 1 ( s 1 ) ! ( ( s 1 ) k ) ! k ! Δ t ( s 1 ) k y k + t s
A little algebraic manipulation gives
( t + Δ t ) p ( t + Δ t ) s t p t s = k = 0 p 1 ( p 1 ) ! ( ( p 1 ) k ) ! k ! Δ t ( p 1 ) k t k k = 0 s 1 ( s 1 ) ! ( ( s 1 ) k ) ! k ! Δ t ( s 1 ) k t k , = k = s p 1 ( p 1 ) ! ( ( p 1 ) k ) ! k ! Δ t ( p 1 ) k t k 0 ,
Which shows that ( t + Δ t ) p ( t + Δ t ) s t p t s 0 for all t Z + . Note that the theorem can be extended to any case of p , s R + . One can prove this using the Newton’s generalized Binomial theorem [56,57] in the form of infinite series rather than an infinite sum such as in the above case of t Z + . □

3.2. Numerical Examples

To obtain better insight regarding the results presented in the previous section we illustrate the concept above by using two different values for the shape parameters, the relatively small value β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and the relatively large value β ˜ = ( p = 2.50 ;   q = 2.75 ;   s   : = 2.80 ) . Here p, q, and s are the TFN components which constitute the TFN defined just the same as a, b, and c in Equation (2). The graphs of these TFNs are shown in Figure 2. For the first method, the number of failures for the shape parameters in Figure 2 at t = 10 is presented in Figure 3 while Figure 4 (top figures) shows the number of failures for t in [0,100] with 10 steps size. Figure 4 (bottom figures) shows the nonlinearity of the failure numbers as a function of t. Similarly, for the second method, the number of failures for the shape parameters in Figure 2 at t = 10 is presented in Figure 5 while Figure 6 shows the number of failures for t in [0,100] with all steps of time. For the finer step size, i.e., 100 steps size, the graph of the number of failures from the second method is presented in Figure 7. Clearly the number of failures in Figure 3 are in triangular forms since the first method assumes that the fuzziness of the shape parameter propagates to the number of failures with the same form of fuzzy number membership, while the number of failures in Figure 5 does not have a triangular form since the fuzziness uncertainty is considered and affects the functional calculation of the number of failures through the α-cut arithmetic. Figure 8 gives the comparisons between these two relatively different shapes. Further, if we plot the numbers of failures over time (see bottom figures in Figure 4), then the curves are non-linear and seem to be “exponentially” increase as expected in the theory. The bottom graphs in Figure 4 actually show the numbers of failures over time for the end points and core of the shape parameter TFNs. To be exact these figures show the graphs of Weibull’s numbers of failures bands, which analytically is given by Equation (8) and comparable to Equations (10) and (14) for the α-cat, hence it has a power curve. This is consistent with the curve for Weibull’s number of failures with crisp parameters [58]. This is also true for the second method (the α-cut approach), but we do not show the graphs here.
The figures show that for both values of fuzzy shape parameters β ˜ , the relatively small value β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and the relatively large value β ˜ = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) , the length of the fuzziness of the resulting number of failures get bigger as the time t increases. This means the increase of the possibilistic uncertainty of the number of failures. This phenomenon also appears in the α-cut method as is shown in the next section.

3.3. Results from the α-Cut Method

The following results are plotted from the calculation of the number of failures using the α-cut method. Recall the α-cut of the triangular fuzzy number A ˜ = ( a ; b ; c ) is given by A ˜ α = [ a 1 α , a 2 α ] = ( b - a ) α + a , ( b - c ) α + c hence for the fuzzy shape parameter β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) we obtain its α-cut is
β α = [ 1.25 + 0.30 α , 1.85 0.30 α ] ,
as the fuzzy number of the shape parameter. By considering the α-cut in Equation (7) and substituting it into Equations (5) and (6) using the fuzzy arithmetic give rise to the cumulative distribution
g ( t ) α = [ 1 - exp ( - t 1.25 + 0.30 α ) , 1 - exp ( - t 1.85 - 0.30 α ) ] ,
and the hazard function
h ( t ) α = [ ( 1.25 + 0.30 α ) t 0.25 + 0.30 α , ( 1.85 0.30 α ) t 0.85 0.30 α ] ,
So that by integrating both sides of Equation (9) we end up with the number of failures, which is given by
N ( t ) α = [ t 1.25 + 0.30 α , t 1.85 0.30 α ] .
When we use the α-cut method, we will have a triangular-like fuzzy number which is comparable (not necessarily the same) to the triangular fuzzy number (p;q;r) defined by:
p = min N ( t ) α = 0 = t 5 / 4 ,
q = N ( t ) α = 1 = t 31 / 20 ,
r = min N ( t ) α = 0 = t 37 / 20 ,
We enumerate the fuzzy number of failures in Table 1 based on the calculation of these formulas for t = 0 to t = 10.
Table 1 also gives the counterpart of the fuzzy number of failure calculated by the first method. Note that in Table 1, TFN is a triangle fuzzy number, while FN is not necessarily a triangle fuzzy number. However, they both have the same core and the same support but the shapes are different (see Figure 8).
Further, to compare the resulting fuzzy number of failures among the methods, we defuzzified them using the generalized mean value defuzzification (GMVD) which is defined by (4) with n = 4. The comparison shows that the defuzzified numbers both from the first method and the second method are exactly the same to the results from the crisp method. Table 2 shows that if n is getting larger, then, the defuzzified number gets closer to the core of the fuzzy number, e.g., for t = 10, with n = 1,000,000 the defuzzified number is 35.4813565346595 which approaches the core of its fuzzy number, i.e., q = 35.48133892. This agrees with Theorem 1. We plot the resulting number of failures for t = 10 in Figure 5 and for t = 0 to t = 10 in Figure 6. The same procedure is done for the relatively large value of the shape parameter β ˜ = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) but the details are not presented here. The plots are presented in the righthand side of Figure 5, Figure 6 and Figure 7.
The time-series plots of the Cumulative Distribution Function, the Hazard Function, and the Number of Failures are presented in Figure 9. The shape parameter on the upper-left of Figure 9 is β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and on the upper-right of Figure 9 is β = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) . The figure shows the plots for a short period of time, up to t = 1.5.

4. Discussions

The analytical results in Theorems 1 to 3 are illustrated by numerical examples to gain visual understanding regarding the analytical finding above by using β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and β ˜ = ( p = 2.50 ;   q = 2.75 ;   s   : = 2.80 ) reflecting a relatively small and a relatively large shape parameter, respectively. Here p, q, and s are the TFN components which constitute the TFN defined just the same as a, b, and c in Equation (1). See Figure 2 for the graphs of these TFNs and Figure 3 for the resulting number of failures at t = 10 from the first method and Figure 5 (for the second method). While Figure 4 (top figures) shows the number of failures for t in [0,100] for the first method and Figure 6 for the second method with 10 steps size, (for the finer step size, i.e., 100 steps size see Figure 7). Clearly the number of failures in Figure 3 are in triangular forms due to the assumption in the first method in which the fuzziness of the shape parameter propagates with the same form of fuzzy number membership to the number of failures, while the number of failures in Figure 5 does not have a triangular form since the fuzziness uncertainty is considered and affecting the functional calculation of the number of failures through the α-cut arithmetic. Figure 8 gives the comparisons between these two relatively different shapes. The time-series plots of the cumulative distribution function, the hazard function, and the number of failures are presented in Figure 9. All curves are familiar in shape as it conform to their crisp parameter of Weibull distribution, but here they form twisted-cumulative band, -hazard band, and -number of failures band instead of single curve, respectively.
Furthermore, if we plot the numbers of failures over time (see bottom figures in Figure 4), then the curves are non-linear and seem to “exponentially” increase as expected in the theory. The bottom graphs in Figure 4 actually show the numbers of failures over time for the end points and core of the shape parameter TFNs. To be exact these figures show the graphs of Weibull’s numbers of failures bands, which analytically is given by Equations (12) and (16), hence it has a power curve shape which conforms to the known curve for Weibull’s number of failures with crisp parameters [58]. This is also true for the second method (the α-cut approach), but we do not show the graphs here.
When considering a Weibull distribution with fuzzy shape parameter to calculate the fuzzy number of failures, usually in such imprecise situations, extension principle approach is often used as one choice of calculation though it could lead to a complex form. Here we have proposed a simple method (the first method) to calculate the number of failures, by assuming that the fuzziness of the shape parameter propagates to the number of failures with the same form of fuzzy number membership, and also proposed an alternative method (the second method) which is the calculation done using the α-cut method. This method could be extended to the Weibull distribution with more parameters to enlarge the applicability to other area [59].

5. Conclusions

In this paper we have discussed the Weibull hazard function by assuming a fuzzy shape parameter to calculate the fuzzy number of failures. Here we have proposed a simple method (the first method) to calculate the number of failures, by assuming that the fuzziness of the shape parameter propagates to the number of failures with the same form of fuzzy number membership, and also proposed an alternative method (the second method) which is the calculation done using the α-cut method.
We have shown that both methods have succeeded in computing the number of failures for the system under investigation. Both methods show that when we consider the function of the number of failures as a function of time then the uncertainty (the fuzziness) of the resulting number of failures becomes larger and larger as the time increases. This indicates the propagation of uncertainty in the shape parameter into the resulting number of failures, in which for large values of t, a small value of uncertainty in the shape parameter will produce a large support to the fuzzy number of failures. In practical implication, one should be aware of these properties when using the resulting number of failures as a base for the further process of decision-making.
In this paper we have used a TFN for the shape parameter and by using the first method, the resulting number of failures has a TFN form. Meanwhile, the resulting number of failures from the second method does not necessarily have a TFN form, but a TFN-like form. Some comparisons between these two methods are presented using the Generalized Mean Value Defuzzification (GMVD) method. The results show that for certain weighting factor of the GMVD, the cores of these fuzzy numbers of failures are identical. We did the comparation between the two methods after we use the GMVD which produces crisp number of failures. This can be regarded as a shortcoming of the study since once we defuzzify the resulting number of failures we lose the information of the uncertainty. Further study can be done by considering the comparation with a method that preserves the uncertainty.
The TFN form and value of the shape parameter used in the Weibull distribution function was taken for granted. For the practical applications this would be not easy. The true form of the fuzzy number should be correctly decided from the available real data and the value should be estimated from the same data. These issues are among the limitations of the methods presented here and could also lead to future direction of research. Other concern is that here we only consider one parameter which has fuzzy value. In realty all of the Weibull parameters could also have imprecise measure or uncertainty. This also will lead to important future venue of research (currently four-parameter Weibull distribution has already around in crisp value application ref). Here we only consider one-parameter Weibull distribution by assuming the scale parameter is assumed to be one. This is sufficient in our context of maintenance modeling if we assume that the average of first failure of the equipment/system under investigation happens within one unit of time. However, in general case this may not be true, so we need to extend the analysis into Weibull distribution having arbitrary values of the scale parameters. Further studies can also be done for different approaches with different forms of fuzzy numbers, different uses of defuzzification methods, and explore the applications of the theory in different related field, such as the number of failures in biological processes (e.g., failure in protecting healthy status (susceptibility) for people who are infected by COVID-19 disease), which currently we are working on.

Author Contributions

H.H.: conceptualization, investigation, project administration, resources, data curation, validation, writing—original draft preparation; A.K.S.: methodology and software, formal analysis, visualization, supervision, writing—review and editing, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by the Ministry of Research, Technology, and Higher Education of the Republic of Indonesia, through the scheme of “Penelitian Hibah Riset dan Teknologi/Badan Riset dan Inovasi Nasional”, with grant number 1207/UN6.3.1/PT.00/202.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the anonymous referees who have read the earlier version of the manuscript and raised many constructive comments that helped to improve the content and the presentation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esogbue, A.O.; Hearnes, W.E. On replacement models via a fuzzy set theoretic framework. IEEE Trans. Syst. Man Cybern. Part C Appl. Rev. 1998, 28, 549–560. [Google Scholar] [CrossRef]
  2. Parchami, A. Calculator for fuzzy numbers. Complex. Intell. Syst. 2019, 5, 331–342. [Google Scholar] [CrossRef] [Green Version]
  3. Husniah, H.; Supriatna, A.K.; Iskandar, B.P. Lease contract with servicing strategy model for used product considering crisp and fuzzy usage rate. Int. J. Artif. Intell. 2018, 18, 177–192. [Google Scholar]
  4. Husniah, H.; Pasaribu, U.S.; Iskandar, B.P. Multi-period lease contract for remanufactured products. Alex. Eng. J. 2020, 60, 2279–2289. [Google Scholar] [CrossRef]
  5. Husniah, H.; Pasaribu, U.S.; Wangsaputra, R.; Iskandar, B.P. Condition-based maintenance policy for a leased reman product. Heliyon 2021, 7, e06494. [Google Scholar] [CrossRef] [PubMed]
  6. Verma, A.K. The use of fuzzy numbers in reliability calculations. IETE Tech. Rev. 2001, 18, 27–31. [Google Scholar] [CrossRef]
  7. Glöckner, I. Fuzzy Quantifiers: A Computational Theory; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  8. Venkatesh, A.; Elango, S. A mathematical model for the effect of TRH using fuzzy reliability analysis. Glob. J. Pure Appl. Math. 2017, 13, 5673–5686. [Google Scholar]
  9. Venkatesh, A.; Subramani, G.; Elango, S. A mathematical model for the effect of gastrin in humans using fuzzy Weibull distribution. Int. J. Pure Appl. Math. 2017, 117, 155–164. [Google Scholar]
  10. Adhikari, S.; Khodaparast, H.H. A spectral approach for fuzzy uncertainty propagation in finite element analysis. Fuzzy Sets Syst. 2014, 243, 1–24. [Google Scholar] [CrossRef]
  11. Khalaf, A.; Djouani, K.; Hamam, Y.; Alayli, Y. Failure-cost mathematical model for medical equipment maintenance. In Proceedings of the Fourth International Conference on Advances in Biomedical Engineering (ICABME), Beirut, Lebanon, 19–21 October 2017; pp. 1–4. [Google Scholar] [CrossRef]
  12. Li, G.; Li, Y.; Zhang, X.; Hou, C.; He, J.; Xu, B.; Chen, J. Development of a preventive maintenance strategy for an automatic production line based on group maintenance method. Appl. Sci. 2018, 8, 1781. [Google Scholar] [CrossRef] [Green Version]
  13. Lee, J.; Kim, B.; Ahn, S. Maintenance optimization for repairable deteriorating systems under imperfect preventive maintenance. Mathematics 2019, 7, 716. [Google Scholar] [CrossRef] [Green Version]
  14. Husniah, H.; Herdiani, L.; Kusmaya, K.; Supriatna, A.K. Fuzzy usage pattern in customizing public transport fleet and its maintenance options. J. Phys. Conf. Ser. 2018, 1013, 012186. [Google Scholar] [CrossRef] [Green Version]
  15. Mohammed, A.; Ghaithan, A.; Al-Saleh, M.; Al-Ofi, K. Reliability-based preventive maintenance strategy of truck unloading systems. Appl. Sci. 2020, 10, 6957. [Google Scholar] [CrossRef]
  16. Nordal, H.; El-Thalji, I. Assessing the technical specifications of predictive maintenance: A case study of centrifugal compressor. Appl. Sci. 2021, 11, 1527. [Google Scholar] [CrossRef]
  17. Husniah, H.; Supriatna, A.K.; Iskandar, B.P. Improving service contract management with availability and negotiation. In Interdisciplinary Behavior and Social Sciences, Proceedings of the 3rd International Congress on Interdisciplinary Behavior and Social Science 2014 (ICIBSoS 2014), Bali, Indonesia, 1–2 November 2014; Gaol, L., Ed.; CRC Press-Taylor and Francis Group: London, UK, 2015; pp. 409–413. [Google Scholar]
  18. Wisniewski, R. Using Weibull Analysis to Guide Preventative Maintenance Strategy. Def. Syst. Inf. Anal. Cent. 2014, 2. Available online: https://dsiac.org/articles/using-weibull-analysis-to-guide-preventative-maintenance-strategy/ (accessed on 2 September 2021).
  19. Pascale, E.; Freneaux, T.; Sista, R.; Sannino, P.; Marmo, P.; Bouillaut, L. Application of the Weibull distribution for the optimization of maintenance policies of an electronic railway signaling system. In Proceedings of the ESREL 2017 European Safety and Reliability Conference, Portoroz, Slovenia, 18–22 June 2017; Available online: https://hal.archives-ouvertes.fr/hal-01521640/file/doc00026876.pdf (accessed on 2 September 2021).
  20. Hong, J. Modified Weibull Distributions in Reliability Engineering. Ph.D. Thesis, National University of Singapore, Singapore, 2010. Available online: http://scholarbank.nus.edu.sg/handle/10635/17099 (accessed on 2 September 2021).
  21. Baloui Jamkhaneh, E. Analyzing system reliability using fuzzy weibull lifetime distribution. Int. J. Appl. Oper. Res. 2014, 4, 81–90. [Google Scholar]
  22. Husniah, H.; Supriatna, A.K. Optimal number of fishing fleet for a sustainable fishery industry with a generalized logistic production function. In Proceedings of the 2015 International Conference on Industrial Engineering and Systems Management (IESM), Seville, Spain, 21–23 October 2015. [Google Scholar] [CrossRef]
  23. Husniah, H.; Pasaribu, U.S.; Supriatna, A.K.; Iskandar, B.P. Optimal number of fleet maintenance contract with policy limit cost. In Proceedings of the 4th International Conference on Control, Decision and Information Technologies (CoDIT), Bandung, Indonesia, 5–7 April 2017. [Google Scholar] [CrossRef]
  24. Zio, E.; Pedroni, N. Possibilistic methods for uncertainty treatment applied to maintenance policy assessment. In Industrial Safety Cahiers; Foundation for an Industrial Safety Culture: Toulouse, France, 2014; no. 2014-07; Available online: https://www.foncsi.org/en/publications/collections/industrial-safety-cahiers/possibilistic-methods-uncertainty-maintenance-policy/CSI-uncertainty-maintenance.pdf (accessed on 2 September 2021).
  25. Basile, O.; Dehombreux, P.; Riane, P. Identification of reliability models for non repairable and repairable systems with small samples. In Faculte Polytechnique de Mons Research Report; Arles: Mons, Belgium, 2004; Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.595.8268&rep=rep1&type=pdf (accessed on 2 September 2021).
  26. Basile, O.; Dehombreux, P.; Riane, P. Evaluation of the uncertainty affecting reliability models. J. Qual. Maint. Eng. 2007, 13, 137–151. [Google Scholar] [CrossRef]
  27. Baraldi, P.; Compare, M.; Zio, E. Uncertainty treatment in expert information systems for maintenance policy assessment. Appl. Soft Comput. 2014, 22, 297–310. [Google Scholar] [CrossRef] [Green Version]
  28. Zadeh, L. The concept of a linguistic variable and its application to approximate reasoning–I. Inf. Sci. 1975, 8, 199–249. [Google Scholar] [CrossRef]
  29. Zadeh, L. Fuzzy sets as a basis for a theory of possibility. Fuzzy Sets Syst. 1999, 100, 9–34. [Google Scholar] [CrossRef]
  30. Rotshtein, A.; Katielnikov, D.; Pustylnik, L. Reliability modeling and optimization using fuzzy logic and chaos theory. Int. J. Qual. Stat. Reliab. 2012, 2012, 847416. [Google Scholar] [CrossRef] [Green Version]
  31. Rezaei, E.; Imani, D.M. A new modeling of maintenance risk based inspection interval optimization with fuzzy failure interaction for two-component repairable system. Indian J. Nat. Sci. 2015, 6, 9003–9017. [Google Scholar]
  32. Garg, H.; Sharma, S.P.; Rani, M. Weibull fuzzy probability distribution for analysing the behaviour of pulping unit in a paper industry. Int. J. Ind. Syst. Eng. 2013, 14, 395–413. [Google Scholar] [CrossRef]
  33. Lee, K.H. First Course on Fuzzy Theory and Applications; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  34. Husniah, H.; Widjajani, W.; Rohmana, R.; Supriatna, A.K. Number of failures for weibull hazard function with a fuzzy shape parameter. J. Phys. Conf. Ser. 2019, 1280, 022034. [Google Scholar] [CrossRef]
  35. Husniah, H.; Supriatna, A.K. Application of fuzzy theory in determining the number of failures for Weibull hazard function. In Proceedings of the ICOCSIM, Lombok, Indonesia, 21 March 2019. [Google Scholar]
  36. Sharma, R.K.; Sharma, P. System failure behavior and maintenance decision making using, RCA, FMEA and FM. J. Qual. Maint. Eng. 2010, 16, 64–88. [Google Scholar] [CrossRef]
  37. Rahmani, A.; Lotfi, F.H.; Rostamy-Malkhalifeh, M.; Allahviranloo, T. A new method for defuzzification and ranking of fuzzy numbers based on the statistical beta distribution. Adv. Fuzzy Syst. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [Green Version]
  38. Buckley, J.J.; Eslami, E.; Feuring, T. Fuzzy Mathematics in Economics and Engineering; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  39. Lai, C.D.; Murthy, D.N.; Xie, M. Weibull distributions and their applications. In Springer Handbook of Engineering Statistics; Pham, H., Ed.; Springer: New York, NY, USA, 2006; pp. 63–78. [Google Scholar]
  40. Kumar, P.; Singh, S.B. Fuzzy system reliability using intuitionistic fuzzy Weibull lifetime distribution. Int. J. Reliab. Appl. 2015, 16, 15–26. [Google Scholar]
  41. Iskandar, B.P.; Jack, N.; Murthy, D.N.P. Two new servicing strategies for products sold with warranty. Asia Pac. J. Oper. Res. 2012, 29, 1240022. [Google Scholar] [CrossRef]
  42. Nakagawa, T. Maintenance Theory of Reliability; Springer-Verlag: London, UK, 2005. [Google Scholar]
  43. Kobbacy, K.A.H.; Murthy, D.N.P. Complex System Maintenance Handbook; Springer Nature: Cham, Switzerland, 2008. [Google Scholar]
  44. Rinne, H. The Weibull Distribution A Handbook, 1st ed.; Chapman and Hall/CRC: New York, NY, USA, 2008. [Google Scholar]
  45. McCool, J.I. Using the Weibull Distribution: Reliability, Modeling, and Inference; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef]
  46. Ebrahimnejad, A.; Verdegay, J.L. Fuzzy Sets-Based Methods and Techniques for Modern Analytics; Springer International Publishing: Cham, Switzerland, 2018. [Google Scholar] [CrossRef]
  47. Ebrahimnejad, A. New method for solving Fuzzy transportation problems with LR flat fuzzy numbers. Inf. Sci. 2016, 357, 108–124. [Google Scholar] [CrossRef]
  48. Jain, S. Close interval approximation of piecewise quadratic fuzzy numbers for fuzzy fractional program. Iran. J. Oper. Res. 2010, 2, 77–88. [Google Scholar]
  49. Panda, A.; Pal, M. A study on pentagonal fuzzy number and its corresponding matrices. Pac. Sci. Rev. B. Humanit. Soc. Sci. 2015, 1, 131–139. [Google Scholar] [CrossRef] [Green Version]
  50. Maturo, F.; Fortuna, F. Bell-shaped fuzzy numbers associated with the normal curve. In Topics on Methodological and Applied Statistical Inference; Di Battista, T., Moreno, E., Racugno, W., Eds.; Springer Nature: Cham, Switzerland, 2016; pp. 131–144. [Google Scholar] [CrossRef]
  51. Thangavelu, K.; Uthra, G.; Shunmugapriya, S. A new approach on the membership functions of fuzzy numbers. Int. J. Pure Appl. Math. 2017, 114, 145–152. [Google Scholar]
  52. Bagheri, N.; Hashemin, S.S. A new bell shape fuzzy number. Int. J. Math. Trends Technol. 2018, 54, 377–382. [Google Scholar] [CrossRef]
  53. van Leekwijck, W.; Kerre, E.E. Defuzzification: Criteria and classification. Fuzzy Sets Syst. 1999, 108, 159–178. [Google Scholar] [CrossRef]
  54. Madau, D.P.; Feldkamp, L.A. Influence value defuzzification method. In Proceedings of the IEEE 5th International Fuzzy Systems, New Orleans, LA, USA, 11 September 1996; pp. 1819–1824. [Google Scholar] [CrossRef]
  55. Jager, B.; Verbruggen, H.B.; Bruijn, P.M. The role of defuzzification methods in the application of fuzzy control. IFAC Proc. Umes 1992, 25, 75–80. [Google Scholar] [CrossRef]
  56. Available online: https://en.wikibooks.org/wiki/Advanced_Calculus/Newton%27s_general_binomial_theorem (accessed on 2 September 2021).
  57. Koh, Y.; Ree, S. The Origin of Newton’s Generalized Binomial Theorem. J. Hist. Math. 2014, 27, 127–138. [Google Scholar] [CrossRef] [Green Version]
  58. Scheideggera, A.; Scholtena, L.; Maurera, M.; Reichert, P. Extension of pipe failure models to consider the absence of data from replaced pipes. Water Res. 2013, 47, 3696–3705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Izadparast, H.A.; Niedzwecki, J.M. Four-parameter Weibull probability distribution model for weakly non-linear random variables. Probabilistic Eng. Mech. 2013, 32, 31–38. [Google Scholar] [CrossRef]
Figure 1. Graphical representation of a triangular fuzzy number (a;b;c)—left figure, and a trapezoidal fuzzy number (a;b;c;d)—right figure.
Figure 1. Graphical representation of a triangular fuzzy number (a;b;c)—left figure, and a trapezoidal fuzzy number (a;b;c;d)—right figure.
Mathematics 09 02858 g001
Figure 2. On the left figure is shown the relatively small shape parameter β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and on the right figure is shown the relatively large shape parameter β ˜ = ( p = 2.50 ;   q = 2.75 ;   s   : = 2.80 ) . The vertical axis is the fuzzy membership µ of the shape parameter’s TFN. The first shape parameter is a symmetrical TFN and the second shape parameter is a nonsymmetrical TFN. These TFNs are used to calculate their respective number of failures in the subsequent figures.
Figure 2. On the left figure is shown the relatively small shape parameter β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and on the right figure is shown the relatively large shape parameter β ˜ = ( p = 2.50 ;   q = 2.75 ;   s   : = 2.80 ) . The vertical axis is the fuzzy membership µ of the shape parameter’s TFN. The first shape parameter is a symmetrical TFN and the second shape parameter is a nonsymmetrical TFN. These TFNs are used to calculate their respective number of failures in the subsequent figures.
Mathematics 09 02858 g002
Figure 3. The left figure is the number of failures for the shape parameter β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) at t = 10—see left figure in Figure 1. The right figure is the number of failures for the shape parameter β ˜ = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) at t = 10—see right figure in Figure 2. Note that the vertical axis indicates the fuzzy membership µ.
Figure 3. The left figure is the number of failures for the shape parameter β ˜ = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) at t = 10—see left figure in Figure 1. The right figure is the number of failures for the shape parameter β ˜ = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) at t = 10—see right figure in Figure 2. Note that the vertical axis indicates the fuzzy membership µ.
Mathematics 09 02858 g003
Figure 4. The description is as in Figure 3 above but with t = 0 to t = 100 and step size of t is 10. The left axis is time, the right axis is the number of failures, and the vertical axis is the fuzzy membership degree of the number of failures (above). The figures in the bottom show the core (black), the lower bound (blue), and the upper bound (red) for the resulting number of failures with small shape parameter (left) and large shape parameter (right).
Figure 4. The description is as in Figure 3 above but with t = 0 to t = 100 and step size of t is 10. The left axis is time, the right axis is the number of failures, and the vertical axis is the fuzzy membership degree of the number of failures (above). The figures in the bottom show the core (black), the lower bound (blue), and the upper bound (red) for the resulting number of failures with small shape parameter (left) and large shape parameter (right).
Mathematics 09 02858 g004
Figure 5. The left figure is the number of failures for the shape parameter β = ( p   = 1.25 ;   q   = 1.55 ;   s   = 1.85 ) at t = 10. The right figure is the number of failures for the shape parameter β = ( p   = 2.50 ;   q   = 2.75 ;   s   = 2.80 ) at t = 10. Both figures are generated by the second method with 20 levels of α, i.e., α0 = 0 as the base to α21 = 1 as the peak.
Figure 5. The left figure is the number of failures for the shape parameter β = ( p   = 1.25 ;   q   = 1.55 ;   s   = 1.85 ) at t = 10. The right figure is the number of failures for the shape parameter β = ( p   = 2.50 ;   q   = 2.75 ;   s   = 2.80 ) at t = 10. Both figures are generated by the second method with 20 levels of α, i.e., α0 = 0 as the base to α21 = 1 as the peak.
Mathematics 09 02858 g005
Figure 6. The description is as in Figure 5 above but with complete steps form t = 0 to t = 10. The left axis is time, the right axis is the number of failures, and the vertical axis is the fuzzy membership degree of the number of failures.
Figure 6. The description is as in Figure 5 above but with complete steps form t = 0 to t = 10. The left axis is time, the right axis is the number of failures, and the vertical axis is the fuzzy membership degree of the number of failures.
Mathematics 09 02858 g006
Figure 7. The plots of the number of failures for the shape parameter β = ( p   = 1.25 ;   q   = 1.55 ;   s   = 1.85 ) and β = ( p   = 0.9 ;   q   = 1.0 ;   s   = 1.5 ) from the second method against time from t = 0 to t = 100 as in Figure 6 but with a finer step size of t (other parameters are the same as in Figure 5 and Figure 6).
Figure 7. The plots of the number of failures for the shape parameter β = ( p   = 1.25 ;   q   = 1.55 ;   s   = 1.85 ) and β = ( p   = 0.9 ;   q   = 1.0 ;   s   = 1.5 ) from the second method against time from t = 0 to t = 100 as in Figure 6 but with a finer step size of t (other parameters are the same as in Figure 5 and Figure 6).
Mathematics 09 02858 g007
Figure 8. The top and bottom figures are plots of the number of failures for β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and β = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) , respectively, with the left hand side is for t = 10 and the right hand side is for t = 100.
Figure 8. The top and bottom figures are plots of the number of failures for β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and β = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) , respectively, with the left hand side is for t = 10 and the right hand side is for t = 100.
Mathematics 09 02858 g008
Figure 9. The time-series plots of the Cumulative Distribution Function, the Hazard Function, and the Number of Failures. The shape parameter on left figure is β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and on the right figure is β = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) .
Figure 9. The time-series plots of the Cumulative Distribution Function, the Hazard Function, and the Number of Failures. The shape parameter on left figure is β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) and on the right figure is β = ( p = 2.50 ;   q = 2.75 ;   s = 2.80 ) .
Mathematics 09 02858 g009aMathematics 09 02858 g009b
Table 1. Number of failures comparisons for β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) . Note that for the α -cut method we use α = 0 to obtain the support (a,c) and α = 1 to find the core b of the resulting fuzzy number so that we have an analogous TFN (a;b;c).
Table 1. Number of failures comparisons for β = ( p = 1.25 ;   q = 1.55 ;   s = 1.85 ) . Note that for the α -cut method we use α = 0 to obtain the support (a,c) and α = 1 to find the core b of the resulting fuzzy number so that we have an analogous TFN (a;b;c).
Time tCrisp
Method
Fuzzy Propagation
Method
Fuzzy α-Cut Method
TFN
(p;q;s)
Defuzzification
(p + 4q + s)/6
TFN-like
(p;q;s)
Defuzzification
(p + 4q + s)/6
0000p = 0
q = 0
s = 0
0
1111p = 1
q = 1
s = 1
1
22.949350275p = 2.378414230
q = 2.928171392
s = 3.605001850
2.949350275p = 2.378414230
q = 2.928171392
s = 3.605001850
2.949350275
35.589852442p = 3.948222039
q = 5.489565165
s = 7.632631956
5.589852442p = 3.948222039
q = 5.489565165
s = 7.632631956
5.589852442
48.824940564p = 5.656854248
q = 8.574187700
s = 12.99603834
8.824940564p = 5.656854248
q = 8.574187700
s = 12.99603834
8.824940564
512.59725950p = 7.476743905
q = 12.11723434
s = 19.63787576
12.59725950p = 7.476743905
q = 12.11723434
s = 19.63787576
12.59725950
616.86728508p = 9.390507480
q = 16.07438767
s = 27.51565232
16.86728508p = 9.390507480
q = 16.07438767
s = 27.51565232
16.86728508
721.60548840p = 11.38603593
q = 20.41277093
s = 36.59581083
21.60548840p = 11.38603593
q = 20.41277093
s = 36.59581083
21.60548840
826.78864158p = 13.45434265
q = 25.10669114
s = 46.85074227
26.78864158p = 13.45434265
q = 25.10669114
s = 46.85074227
26.78864158
932.39780510p = 15.58845727
q = 30.13532570
s = 58.25707056
32.39780510p = 15.58845727
q = 30.13532570
s = 58.25707056
32.39780510
1038.41712138p = 17.78279410
q = 35.48133892
s = 70.79457844
38.41712138p = 17.78279410
q = 35.48133892
s = 70.79457844
38.41712138
Table 2. The illustration of Theorem 1 of the GMVD for non-symmetrical fuzzy number in Figure 5 (left). The fuzzy number is (17.782794100; 35.481338920; 70.794578440).
Table 2. The illustration of Theorem 1 of the GMVD for non-symmetrical fuzzy number in Figure 5 (left). The fuzzy number is (17.782794100; 35.481338920; 70.794578440).
nGMVDnGMVD
044,28868627637,68317576
141,35290382737,43852722
239,88501260837,24280839
339,00427786937,08267480
438,417121371036,94923015
537,9977238810,000,00035,48135653
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Husniah, H.; Supriatna, A.K. Computing the Number of Failures for Fuzzy Weibull Hazard Function. Mathematics 2021, 9, 2858. https://doi.org/10.3390/math9222858

AMA Style

Husniah H, Supriatna AK. Computing the Number of Failures for Fuzzy Weibull Hazard Function. Mathematics. 2021; 9(22):2858. https://doi.org/10.3390/math9222858

Chicago/Turabian Style

Husniah, Hennie, and Asep K. Supriatna. 2021. "Computing the Number of Failures for Fuzzy Weibull Hazard Function" Mathematics 9, no. 22: 2858. https://doi.org/10.3390/math9222858

APA Style

Husniah, H., & Supriatna, A. K. (2021). Computing the Number of Failures for Fuzzy Weibull Hazard Function. Mathematics, 9(22), 2858. https://doi.org/10.3390/math9222858

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