MPEMathematical Problems in Engineering1563-51471024-123XHindawi10.1155/2021/55594755559475Research ArticleReliability Analysis of Dynamic Fault Tree Based on Binary Decision Diagrams for Explosive Vehiclehttps://orcid.org/0000-0002-2095-7120JiangGuang-Jun12LiZong-Yuan1https://orcid.org/0000-0002-3386-5831QiaoGuan12ChenHong-Xia12LiHai-Bin1SunHong-Hua12YazdiMohammad1School of Mechanical EngineeringInner Mongolia University of TechnologyHohhot 010051Inner MongoliaChinaimut.edu.cn2Inner Mongolia Key Laboratory of Advanced Manufacturing TechnologyHohhot 010051Inner MongoliaChina20211042021202126220212432021293202110420212021Copyright © 2021 Guang-Jun Jiang et al.This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Dynamic fault tree is often used to analyze system reliability. The Markov model is a commonly used method, which can accurately reflect the relationship between the state transition process and the dynamic logic gate transfer in the dynamic fault tree. When the complexity or scale of system is increasing, the Markov model encountered a problem of state space explosion leading to increase troubles. To solve the above problems, a modular approach is needed. Based on the modular approach, a hybrid fault module was researched in this paper. Firstly, the stackable fault subtree containing complex static/dynamic logic gate is transformed into four common combinational logic gates through preprocessing of the dynamic gate in the module. Then, the complexity of the model was reduced by incorporating four common combinational logic gates and using the binary decision graph to solve variable ordering in the calculation of failure probability of static subtree. Moreover, the calculating process of complex mixed logic gate fault tree can be simplified. An example of the ammonium nitrate/fuel explosive production system for BCZH-15 explosive vehicle was used to verify the feasibility of the presented method.

National Natural Science Foundation of China71761030Inner Mongolia University of TechnologyYJG2017013
1. Introduction

Fault tree analysis (FTA) is a common method for the reliability modeling and evaluating large safety-critical systems as discussed in [1, 2]. The classic fault tree is a static tool, and its primary objective is to improve the safety for the probabilistic safety assessment as discussed in [3, 4]. Dynamic fault tree is an extension of common fault tree, which adds many special logic gates to express the existing problems in specific systems. Mi et al.  studied system action of mechatronics systems through taking the advantages of the dynamic fault tree (DFT). Fault tree analysis is a risk assessment method widely used in the process industry. Yazdi et al.  improved the traditional fault tree analysis, mainly for the study of uncertainty handling in qualitative and quantitative risk assessment procedures. The Markov model is appropriate for modeling complex systems that have traits of timing, sequencing, repair, redundancy, and fault tolerance as discussed by Ericson . However, the Markov model becomes unwieldy as the system size increases. It is difficult to construct and solve Markov models for the large systems as discussed in . Li  used the Bayesian network method to avoid dealing with Markov state explosion problem. Li et al.  used multiple methods such as multisource information fusion and Bayesian networks in the process of studying the reliability of offshore wind turbines. Dugan et al.  utilized dynamic gates to discern independent dynamic subtrees and used the different Markov model for each subtree through the process of modularization. As discussed by Jiang et al. , in solving the problem of reliability analysis and analysis of dynamic systems, fuzzy set theory is used. Huang et al.  proposed an improved decomposition scheme, where the dynamic subtree can be further modularized (if there exist some independent subtrees). A huge fault tree may contain dynamic subtrees and static subtrees. Use the binary decision diagram (BDD) method to calculate the failure rate of the top event of the fault tree. BDD was applied to efficiently realize synthesis of reversible functions as discussed by Mathias Soeken et al. . The modularization for large fault trees before analysis can improve computational efficiency. The approach of dynamic innovative fault tree (DIF tree) as discussed in  exploited the modularization by combining static and dynamic analysis. Ge et al.  analyzed DFT using improved the sequential binary decision diagram (SBDD) method. Bandyopadhyay et al.  improved a design technique, where BDDs are used to construct optical circuits.

However, after the modularization, the applicability of the model should be considered for relatively complex systems. Zhang et al.  proposed a generalized dynamic reliability model to calculate reliability of system under complex load. In dynamic systems, the semi-Markov process (SMP) was applied to solve the problem that the lifetime of components follows nonexponential distributions as discussed by Li et al. . Abimbola et al.  presented a risk assessment methodology for analyzing the safety of critical components based on BN as discussed in . This paper studies how to use modular methods to solve system reliability problems.

