# Fuzzy Dynamic Reliability Models of Parallel Mechanical Systems Considering Strength Degradation Path Dependence and Failure Dependence

Authors: Peng Gao; Liyang Xie
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Fuzzy dynamic reliability models of mechanical parallel systems with respect to stress parameters and strength parameters are developed in this paper. Strength degradation path dependence (SDPD) and failure dependence of components in the system are two main problems to be addressed in developing fuzzy dynamic reliability of mechanical systems, which are taken into account in the proposed reliability models. In addition, the SDPD sensitivity function and the failure dependence sensitivity function are defined to analyze the influences of the fuzzy characteristics of stress on the effects of SDPD and failure dependence of mechanical components on dynamic reliability. The bolted joint systems for connection between launch vehicle and satellite are chosen as illustrative examples to demonstrate the proposed model. Moreover, Monte Carlo simulations are carried out to validate the effectiveness of the proposed models. The results show that SDPD and failure dependence of components have significant impacts on fuzzy dynamic reliability of mechanical systems. Besides, the influences of the fuzzy characteristics of stress on SDPD sensitivity and that on failure dependence sensitivity are different.


---

## Body

## 1. Introduction

Over the last few decades, reliability models of mechanical systems have been well reported in the probability context [1–5]. The basic idea for reliability analysis of mechanical systems is to calculate the system reliability according to the component reliability in the system and structure function (SF). For mechanical components, the well-known stress-strength interference (SSI) model is the most important theory for reliability assessment, where both stress and strength are modeled as random variables. In the SSI model, the failure of mechanical components is defined as the probability that the stress is higher than strength as shown in Figure 1. The interference area between the probability density function (pdf) curve of stress and the pdf curve of strength is an indication of the possibility that the failure of the component might occur.Figure 1
Schematic SSI model.However, for mechanical systems in aerospace products, failures are relatively rare events. It is a costly and difficult undertaking to collect enough data to calculate the probability of failure. Besides, due to the complexity of the operational environment, it is difficult to obtain a large number of samples to precisely describe the statistical characteristics of stress caused by environmental load on the mechanical systems in aerospace products, especially in the preliminary design stage where many factors have to be determined by the judgment of engineers. In this case, it can only increase the uncertainty in reliability estimation to adopt the imprecise pdf of stress. Alternatively, fuzzy theory can be used to deal with this problem by modeling stress in an analytical form via integrating the limited available data and expert judgments.As a matter of fact, fuzzy variables are always used to deal with the uncertainty brought by insufficient information. In this paper, we model the stress at each load application as a fuzzy variable. We assume that the uncertainty of stress comes from the insufficient information about the stress. When the amount of available information increases, which is reflected via the enhancement of membership degree, the uncertainty of stress decreases and the stress tends to be a deterministic number.Since Zadeh introduced the fuzzy theory in 1965 [6], great efforts have been made to develop fuzzy reliability models. Cai et al. classified the fuzzy reliability models into three basic classes: the robust reliability models, the posbist reliability models, and the posfust reliability models [7, 8]. Besides, Cai et al. proposed fuzzy reliability models of various typical systems based on the possibility assumption and the fuzzy state assumption [9]. Wu discussed the relationship between the fuzzy reliability of a system and the fuzzy reliability of its components via the proposed fuzzy reliability models [10]. Utkin proposed the method for fuzzy reliability analysis and fuzzy reliability optimization [11]. Huang et al. presented an approach for reliability estimation by using the Bayesian method based on fuzzy lifetime data [12]. Bing et al. extended the conventional SSI model to the fuzzy stress—random strength interference model to analyze the fuzzy reliability of mechanical components [13]. These fuzzy reliability models establish the basis and framework for reliability assessment in the possibility context. However, for mechanical systems, it is difficult to calculate reliability by directly employing these models because of the problems as follows.(1)
Strength degradation exists commonly for mechanical components in practical engineering. The interaction between the stress and the strength is a dynamic process. Most of current fuzzy reliability models of mechanical components and mechanical systems are static, which are only suitable for the reliability calculation at a specific moment as shown in Figure2. As mentioned in the authors’ former research in [14], when considering the strength degradation process of mechanical components, the strength degradation path dependence (SDPD) is a main problem to be addressed in developing dynamic reliability of mechanical components, which also exists in establishing fuzzy reliability models. As a matter of fact, fuzzy dynamic reliability models of parallel mechanical systems considering SDPD is seldom reported. Besides, the fuzzy characteristics of stress have great influences on SDPD, which is also seldom reported in current literatures. In this paper, impacts of fuzzy characteristics of stress on the effects of SDPD on reliability will be analyzed via the proposed fuzzy dynamic reliability models.(2)
Due to the unique operating mode of mechanical product, failure dependence exists commonly in mechanical systems and the components, under the stress with a common source of load, in the system are statistically dependent on each other. Fuzzy dynamic reliability models of parallel mechanical systems with respect to stress and material parameters, which consider failure dependence between different components in a system, are seldom reported. Correspondingly, the influences of fuzzy characteristics of stress on failure dependence need further investigation. In this paper, the sensitivity analysis of dynamic system reliability with respect to fuzzy parameters of stress will be performed in order to present how the fuzzy characteristics of stress affect the failure dependence behavior of components in mechanical systems.Figure 2
Schematic stress process and strength process.The remainder of this paper is organized as follows. In Section2, fuzzy dynamic reliability models of parallel mechanical systems are developed, which take SDPD and failure dependence into consideration. In Section 3, Monte Carlo simulations are carried out to verify the reliability models. Furthermore, the influences of fuzzy characteristics of stress on both the effects of SDPD and the effects of failure dependence on reliability are analyzed via the numerical examples.