In recent years, the focus on field mixed explosive vehicle research has been performed to improve or design new electrical control systems to enhance the safety of equipment; to study the advantages of equipment in special geographic locations; to improve the blasting efficiency by changing the composition of the mixed explosive. There are only a few literatures on the reliability of the equipment, and only the reliability of the equipment can meet the operating requirements and can the safety, superiority, and efficiency of the equipment be demonstrated. With the dynamic characteristics becoming more and more obvious, it is urgent to evaluate equipment reliability considering the dynamic characteristics.

The remainder of the paper is organized as follows: Section 2 introduces division method of fault tree modular and illustrates how to preprocess the dynamic fault tree’s logical gates. Section 3 investigates the static subtree BDD theory solution method. The dynamic fault tree Markov theory solution method is described in Section 4. The proposed method was applied to calculate the reliability of the explosive production system for the BCZH-15 explosive vehicle in Section 5. Finally, the conclusion of the article is presented.

2. Fault Tree Modular Analysis Method2.1. Modular Division for Fault Tree

Rauzy et al.  presented the fault tree depth traversal method, which can modularize the fault tree. This method is concise and highly effective to deal with fault tree modularization problem. The deep traversal node method is used to divide module for the fault tree, which is easy to discover modules and zone. The procedure is as follows:

Calculate the traverse node for the fault tree and list the results

Sort the root node, leaf nodes, and the intermediate nodes and perform depth-first leftmost traversal for them

Collect the information of the internal event V for the third traversal node and collect the first and the last period of each subtime to provide strong evidence for modularizing the event.

Module division for the fault tree is shown in Figure 1, and the results are shown in Tables 1 and 2.

The fault tree.

First traversal result.

S12345678
VTM1M2M4A1A2M4M5
S910111213141516
VA3A4M5M2M3M5M6A5
S171819202122
VA6M6M3M1A7T

Final traversal result.

Module nameTM1M2M3M4M5M6
Min23485916
Max2119141861017
Y/NYYNNYYY

In Table 1, S indicates that the fault tree executes the left traversal node rule and then passes through the number of the node in turn, and V represents the name of each node in the fault tree. The sequence is given in Table 1.

Table 1 can summarize the target node numbers that each module can reach after executing the leftmost traversal node method, and then push out whether the modules are independent. In Table 2, Min is the node number that the module first arrives, and Max is the last module. Y means that the module is a separate module, and N means a nonindependent module. The division results are shown in Table 2.

From the operations given in Table 2, it can be concluded that the part of the fault tree modulo can be divided into T, M1, M4, M5, and M6.

2.2. Preprocessing of Logic Gate

After modularizing the fault tree, it is necessary to further determine the module state (dynamic or static). The BDD method is used to calculate the subtree for the static submodule. The Markov model is applied to mold for the dynamic submodule. Since the system is divided into small modules, state explosion problem will not occur during the model processing. Due to logic gate nesting, the static submodule may contain dynamic submodule after the module division in a complex fault tree. In this case, the corresponding algorithm cannot be applied to the fault tree submodule. Preprocessing of dynamic logic gate is required for the basic fault tree, and the complex hybrid dynamic logic gate is transformed into an easy solution model, which can greatly simplify the subsequent calculation.

2.2.1. Processing of Simple Logic Gate

Pretreatment of the logic gate includes the following: copretreatment of AND gate by OR gate, pretreatment of AND gate or PAND gate, pretreatment of two kinds of PAND gate stacking, and pretreatment of FDEP gate (Figures 25).

AND gate and OR gate preprocessing.

PAND gate and OR gate preprocessing.

Priority and gate preprocessing.

FDEP gate preprocessing.

The above is the basic preprocessing process of logic gates.

3. BDD Subtree Calculation

The BDD method is generally applied to solve the problem of the static fault tree. The BDD method plays a great role in promoting the analysis of the static subtree. The BDD method was proposed by the American scholar Akers in 1978, as discussed in . The method can transform the fault tree into a BDD only containing bottom events. The advantage of the BDD method is that the fault tree can be analyzed only using Weibull function expression without depending on intermediate events. Based on the BDD, disjoint Weibull function can be written down directly, and the minimum cut set can be obtained for qualitative analysis. The calculating amount of the method is linearly related to the scale of the graph. Therefore, it is an effective way to solve the problem that the calculating amount of the traditional FTA method increases exponentially with the scale of the fault tree.