## 2. Fuzzy Dynamic Reliability Models of Parallel Mechanical Systems

Denote the pdf of a random variablex by f(x) and the membership function (MF) of x belonging to a fuzzy set X~ by μX~(x). Then, the fuzzy probability that x belongs to X~ can be calculated as [6] (1)Px=∫-∞∞μX~xfxdx.Meanwhile, the fuzzy probability in (1) can also be calculated by employing the α-cut of fuzzy set as follows:(2)Px=∫01∫X~αLX~αUfxdxdα.In fact, the expression of fuzzy probability in (2) provides the basis for reliability evaluation with multiple format of variables, such as the random variables and fuzzy variables, involved in the calculation. As mentioned above, in some cases, it is difficult to obtain the accurate pdf of stress on components in mechanical systems because of insufficient load history samples, while most pdf of material parameters can be acquired from handbooks for mechanical design or experiment. Therefore, in this paper, the stress is modeled as a fuzzy variable and strength is modeled as a random variable. As a matter of fact, the triangular membership function (TMF) is always adopted to express the fuzzy characteristics of stress and the parameters in TMF of stress s can be evaluated by using the method in [13]. The mathematical expression of TMF of s can be given by(3)μss=s-ab-a,a≤s<b,1,s=b,s-cb-c,b≤s<c,0,otherwise.The TMF can also be expressed as a triplet(a,b,c) for descriptive convenience. In addition, appoint the MF to be equal to a specific value of α. Then, the elements with degree of membership larger than α constitute an α-cut. Through the α-cut, the computation under fuzzy uncertainty can be transformed into the computation under the interval uncertainty. Correspondingly, the calculation with respect to fuzzy numbers can be reduced to the calculation with respect to interval numbers. To consider the weight of elements on a specified interval and for the purpose of computational convenience in case of both random variables and interval variables existing, a certain distribution is sometimes assumed for the elements on the interval, such as the uniform distribution [15, 16] or the linear distribution [17]. As pointed out by Nguyen et al. when the distribution is unknown on a given interval, the maximum entropy approach is always used to determine an appropriate distribution for the elements on the interval. In the case where the uniform distribution is adopted, the entropy reaches its maximum value [18]. In this paper, we assign uniform distribution to the elements in an α-cut and denote the fuzzy stress at the ith load application by (ai,bi,ci). Then, the distribution function of the α-cut can be given by(4)Fiαsi=0,si<ai+bi-aiα,si-ai-bi-aiα1-αci-ai,ai+bi-aiα≤si<ci-ci-biα,1,si≥ci-ci-biα.According to the reliability theory of systems, a parallel system fails to operate when all the components in the system fail to work. In this paper, we assume that there arek identical components with the same structure and with the strength pdf of fr(r) in a parallel system. When the components in a parallel system are subjected to the same load, the stress of s with the pdf of fs(s) on each component is mutually identical. Then, the probability that the parallel system does not fail to work is equal to the probability that the maximum strength is higher than s, which can be expressed by using the order statistics theory as follows:(5)Prmax⁡>s=∫s∞k∫-∞rfrrdrk-1frrdr.Therefore, the reliability of the parallel system can be given by(6)Rparallel=∫-∞∞fss∫s∞k∫-∞rfrrdrk-1frrdrds.Equation (6) is similar to the SSI model in mathematical expression. Hence, the equivalent strength of the parallel system can be expressed as (7)fparallelr=k∫-∞rfrrdrk-1frr.In the failure mode of fatigue, strength could degrade under the application of random load which causes damage to the components in the system. The remaining strength is the function of the load application times and the magnitude of stress, which is generally expressed as follows [14]:(8)rn=r01-Dnγ,where r0 and γ are initial strength and material parameter, respectively. D(n) is the cumulative damage induced by the applied load. In addition, the S-N Curve, which is usually adopted to compute the lifetime of mechanical components under the stress with a specific magnitude, can be written in an analytical format as follows:(9)smN=C,where m and C are material parameters. When considering the effects of SDPD, the equivalent remaining strength under the stress at the level of α can be obtained according to the rule of damage accumulation and the rule of damage equivalence as follows [14]:(10)rαn=r01-∑i=1nci-aiC-1ci-ci-biαm+1-ai+bi-aiαm+1×m+11-αci-aiC-1∑i=1nγ.When the initial strength r0 is deterministic, the reliability at the ith load application can be given according to the distribution of the fuzzy stress on the interval of α-cut as follows:(11)Ri=0,rαn-1<ai+bi-aiα,rαn-1-ai-bi-aiα1-αci-ai,ai+bi-aiα≤rαn-1<ci-ci-biα,1,rαn-1≥ci-ci-biα.Therefore, the probability that the system operates normally when load applies for n times can be expressed as (12)Rαn=∏i=1nRi.When considering the distribution of the equivalent strength of the system characterized by the pdf shown in (7), the fuzzy reliability of the parallel system at the level of α can be obtained according to the total probability formula for continuous variables as follows:(13)Rαn=∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0.According to the decomposition theorem, the fuzzy dynamic reliability of the parallel system can be given by(14)Rn=∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0dα.Hence, the failure rate of the parallel system can be expressed as follows:(15)λn=∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0dα-∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1n+1Ridr0dα×∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0dα-1.It should be noted that (14) is derived with the effects of SDPD taken into consideration. For computational convenience, conventional dynamic reliability models always model the strength degradation process as a stochastic process with the distribution of strength at each load application determined by the specified stochastic process and calculate the reliability based on the strength at each load application. As pointed by the authors’ former study, it could result in large computational error due to neglecting the effects of SDPD [14]. When the fuzzy system reliability is calculated according to the distribution of component strength at each load application, which can be obtained from the distribution of initial strength and (10), the fuzzy reliability can be calculated as follows:(16)R1n=∫01∏i=1n∫-∞∞fiαsi1-∫-∞sfriridrikdsdα,where fiαsi, which can be obtained via derivation of Fiα(si) in (4) with respect to si, is the pdf of stress caused by the ith load application and friri is the pdf of strength at the ith load application.As a matter of fact, the fuzzy characteristics of stress are mainly determined by the parameters ofa, b, and c. b is the center parameter with the largest membership degree. When more information is available, which means more investment in experiments, the parameters of a, c will be closer to b and the fuzzy number tends to be a deterministic number. In order to analyze the influences of the parameters of a and c on the effects of SDPD on reliability, the SDPD sensitivity functions with respect to a and c are defined as follows:(17)ISDPDan=∂Rn-R1n∂a=∂∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0dα∂a-∂∫01∏i=1n∫-∞∞fiαsi1-∫-∞sfriridrikdsdα∂a,(18)ISDPDcn=∂Rn-R1n∂c=∂∫01∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0dα∂c-∂∫01∏i=1n∫-∞∞fiαsi1-∫-∞sfriridrikdsdα∂c.In addition, when the components in the system are mutually independent, the fuzzy reliability can be computed as(19)R2n=1-1-∫01∫-∞∞fr0r0∏i=1nRidr0dαk.In order to evaluate the impacts of the parameters ofa, c on failure dependence, we define the failure dependence sensitivity functions with respect to a and c as follows:(20)IFDan=∂Rn-R2n∂a=∂∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0∂a-k1-∫-∞∞fr0r0∏i=1nRidr0k-1×∂∫01∫-∞∞fr0r0∏i=1nRidr0dα×∂a-1,IFDcn=∂Rn-R2n∂c=∂∫-∞∞k∫-∞r0fr0r0dr0k-1fr0r0∏i=1nRidr0∂c-k1-∫-∞∞fr0r0∏i=1nRidr0k-1×∂∫01∫-∞∞fr0r0∏i=1nRidr0dα×∂c-1.