3.1. Static Fault Tree Transformation

In the BDD conversion process, the main operation is based on the if-then-else (ite) structure of the Shannon decomposition, shown in the following equations:(1)fX=xifi+xi¯f0,(2)f1=f1x1,x2,,xi1,1,xi+1,,xn,f0=f0x1,x2,,xi1,0,xi+1,,xn,where the parent node is xi and the child nodes are f0 and f1.

3.1.1. Validation

Simple coding and index ordering are performed using ite for a simple fault tree. The influence of different index sorting on BDD generation is analyzed.

The fault tree is shown in Figure 6. The following equations represent the relationships among events:(3)T=x1+G1,G1=x2x3,=itex2,1,0itex3,1,0,=itex2,itex3,1,0,0,(4)T=itex1,1,0+itex2,itex3,1,0,0,=itex1,1,itex2,itex3,1,0,0.

The fault tree transformation diagram.

The BDD conversion graph of the fault tree is drawn by the above relationship.

3.1.2. Impact Analysis of Variable Ordering Results

At present, the optimal ordering of BDD is still under study. The mainstream BDD optimization algorithms are precise sorting algorithm and dynamic heuristic algorithm. The optimal index ranking is shown in Figure 6). When the variable order becomes index (x2) > index (x3) > index (x1), BDD would change (Figure 7).

Reordering BDD.

Since the different ordering, the original six sides change to twelve sides. The calculation process is increased from two sets of equations to four groups, and variable ordering is especially critical when dealing with more complex models. It is shown that a good variable ordering is related to the complexity of the BDD.

3.2. BDD Probability Calculation

After sorting the variables, the path pointing to 0 or 1 can be obtained, where the point 1 represents the top event occurring and the point 0 represents not occurring.

3.2.1. Calculation of Top Event Probability

Determining all paths point to 1, the path is recorded as Pii=1,2,3,,m, where m is the number of paths. Pi is shown in the following equation:(5)Pi=i=1nixij',xij'xij,xij¯,xijx1,x2,,xn.

Disjoint expressions of the fault tree can be represented as follows:(6)PT=Pi=1mPi=j=1mj=1nixij'.

The probability of occurrence of the top event can be calculated using probability formula of the mutual exclusion event as below:(7)PT=Pi=1mPi=j=1mj=1niPxij'.

The probability of failure of the bottom event can be expressed as Pxij.

3.3. A Case of Study

The above content is verified by an example. A, B, C, D, and E are bottom events (Figure 8). The probability of failure of the bottom event is shown in Table 3.

Modular fault tree.

Bottom event failure probability.

Bottom event codeABCDE
Failure probability0.20.10.20.30.15

After determining the probability of the leaf node, the fault tree can be modularly divided. The fault tree is divided into three parts: M1, M2, M3, where the submodules are M1 and M2. The submodule M1 is converted into BDD, as shown in Figure 9. Combining the transfer graph with the failure probability of variables B and C, probability of occurrence of module M1 can be obtained, which is 0.28.

Model M1 and model M2.

The probability of occurrence of module M2 is calculated based on the following BDD (Figure 9). Combined with failure probability of variables D and E, the probability of occurrence of module M2 can be obtained, which is 0.405.

After obtaining the probability of occurrence of the modules M1 and M2, probability of occurrence for M3 can be further calculated (Figure 10). The BDD model of the fault tree, defined as module M3, is a combination of M1, A3, and M2.

Model M3.

Obtain the probability of M1, M2, and A, and the probability of M3 can be calculated as 0.657.

4. Solution of Markov Model by Dynamic Fault Tree4.1. Markov Model Theory

Continuous-time Markov decision processes provide a very powerful mathematical framework to solve widely used decision problems, as discussed by Bartocci . The process can be defined as the following random process:(8)Xt,tT,ift0<t1<tn<tn+1T,PXtn+1xn+1|Xt0=x0,Xt1=x1Xtn=xn,=PXtn+1xn+1|Xtn=xn,where A is a random variable event, which occurs before time tn, and event B is an event that occurs afterwards. The Markov chain is shown in the following equation:(9)PB|A,Xtn=xn=PB|Xtn=xn,PAB|Xtn=xn=PA|Xtn=xnPB|Xtn=xn.

The sequence is shown as follows:(10)PXn+1=in+1|X0=i0,X1=i1Xn=in,=PXn+1=in+1|Xn=in.

4.2. Markov and Dynamic Logic Gate Transformation

In order to make the fault tree model better deal with the sequential logical relationship between the parts in the dynamic system, Dugan et al.  proposed the dynamic gate to represent the dynamic working relationship of the system in 1992. The Markov chain can represent the conversion relationship between system states using graphics. This graphical structure provides a visual tool to solve problems. Dynamic fault trees (DFTs) introduce four basic (dynamic) gates: the priority AND (PAND), the sequence enforcing (SEQ), the standby or spare (SPARE), and the functional dependency (FDEP), as discussed by Dugan et al. . The Markov transformation process of specific dynamic logic gates will be introduced.

4.2.1. Dynamic Logic Gate Transformation

The process of converting the four common dynamic logic gates into Markov chain will be described. The four common dynamic logic gates are PAND gate, SEQ, SPARE, and FDEP, and the relevant transition diagrams are shown in Figures 1114.

PAND logical gates transform to MC.

SEQ logical gates transform to MC.

CSP logical gates transform to MC.

FDEP logical gates transform to MC.

4.3. Quantitative Analysis of Markov chain

The quantitative analysis method of the Markov model generally consists of five steps: first, define the system state; second, the transition probability matrix construction; third, solve the spherical transition probability matrix; fourth, solve the differential equations; fifth, solve the fault probability state at any time.

The system state set, the fault state set, and the working state set are defined as follows: M=0,1, F=1, and W=0.

When defining the random process, the corresponding time point t should be set as shown below:(11)Xt=1,system failure,0,system working.

Defining A=P_E, where A is the transition matrix and E is equal to the identity matrix of P. The equation of state is shown below:(12)P=AP,P0=E.

The derivative column of the state probability and the column vector of the probability derivative are, respectively, expressed as Pand P. Equation (13) is C–K equation:(13)Pijt=qijpijt+k=jqikpkjt,pijt=1,i=j,pijt=0,ij.

In the process of the Markov chain transfer, the complexity of state transition usually increases with the increase of chain length. The transition probability from i to j can be expressed as Pij as follows:(14)Pt=P00tP01tP0NFtP10tP11tP1NFtPNF0tPNF1tPNFNFt.

The transfer probability matrix can be written as follows:(15)Q=λ0,1+λ0,NFλ0,1λ0,NF000000,when the chain length is 1, the state transition process can be represented as follows:(16)PT1=P01t=λ0,1λ0,1+λ0,NF1eλ0,1+λ0,NFt,=λ0,11λ0,1+λ0,NF1λ0,1+λ0,NFeλ0,1+λ0,NFt,λ0,10λ0,NF0.

The state transition process of chain length 2 can be deduced by the formula of chain length 1. When the event from 0 to j requires two transfer processes, this process is defined as T2, and PT is the ultimate transition probability of each edge on T2. This process can be seen as follows:(17)PT2=P02t=λ0,1λ1,21λ0,1+λ0,NF+λ1,2+λ1,NF1λ0,1+λ0,NFλ0,1λ0,NF+λ1,2+λ1,NFeλ0,1+λ0,NFt1λ1,2+λ1,NFλ1,2λ1,NF+λ0,1+λ0,NFeλ1,2+λ1,NFt.

The following equation represents the state transition of chain length is n:(18)PTn=P0nt=i=1nλi1,ii=1n1λi1,i+λi1,NFi=1n1λi1,i+λi1,NFj=1,jinλi1,iλi1,NF+λj1,j+λj1,NFeλi1,i+λi1,NFt,where λi1,i>0,λi1,NF>0.

5. Fault Tree Modular Analysis

The main work of the pharmaceutical system of explosive vehicle is to manufacture explosives. The working environment of the explosive vehicle is harsh, and the failure rate is high. Through existing data, the failure probability of bottom event was obtained for the explosive vehicle.

5.1. Fault Tree Establishment and Module Division

Figure 15 shows the dynamic fault tree. It can be seen that the fault tree established includes several static and dynamic modules. The system is mainly divided into four types: flow display module, storage tank module, pump body module, and sensitizing agent system module. The storage tank module is mainly responsible for storing various liquids to ensure that the catalysts used in the preparation process are kept in a stable state. The main function of the chemical agent module is to filter to ensure that no impurities are mixed in during the configuration process. The pump body module controls the stable output of the ingredients during the preparation of explosives. The flow display module ensures that the flow output process can be observed by the outside world.