## 3. Illustrative Examples

The explosive bolts are comprehensively used for connection and separation in aerospace products. For example, the bolted joints for connection between launch vehicle and satellite are shown in Figure3. The redundancy design of the explosive bolt system largely improves the reliability of the success launch of satellites. In this section, the bolted joint system is chosen as an illustrative example to analyze the influence of the fuzzy characteristics of stress on SDPD and failure dependence of parallel mechanical systems. In the bolted joint system, the explosive bolts are mutually identical with the material parameters [14] and stress parameters listed in Table 1.Table 1
Stress parameters and material parameters of explosive bolts.
μ
(
r
0
) [MPa]
σ
(
r
0
) [MPa]
a [MPa]
b [MPa]
c [MPa]
m
γ
C [MPa2]
400
5
250
300
350
2
1
10
8Figure 3
The bolted joint system.In this section, three cases will be analyzed in order to investigate the following problems:(1)
verify the proposed model via Monte Carlo simulation;(2)
analyze the influences of the parameters ofa, c on the effects of SDPD on reliability;(3)
investigate the influences of the parameters ofa, c on the failure dependence.Case 1. Consider a system with three dependent components. The flowchart for the Monte Carlo simulation is shown in Figure 4. In the flowchart, k represents the number of components in the system. N is the total number of simulation trials. n is the number of load application. b is the total number of failure of the system in the N trials. It should be noted that in the Monte Carlo simulation, the strength degradation process is only dependent on the randomly generated stress and (8) for remaining strength calculation under deterministic stress. Therefore, the simulation is essentially identical with physical experiment and does not rely on any analytical dynamic reliability model. The comparison between the result of the Monte Carlo simulation and the result from the proposed models is shown in Figure 5.Figure 4
Flowchart of the Monte Carlo simulation for parallel subsystem.Figure 5
Reliability from Monte Carlo simulation and reliability from proposed models.From Figure5, it can be learnt that the result from the proposed reliability method shows good agreement with the reliability from Monte Carlo simulations. The reliability decreases rapidly with the load application.Case 2. Consider a system with three dependent components. The system reliability considering SDPD and the system reliability without SDPD taken into account are shown in Figure 6. In addition, the SDPD sensitivity with respect to a and c is shown in Figures 7 and 8, respectively.Figure 6
Reliability considering SDPD and reliability without SDPD taken into consideration.Figure 7
SDPD sensitivity with respect toa.Figure 8
SDPD sensitivity with respect toc.From Figure6, it can be seen that SDPD have great influences on fuzzy system reliability. Although it could facilitate the reliability calculation by directly using the strength distribution at each load application, reliability could be underestimated due to neglecting the effects of SDPD. Moreover, from Figures 7 and 8, it can be learnt that the left parameter a and right parameter c have different impacts on the effects of SDPD on reliability. The effects of SDPD on reliability are more sensitive to variation of c.For a deterministica, the local maximum of the error caused by SDPD arises at middle stage and the end of the operational duration of the system. The reliability calculated considering SDPD is higher than that without SDPD taken into account. Besides, from (17) and (18), we can learn that, for the SDPD sensitivity, a positive value means that the error increases with the increase of the parameter considered (a or c), while a negative value means that the error decreases with the increase of the parameter considered (a or c). Moreover, the reliability calculated considering SDPD is always larger than the reliability calculated without SDPD taken into account. Therefore, the SDPD sensitivity with a negative value means that the effects of SDPD on reliability are weakened. In mechanical design, engineers are more interested in the sensitive areas where a small variation in the parameter considered could lead to a rapid increase in the error due to neglecting SDPD. Thus, in this paper, we are concentrated on the characteristics of the positive peak value of SDPD sensitivity. From Figure 7, we can see that two evident positive peak values of the SDPD sensitivity appear in the middle stage of the system operational duration, which represent the most sensitive time instant. Therefore, in dynamic system reliability analysis, attention should be paid to the influences of the change in a on the effects of SDPD at these two sensitive “time instants” associated with their adjacent time interval. In addition, the positive peak value of the SDPD sensitivity comes earlier when a increases.For a deterministicc with a low value, only one positive peak value of the SDPD sensitivity appears in the middle stage of the system operational duration. However, when c is large, two positive peak values arise in the middle stage of the system operational duration with one peak value evidently larger than the other one. Therefore, we should pay attention to the variation of c at the moment where the positive peak value of the SDPD sensitivity appears in the reliability estimation when considering the effects of SDPD. Similar to the case in the sensitivity analysis with respect to a, the maximum of the SDPD sensitivity comes earlier when c increases.Case 3. In this section, we are aimed at analyzing the influences of the failure dependence of components on system reliability and influences of fuzzy parameters of stress on failure dependence sensitivity function. The system reliability considering failure dependence of components can be calculated according to (14), while the system reliability under the assumption that the components are independent of each other can be calculated according to (19). The failure dependence sensitivity function with respect to a and c can be obtained through (20). The reliability of a system with three dependent components and the reliability of a system with three independent components are shown in Figure 9. In addition, the failure dependence sensitivity with respect to a and c is shown in Figures 10 and 11, respectively.Figure 9
Reliability of dependent system and reliability of independent system.Figure 10
Failure dependence sensitivity with respect toa.Figure 11
Failure dependence sensitivity with respect toc.From Figure9, it can be seen that the failure dependence makes the parallel system less reliable. Thus, from (18) and (19), Figures 10 and 11, it can be seen that the negative peak value of the failure dependence sensitivity means that the error caused by failure dependence increases at the highest speed, while the positive peak value of the SDPD sensitivity means that the error decreases at the highest speed. Therefore, we are more concentrated on the negative peak value of the failure dependence sensitivity. A single negative peak value of the failure dependence sensitivity appears in the middle stage of the system operational duration in both the case of failure dependence sensitivity analysis with respect to a and the case of failure dependence sensitivity analysis with respect to c. We should pay attention to the variation of a and c at the moment where the negative peak value of the failure dependence sensitivity appears in the reliability assessment when considering the effects of failure dependence. Besides, the negative peak value of the failure dependence sensitivity comes earlier when a or c increases. In general, the effects of the failure dependence of components in the parallel system are more sensitive to the change of c.