The model of dynamic fault tree.

5.1.1. Preprocessing the Fault Tree

The FDEP gate exists under the OR gate of module N1. The FDEP gate transformation process is as follows: the subordinate modules of N2 are all static gates. The module M1 under N2 can be combined due to a static module. Fault tree after preprocessing is shown in Figure 16.

Dynamic fault tree module division.

The fault tree can be divided into four basic modules S1, S2, S3, and S4. In module S1, there is a submodule C consisting of static logic gates, which can be divided into static module. Modules S2 and S3 are static modules. Module S4 is divided into dynamic submodule because its subtree root node is the spare parts gate.

5.2. Fault Tree Modular Calculation

According to BDD, its nodes are firstly divided. Sorting the bottom events uses a heuristic sorting method and defined the order indexx1>indexx2>indexx3. Get BDD (Figure 17).

The binary decision diagrams.

For the transformation of the static module S2, the above method should be used. First, the heuristic algorithm is applied to determine the subunit sequencing:(19)in  de  xx7>in  de  xx8>in  de  xx9>in  de  xx10>in  de  xx11>in  de  xx12>in  de  xx13.

The BDD of module S2 is shown in Figure 18.

The binary decision diagrams.

Similarly, module S3 was converted into BDD, where the sorting is defined as(20)in  de  xx15>in  de  xx16>in  de  xx17>in  de  xx18.

The BDD of transformation is shown in Figure 19.

The binary decision diagrams.

The module S4 contains the submodules C4 and C5, where C5 is the spare part of C4. BDD was established to solve the probability of C4 (Figure 13). The order has no succession, and Markov chain is shown in Figure 20.

The BDD and Markov chain.

The failure probability of parts in the bottom event is shown in Table 4. The following data sources are from the maintenance records of the BZ15 heavy-duty oil ammonia explosive vehicle and related experimental data in the factory.

Bottom event failure probability.

EncodingEventProbability (104h1)
X1Line interface failure0.5
X2Fuel flow-meter failure0.4
X3Sensitizer flow-meter failure0.2
X4Solenoid valve failure0.3
X5Connection line failure0.2
X6System display failure0.1
X7The leakage of fuel tank0.1
X8The leakage of oil pipe joints0.5
X9Oil pipe failure0.1
X10Oil filter element failure0.8
X11Oil-pressure meter failure0.2
X12One-way valve failure0.2
X13Stator failure0.45
X14Bearing failure0.15
X15Emulsion matrix pump failure0.3
X16Flexible pipe failure0.7
X17Flexible pipe joint failure0.2
X18Butterfly valve failure0.4
X19Sensitizer pump failure0.1
X20Water tube failure0.5
X21Sensitizer filter element failure0.2
X22Sensitizer tank failure0.2

The probability of failure of each module at 1000 hours is obtained by the bottom event probability, as shown in Table 5.

Failure probability of each module.

EncodingModelProbability
CFuel flow module0.1008
S1Flow display module0.1537
S2Tank module0.2261
S3Pump body module0.1513
S4Sensitizing agent system module0.3333

After obtaining the data failure probability of each module, the data are further integrated to obtain failure probability of the top event.

According to the integration, the reliability of the system after one thousand hours of running time can be obtianed. The integration diagram is shown in Figure 21.

Integration diagram.

Finally, combining the failure probability of each module to evaluate the system reliability and failure rate can be obtained as 0.6292 and 0.3708, respectively.

6. Discussion

The Markov model is usually used to solve the problem of dynamic fault tree. Complex systems are difficult to build using traditional Markov models. Due to the dynamic fault tree containing static subtrees, the BDD can be used to solve the problem of static subtree. The logic gate combination is simplified using the pretreatment method. For the simplified dynamic fault tree , BDD is used to solve the problem. The complexity of the solution is closely related to the construction of the binary decision graph model. The complexity of model construction is affected by the ordering of variables. The complex structure of the fault tree is improved, and the difficulty of sorting the BDD variables is reduced by preprocessing the logic gate. Taking the subtree N1 in Figure 16 as an example, the FDEP gate is connected under the OR gate. If the variables are sorted directly without preprocessing, it is difficult to reasonably sort due to the various sorting of dynamic logic gate. The preprocessed module is shown in Figure 17. The preprocessed module S1 only contains two kinds of static gates. Sorting the variables is becoming easier, and difficulty of solving is reduced. The binary decision method can be directly used to solve the solution, and the result is as shown in Table 5. Using this method, the failure probability of each subtree can be obtained at the same time, and the dynamic fault tree can be rigorously analyzed.

The information obtained from the current research results is still very limited. Because of its influence by many factors, the accuracy of the failure probability still needs further research. Therefore, the subsequent researches should go along diversification. Since the amount of information that can be obtained using only one data source is small, in the future, a multisource information fusion method will be used to more accurately evaluate the system reliability.

7. Conclusion

This article introduces how to use the modular approach to solve the probability density of basic events in the dynamic fault tree when the system has certain complexity and multiple subsystems. The module pretreatment method is adopted to simplify the submodule of the fault tree due to the excessive mixture and overlap of logic gates. This modular approach has two advantages: first the variables ordering for the module becomes easy to reduce the complexity of BDD and solve difficulty by pretreatment, and second the failure probability of the subsystem can be obtained for evaluating the reliability of the subsystem. The feasibility of the method was verified by analyzing dynamic fault tree of explosive production system for the BCZH-15 explosive vehicle. At present, due to the difficulty of individual data collection in the data collection process and the lack of a large amount of experimental data support, the failure rate of the collected products is not accurate enough. It may have an impact on the accuracy of the result evaluation. In the future, when the data are insufficient, fuzzy theory will be used to define the failure rate of the product, and the rationality of the final output result will be determined through the comparison of expert experience.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

There are no conflicts of interest.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China under the contract no. 71761030 and the Graduate Teaching Program of Inner Mongolia University of Technology under the contract no. YJG2017013.