## 4. Conclusion

In this paper, fuzzy dynamic reliability models of mechanical parallel systems with respect to stress parameters and strength parameters are established. In the proposed reliability models, the SDPD and failure dependence of components in a mechanical parallel system are taken into account. For computational convenience, strength distribution of components at each load application is always adopted in conventional dynamic reliability models, which is determined by the assumed stochastic process. However, the results show that SDPD has considerable influences on fuzzy dynamic reliability of parallel systems. In addition, the failure dependence of components in a system also significantly affects the system fuzzy dynamic reliability.By defining the SDPD sensitivity function and the failure dependence sensitivity function in this paper, it is found that the fuzzy characteristics of stress have great impacts on the effects of SDPD and failure dependence on fuzzy dynamic reliability. In general, the effects of SDPD and failure dependence of mechanical components in a parallel system on reliability are more sensitive to variation of the right parameter ofc than the left parameter of a. Moreover, the influences of a and c on SDPD sensitivity and that on failure dependence sensitivity are different. The proposed models provide the basis for analytically determining the moment when the peak values of SDPD sensitivity and failure dependence sensitivity appear, which are important to the dynamic reliability evaluation and reliability-based design of mechanical parallel systems. In addition, Monte Carlo simulation is carried out to validate the effectiveness of the proposed models.



---

*Source: 649726-2015-05-18.xml*