HenleyE. J.KumamotoH.Reliability Engineering and Risk Assessment1981Englewood Cliffs, NJ, USAPrentice-HallLevesonN. G.Safeware: System Safety and Computers1995Reading, MA, USAAddison-WesleyRenY.DuganJ. B.Optimal design of reliable systems using static and dynamic fault treesIEEE Transaction on Reliability199847234244SiuN.Risk assessment for dynamic systems: an overviewReliability Engineering and System Safety1994431437310.1016/0951-8320(94)90095-72-s2.0-0028197894MiJ.LiY. F.YangY. J.PengW.HuangH. Z.Reliability assessment of complex electromechanical systems under epistemic uncertaintyReliability Engineering and System Safety2016152115YazdiM.KabirS.WalkerM.Uncertainty handling in fault tree based risk assessment: state of the art and future perspectivesProcess Safety and Environmental Protection20191318910410.1016/j.psep.2019.09.0032-s2.0-85072265041YazdiM.NikfarF.NasrabadiM.Failure probability analysis by employing fuzzy fault tree analysisInternational Journal of System Assurance Engineering and Management2017821177119310.1007/s13198-017-0583-y2-s2.0-85015420795YazdiM.NedjatiA.ZareiE.AbbassiR.A novel extension of DEMATEL approach for probabilistic safety analysis in process systemsSafety Science202012111913610.1016/j.ssci.2019.09.006EricsonC. A.Hazard Analysis Techniques for System Safety2015Hobokrn, NJ, USAJohn Wiley & SonsTorres-EcheverriaA. C.Modeling and optimization of proof testing policies for safety instrumented systemsReliability Engineering and System Safety200994838854LisnianskiA.Extended block diagram method for a multi-state system reliability assessmentReliability Engineering and System Safety20079216011607DutuitY.InnalF.RauzyA.SignoretJ.-P.Probabilistic assessments in relationship with safety integrity levels by using Fault TreesReliability Engineering and System Safety200893121867187610.1016/j.ress.2008.03.0242-s2.0-49949089666LiY. F.MiJ.LiuY.YangY. J.HuangH. Z.Dynamic fault tree analysis based on continuous-time Bayesian networks under fuzzy numbersProceedings of the Institution of Mechanical Engineers, Part O, Journal of Risk and Reliability20052296530541LiH.DiazH.Guedes SoaresC.A developed failure mode and effect analysis for floating offshore wind turbine support structuresRenewable Energy202116413314510.1016/j.renene.2020.09.033LiH.TeixeiraA. P.Guedes SoaresC.A two-stage Failure Mode and Effect Analysis of offshore wind turbinesRenewable Energy20201621438146110.1016/j.renene.2020.08.001LiH.SoaresC. G.HuangH. Z.Reliability analysis of a floating offshore wind turbine using Bayesian NetworksOcean Engineering2020217107827DuganJ. B.BavusoS. J.BoydM. A.Dynamic fault-tree for fault-tolerant computer systemsIEEE Transaction on Reliability1992413363376DuganJ. B.SullivanK. J.CoppitD.Developing a low-cost high-quality software tool for dynamic fault-tree analysisIEEE Transactions on Reliability2000491495910.1109/24.8555362-s2.0-0034155553MeshkatL.DuganJ. B.AndrewsJ. D.Dependability analysis of systems with on demand and active failure modes using dynamic fault treesIEEE Transaction Reliability200251324025110.1109/tr.2002.10115312-s2.0-0036599709JiangG.-J.GaoL.Reliability analysis of martial arts arena robot systems based on fuzzy set theoryJournal of Mechanical Science and Technology201832115069507710.1007/s12206-018-1003-12-s2.0-85057984015HuangC.-Y.ChangY.-R.An improved decomposition scheme for assessing the reliability of embedded systems by using dynamic fault treesReliability Engineering and System Safety200792101403141210.1016/j.ress.2006.09.0082-s2.0-34250156732SoekenM.TagueL.DueckG. W.DrechslerR.Ancilla-free synthesis of large reversible functions using binary decision diagramsJournal of Symbolic Computation20167312610.1016/j.jsc.2015.03.0022-s2.0-84939651683DuganJ. B.DIF tree: a software package for the analysis of dynamic fault tree modelsProceedings of the Reliability and Maintainability SymposiumJanuary 1997Philadelphia, PA, USA6470GeD.LinM.YangY.ZhangR.ChouQ.Quantitative analysis of dynamic fault trees using improved Sequential Binary Decision DiagramsReliability Engineering and System Safety201514228929910.1016/j.ress.2015.06.0012-s2.0-84934964324BandyopadhyayC.DasR.WilleR.DrechslerR.RahamanH.Synthesis of circuits based on all-optical mach-zehnder interferometers using binary decision diagramsMicroelectronics Journal201871192910.1016/j.mejo.2017.11.0082-s2.0-85034829660ZhangX.GaoH.HuangH.-Z.LiY.-F.MiJ.Dynamic reliability modeling for system analysis under complex loadReliability Engineering and System Safety201818034535110.1016/j.ress.2018.07.0252-s2.0-85051279543LiX.-Y.HuangH.-Z.LiY.-F.Reliability analysis of phased mission system with non-exponential and partially repairable componentsReliability Engineering and System Safety201817511912710.1016/j.ress.2018.03.0082-s2.0-85043987673AbimbolaM.KhanF.KhakzadN.ButtS.Safety and risk analysis of managed pressure drilling operation using Bayesian networkSafety Science20157613314410.1016/j.ssci.2015.01.0102-s2.0-84925067459RauzyY. D.A linear-time algorithm to find modules of fault treeIEEE Transaction on Reliability1996453AkersS. B.Binary decision diagramsIEEE Transaction on Computers1978C-27509516BartocciE.BortolussiL.BrázdilT.MiliosD.SanguinettiG.Policy learning in continuous-time Markov decision processes using Gaussian ProcessesPerformance Evaluation20171168410010.1016/j.peva.2017.08.0072-s2.0-85029590224ReedS.An efficient algorithm for exact computation of system and survival signatures using binary decision diagramsReliability Engineering and System Safety201716525726710.1016/j.ress.2017.03.0362-s2.0-85018500963ToppilaA.SaloA.Binary decision diagrams for generating and storing non-dominated project portfolios with interval-valued project scoresEuropean Journal of Operational Research20172601244254MiJ.LiY.-F.PengW.HuangH.-Z.Reliability analysis of complex multi-state system with common cause failure based on evidential networksReliability Engineering and System Safety2018174718110.1016/j.ress.2018.02.0212-s2.0-85042370212LiH.HuangH.-Z.LiY.-F.ZhouJ.MiJ.Physics of failure-based reliability prediction of turbine blades using multi-source information fusionApplied Soft Computing20187262463510.1016/j.asoc.2018.05.0152-s2.0-85047067866