Next Article in Journal
Development of a Textile Nanocomposite as Naked Eye Indicator of the Exposition to Strong Acids
Next Article in Special Issue
A Robust Vehicle Localization Approach Based on GNSS/IMU/DMI/LiDAR Sensor Fusion for Autonomous Vehicles
Previous Article in Journal
Response Characteristics and Experimental Study of Underground Magnetic Resonance Sounding Using a Small-Coil Sensor
Previous Article in Special Issue
Fuzzy Risk Evaluation in Failure Mode and Effects Analysis Using a D Numbers Based Multi-Sensor Information Fusion Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of BN Lifetime Prediction of a System Based on Integrated Multi-Level Information

1
School of Reliability and Systems Engineering, Beihang University, Beijing 100191, China
2
Unmanned System Institute, Beihang University, Beijing 100191, China
3
Key Laboratory of Advanced Technology of Intelligent Unmanned Flight System of Ministry of Industry and Information Technology, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(9), 2123; https://doi.org/10.3390/s17092123
Submission received: 11 August 2017 / Revised: 13 September 2017 / Accepted: 13 September 2017 / Published: 15 September 2017

Abstract

:
Predicting system lifetime is important to ensure safe and reliable operation of products, which requires integrated modeling based on multi-level, multi-sensor information. However, lifetime characteristics of equipment in a system are different and failure mechanisms are inter-coupled, which leads to complex logical correlations and the lack of a uniform lifetime measure. Based on a Bayesian network (BN), a lifetime prediction method for systems that combine multi-level sensor information is proposed. The method considers the correlation between accidental failures and degradation failure mechanisms, and achieves system modeling and lifetime prediction under complex logic correlations. This method is applied in the lifetime prediction of a multi-level solar-powered unmanned system, and the predicted results can provide guidance for the improvement of system reliability and for the maintenance and protection of the system.

1. Introduction

Life expectancy is of great significance to ensure safe and reliable operation of products, and demands for superior product performance are increasing. Devices and components with multiple lifetime characteristics are associated by complex logic to form a multi-level, complex system whose state features are difficult to quantify and predict directly. Therefore, the integration of multi-sensor information is a necessary path to achieving lifetime prediction for such systems. In fact, research on the integration of multi-sensor information is one of the trends in comprehensive utilization of data and monitoring technology innovation in the international research community. Through signal processing from multiple sensors, combined with optimization theory [1], pattern recognition [2], filtering technique [3,4], neural networks [5], artificial intelligence [6], etc., the goal is to achieve the correlation and combination of information, and provide more accurate monitoring and evaluation of products.
In order to effectively integrate equipment information provided by sensors in a complex, multi-level system to achieve evaluation and prediction of the lifetime of a system, the following requirements must be met: first, sensor data are required to carry out evaluation at the equipment level for designing prediction models. Second, the prediction framework that can integrate these low-level prediction models must be provided. A variety of prediction methods have been developed at the equipment level, including statistical methods that are based on Wiener diffusion [7], the gamma process [8], and the Markov process [9], as well as artificial intelligence methods such as data-driven neural networks [10,11] and filtering techniques [12,13], and other methods based on physical failure [14]. For instance, Nguyen considered the predictive reliability, dependencies of components and the hierarchy configuration to study the multi-level system [15]. Zheng combined unscented Kalman Filter and Relevance Vector Regression to predict the remaining useful life of lithium-ion batteries [16]. Hamed compared the stochastic simulation and the inverse First Order Reliability Method in remaining useful life prediction under different forms of uncertainty [17]. Actually, owing to evolving sensor technology, studies on prediction models for lower level objects and simple products are prevalent, and data-based lifetime prediction technology is relatively mature. However, the construction of the system prediction framework should take into account the uniform quantification of lifetime indicators and the expression of complex logic. Coolen extended the sensor status information obtained by a sensor to the probability level, and carried out the uncertainty prediction research for many types of components [18]. Jackson and Mosleh also studied the predictive modeling of uncertainty for multilayer overlapping data in complex systems [19]. The above studies affirm the advantage of a probabilistic description of complex system states and logical association uncertainties. On the other hand, George [20] and Meineri [21] adopted directed graphs in their research on complex system topology modeling and multi-factor complex correlation, which inspired the idea of modeling logical relationships of multi-level system interactions.
The Bayesian network (BN), proposed by Pearl [22], is a reasoning model based on Bayesian theory and graph theory that integrates probability and graphitization. It can be applied to make up the deficiency of traditional reliability analysis method. For an example, with the Fault Tree Analysis (FTA) method it is difficult to model systems with polymorphism problems, random uncertainties and dependent events, and some calculations are impossible to achieve due to the computation complexity under the disjointed algorithm scenario, while BN is able to describe the multi-state, uncertainties and correlations to conduct two-way inference due to probability basis, conditional independent principle and information integration skills, as well as generate results conveniently according to the available FT and simplify modeling problems. Therefore BN has a great advantage in complex system modeling applications such as disease diagnosis, financial risk analysis, and wireless sensor network and system reliability analysis, and provides the application basis for integration of multi-level information [23], which can be used to carry out the construction of lifetime prediction models based on state probability. A BN qualitative network topology and quantitative conditional probability description are endowed with strong ability to express large complex systems that have a large number of subsystems [24], which is applicable for the modeling of complex relationships such as system reliability diagrams and fault trees that are difficult to describe using traditional models [25,26]. At present, BN is widely used in the field of reliability, for system fault diagnosis [27,28], safety analysis [29], and optimization of maintenance strategy [30]. For instance, Cai et al. proposed a BN-based data-driven fault detection and diagnosis methodology which is equipped with a good toleration for sensor noise and bias of PMSM drive system [31]. Petek conducted an evaluation of failure conditional probability with a multi-level system by BN method [32]. Hu et al. developed an integrated safety prognosis model with BN to study the propagation mechanisms of faults in complex system [33]. However, there are few reports on the prediction method for a complex, multi-level system [34]. Kabir et al. investigated the failure prediction method of water mains with BN by considering the uncertainties from multi-source data and human interpretations with different credibility [35]. Therefore, system lifetime prediction based on BN to achieve the integration of equipment information from multiple sensors has substantial research value and exploration significance.
In order to solve the problem of lifetime prediction of multi-level systems under complex logic relationships and guide the design, improvement, and maintenance of systems, this paper studies the BN lifetime prediction method based on the integration of multi-level information. The system life expectancy is estimated by considering different failure mechanisms and the logical association of its mutual coupling, and the integration of equipment information with different lifetime characteristics. This paper is organized as follows: in Section 2, we present an overview of the basic concept, construction, and estimation algorithm of the BN model. In Section 3, we describe the modeling process of BN lifetime prediction and its solution method based on the prediction model at equipment level. In Section 4, we present the verification of the prediction model and in Section 5 a case study of a solar unmanned system. Finally, a summary and outlook are given in Section 6.

2. Overview of BN Model

A BN is a graphical network model based on the probabilistic reasoning of Bayesian theory. It consists of a Directed Acyclic Graph (DAG) and a Conditional Probability Table (CPT). The former is a graphical structure composed of directed edges that connect the node variable X = { X 1 , X 2 , , X n } based on causality. The latter is the quantitative expression of the logical relationship of variables. The directed edge in a DAG always has the parent node pointing to the child node, while the variable with no parent node is the root node and the variable with no child node is the leaf node, and the rest are the intermediate nodes.
The probability calculation of a BN is based on the conditional independence assumption, i.e., that the probability of the child node depends only on the parent node and is independent from the other child nodes of the parent node, as shown in Equation (1). Therefore, when solving the probability of a multi-node joint, a BN only needs to consider the correlation of variables, thus reduces the solving complexity:
P ( X i | X p i , X p a i ) = P ( X i | X p i )
where X p i is the parent node of node X i and X p a i the child node of X p i other than X i . Applying conditional independence to chain rules enables computation of the joint probability, as follows:
P ( X 1 , X 2 , , X n ) = i = 1 n P ( X i | X p i )
The following three elements must be completed for a BN construction:
  • Determine node variables and their interpretation.
  • Create a DAG with a directed edge connecting node variables.
  • Create a CPT for non-root nodes.
The CPT of node G is established for the seven-node DAG of Figure 1 according to the series logic, as shown in Table 1.
In the precision reasoning and approximate reasoning algorithm, the joint tree (JT) algorithm based on clique tree propagation is widely used because of its unique advantages of high search efficiency, capability of returning logarithmic results, and dual-channel transmission. The solution is shown in Figure 2. The JT is obtained by the steps of the transformation and triangulation of the DAG. After initialization, the nodes can absorb the information and update the distribution function ϕ C of the partition nodes to achieve information transmission. When the single-transfer process of the JT information satisfies the globally consistent steady state, the distribution of V can be obtained according to P ( V ) = C { V } ϕ C for any desired variable V . When new evidence e is added, the conditional probability distribution of the variable V is solved as follows:
P ( V | e ) = P ( V , e ) P ( e ) = P ( V , e ) V P ( V , e ) .

3. Lifetime Prediction Method of Multi-level System Based on BN

3.1. Prediction Model at Equipment Level

From the prediction tools and means, the prediction model at the equipment and component levels can be divided into three categories: stochastic process based, data based, and physical failure based. The optimal prediction method is chosen to establish the lifetime prediction model for root nodes of the BN, and the lifetime prediction of intermediate nodes and leaf nodes can then be deduced. Therefore, the soundness and accuracy of the model and parameters of the equipment can affect the prediction of higher-level nodes.
The system of a complex product often involves performance degradation and accidental failure. Here, the equipment whose performance degradation obeys the Wiener process is taken as an example to devise a prediction model. Assuming that the performance parameter W is a key indicator of product lifetime and is sensitive to stress S , the parameter then follows the Wiener degradation process as follows:
W ( t ) = μ ( s ) t + σ B ( t ) + W 0
where W ( t ) is the product performance at time t , and μ ( s ) is the drifting coefficient reflecting the performance degradation rate, which is a function of stress and time. In an accelerated model, μ ( s ) = exp [ β 0 + β φ ( s ) ] . Constant σ is the diffusion coefficient that is irrelevant with respect to environment and time. B ( t ) ~ N ( 0 , t ) is the standard Brownian motion and W 0 is the initial value of the parameter.
The degradation amount within the time Δ t from the properties of the Wiener process is Δ W ~ N ( μ ( s ) Δ t , σ 2 Δ t ) . L is defined as the failure threshold of performance W , and then the time t that the performance parameter value first passed through L satisfies the inverse Gaussian distribution. The distribution function is the unreliability function of the product, and the corresponding probability density function is given by:
f ( t ; W 0 , L ) = L W 0 2 π σ 2 t 3 exp { [ ( L W 0 ) μ ( s ) t ] 2 2 σ 2 t }
The corresponding reliability function is the prediction model of equipment lifetime, as follows:
R ( t ) = Φ ( L W 0 ) μ ( s ) t σ t exp { 2 μ ( s ) L 2 σ 2 } Φ ( L W 0 ) + μ ( s ) t σ t ,
where Φ   ( * ) is the cumulative distribution function of the standard normal distribution.
Other prediction methods and models based on the data and physical failure can also provide the corresponding prediction information of product lifetime under the timescale satisfying prediction accuracy.

3.2. BN Prediction Modeling and Inference at System Level

Differences exist in the lifetime characteristics of different products. Therefore, the state probability of the unified lifetime based on quantitative indicators is used in this paper. For example, for the case of accidental hardware failure, reliability and cumulative probability of failure can describe the probability in the “normal” and “fault” states. However, for a degradation mechanism in which the lifetime characteristics cannot be directly described by the reliability, the complete state set must be customized before modeling to describe the “intact degree” and “failure degree” as two opposing events, to ensure that the sum of the state probability is 1. For ease of expression, in this work we use to characterize the “intact” state of all nodes, with corresponding probability R , and represent the corresponding “failure” state with , with a probability of F .
The lifetime prediction sequence of state probability reflecting lifetime information given by each equipment prediction model is taken as the prior probability of the corresponding root node in the system BN, and the probability of unknown nodes is deduced to achieve the integration and prediction of the lifetime information from the multi-sensor (process shown in Figure 3). The specific model construction proceeds via the following steps:
  • Obtain expert knowledge and structure and function information of similar products. Analyze the failure mode and mechanism through the system function-level method. Fault tree analysis is used for key faults to determine the equipment and mechanisms that affect system lifetime.
  • Deploy sensors for key performance parameters of each device. Collect and process data.
  • Based on the analysis of sensor data, the prediction model is established for each piece of equipment, and the prediction value of lifetime-related state probability is given.
  • Combined with the system failure mechanism, the prediction value of state probability for the equipment involved is used as the prior probability for the root node.
  • Combine system logic to form a DAG, and establish a non-root node CPT.
  • The JT estimation algorithm is used to solve the joint probability of relevant nodes, to update the conditional probability values of each node, and to achieve the deduction of state probability of system nodes to complete the system prediction.
  • The BN prediction model is still applicable, along with the updating of sensor data and correction of prediction model at equipment level. If the failure mechanism changes, the DAG and CPT should be corrected for the updated logical relationship; proceed to Step 3.
In the above BN model that interprets node variable with state probability, S X i a i ( t ) ,   S Y j b j ( t ) , and S L k ( t ) are used to represent the a i ,   b j ,   k   th state of root node X i   ( i = 1 , 2 , , p ) , the intermediate node Y j   ( j = 1 , 2 , , q ) , and the leaf node L at time t respectively, i.e., a i = b j = k = 1 , 2 .
For a number of p root nodes, the probability of the state is solved by the prediction model for each respective piece of equipment and discretized according to the unsupervised equal-width interval method with a time sequence to achieve the state prediction in the future T time; that is, the p × n -order state probability prediction matrix:
R p × n = ( R 1 , 1 R 1 , n R p , 1 R p , n ) ,
with a certain time sequence T = ( t + ζ , t + 2 ζ , , t + n ζ ) . The elements R i , τ = R ( X i ( τ ) )   ( i = 1 , 2 , , p ;   τ = 1 , 2 , , n ) represent the probability of the device X i to be in the state R . If ζ is taken as the unit time, the probability set R ( t ) = ( R 1 , t , R 2 , t , , R p , t ) T corresponds to the probability of p nodes to be in set X at state . Correspondingly, the probability of occurrence of state F is given by F ( t ) = 1 R ( t ) .
When the abovementioned multi-sensor information is used to deduce the lifetime prediction based on the BN, the probability of the root node is firstly assigned according to the probability prediction matrix. The probability of the root node X i at time t is P ( S X i ( t ) ) = R ( X i ( t ) ) , i = 1 , 2 , , p . For the solution of the state probability of intermediate nodes, it is assumed that the parent-node set X = { X 1 , X 2 , , X i } exists for the node Y j . According to the assumption of independent conditions, the probability prediction of the intermediate nodes at time t can be solved based on:
P ( Y j ( t ) ) = X P ( S Y j ( t ) , S X ( t ) ) = X P ( S Y j 2 ( t ) | S X ( t ) ) P ( S X 1 ( t ) ) P ( S X 1 ( t ) )   ,
The complexity of logical associations of the BN will increase with the number of nodes contained. Based on the advantages of conditional independence, the complex topological equivalence can be divided into simple structures, and the construction and inference of three kinds of basic structural prediction models are taken as examples as follows:
(1) The parent node set X = { X 1 , X 2 , X 3 } of node Y j contains only the root node (Figure 4).
Each state of the parent node is independent of each other at any time, and then the predicted probability of Y j in state at time t is as follows:
R ( Y j ( t ) ) = X 1 , X 2 , X 3 P ( S Y j ( t ) , S X 1 ( t ) , S X 2 ( t ) , S X 3 ( t ) ) = X 1 , X 2 , X 3 P ( S Y j 2 ( t ) | S X ( t ) ) P ( S X 1 ( t ) ) P ( S X 2 ( t ) ) P ( S X 3 ( t ) ) = P ( S Y j 2 ( t ) | S X 1 2 ( t ) , S X 2 2 ( t ) , S X 3 2 ( t ) ) P ( S X 1 2 ( t ) ) P ( S X 2 2 ( t ) ) P ( S X 3 2 ( t ) )   .
(2) The parent node set X = { X a , Y 1 } of node Y j contains both the root node and the intermediate node (Figure 5).
The predicted probability of node Y j at time t can be solved according to:
R ( Y j ( t ) ) = Y 1 , X a P ( S Y j 2 ( t ) | S Y 1 2 ( t ) , S X a 2 ( t ) ) P ( S X a 2 ( t ) ) X 1 , X 2 P ( S Y 1 2 ( t ) | S X 1 ( t ) , S X 2 ( t ) ) P ( S X 1 ( t ) ) P ( S X 2 ( t ) )   = P ( S Y j 2 ( t ) | S Y 1 2 ( t ) , S X a 2 ( t ) ) P ( S X a 2 ( t ) ) · [ P ( S Y 1 2 ( t ) | S X 1 2 ( t ) , S X 2 2 ( t ) ) P ( S X 1 2 ( t ) ) P ( S X 2 2 ( t ) ) + P ( S Y 1 2 ( t ) | S X 1 1 ( t ) , S X 2 2 ( t ) ) P ( S X 1 1 ( t ) ) P ( S X 2 2 ( t ) ) + P ( S Y 1 2 ( t ) | S X 1 2 ( t ) , S X 2 1 ( t ) ) P ( S X 1 2 ( t ) ) P ( S X 2 1 ( t ) ) ]   .
(3) There is a case in which the same parent node points to multiple child nodes (see Figure 6).
The system model that has interconnected logical correlation mostly has network structure. The same parent node pointing to multiple sub-nodes is the basis for the composition of the network structure. This situation can be predicted according to:
R ( Y j ( t ) ) = X 2 , Y 1 P ( S Y j 2 ( t ) | S X 2 ( t ) , S Y 1 ( t ) ) X 1 , X 2 P ( S Y 1 ( t ) | S X 1 ( t ) , S X 2 ( t ) ) P ( S X 1 ( t ) ) P ( S X 2 ( t ) ) = P ( S Y j 2 ( t ) | S X 2 2 ( t ) , S Y 1 2 ( t ) ) [ P ( S Y 1 2 ( t ) | S X 1 2 ( t ) , S X 2 2 ( t ) ) P ( S X 1 2 ( t ) ) P ( S X 2 2 ( t ) ) + P ( S Y 1 2 ( t ) | S X 1 1 ( t ) , S X 2 2 ( t ) ) P ( S X 1 1 ( t ) ) P ( S X 2 2 ( t ) ) ]   .
Based on the probability of the root node and intermediate node, the predicted probability of the leaf node in state R can be further solved according to:
R L ( t ) = P ( S X 1 ( t ) , S X 2 ( t ) , , S X p ( t ) , S Y 1 ( t ) , S Y 2 ( t ) , , S Y q ( t ) , S L 2 ( t ) ) = P a ( L ) P ( S L 2 ( t ) | S P a ( L ) ( t ) ) P a ( Y 1 ) P ( S Y 1 ( t ) | S P a ( Y 1 ) ( t ) ) P a ( Y q ) P ( S Y q ( t ) | S P a ( Y q ) ( t ) ) P ( S X 1 ( t ) ) P ( S X p ( t ) )   ,
where P a ( ) is the parent node of the node “ ( ) ”.
According to the above formula, the JT estimation algorithm traverses the DAG, and the state and lifetime of the system node L can be predicted. By the probability prediction matrix R p × n = ( R ( t ) , R ( t + ζ ) , R ( t + 2 ζ ) , , R ( t + 2 ζ ) ) of the root node, the corresponding prediction sequence of probability at system level will be obtained to achieve continuous prediction of the lifetime.

4. Modeling, Simulation, and Verification of System-Level BN Prediction

The seven-node DAG shown in Figure 1 contains the three basic structures described above, and the BN prediction method is simulated and verified as an example. The three-tier system G consists of parallel subsystems C and F , comprising four devices A , B , D , and E . In the subsystem C ( F ) , A   ( D ) and B   ( E ) conduct different functions, respectively, which are in series logic. In addition, according to the early collection of product information, the logic that “if the state of subsystem C is abnormal, there is a 60% probability of failure for subsystem F ” exists between the subsystems. The key performance parameters of each piece of equipment were determined, monitoring devices were arranged, and the sensor signals were extracted, processed, and analyzed. The lifetime prediction model of each piece of equipment was then obtained.
(1) The lifetime characteristics of equipment A and B meet the accidental failure, while D and E meet the degradation failure, and the corresponding prediction models are shown in Table 2.
Based on the prediction model above, the BN model with both accidental failure and degeneration failure is established to predict the system lifetime of the next 2000 h. The state probability of system G can be solved by:
R G ( t ) = P ( S A ( t ) , S B ( t ) , S C ( t ) , S D ( t ) , S E ( t ) , S F ( t ) , S G 2 ( t ) ) = C , F P ( S G 2 ( t ) | S C ( t ) , S F ( t ) ) A , B P ( S C ( t ) | S A ( t ) , S B ( t ) )   C , D , E P ( S F ( t ) | S D ( t ) , S E ( t ) , S C ( t ) ) P ( S A ( t ) ) P ( S B ( t ) ) P ( S D ( t ) ) P ( S E ( t ) ) .
Due to the efficient two-way accurate-inference ability, “engine = jtree _ inf _ engine (bnet)” of the Bayes Net Toolbox (BNT) for MATLAB developed by Murphy [36] is used to conduct inference, and the probability of each group of parent and child nodes in state R is shown in Figure 7.
The results show that in the first 1000 h the subsystem F shows better performance than C , and the state probability of system G is the same as that of the former. After 1030 h, however, with the performance degradation, the failure probability of F suddenly increases, and the probability of system G being in state R decreases. After 1090 h, subsystem F completely failed, and the system state was completely determined by subsystem C . Therefore, the prediction curve trends of C and G after this point are basically the same. The lifetime prediction of the system is estimated and evaluated: the median life of the system t 0.5 = 1098   h ; if the system in good condition has a probability threshold of 0.45, the remaining life is approximately 1226 h.
In addition, the performance degradation simulation curve of E (as shown in Figure 8) shows that from 1010 to 1090 h the performance parameters of all 15 samples degenerated below the threshold of 80, and that the 15 samples of device D did not exceed the performance parameter threshold in the first 1200 h (a simulated degradation curve is shown in Figure 9). Through the analysis of the subsystem F , we can see that the main reason for the sudden deterioration of the system is the performance degradation of device E under cyclic stress, which caused it to exceed the performance threshold with considerable probability.
The dotted lines parallel to the horizontal axis of Figure 8 and Figure 9 represent the device performance threshold.
(2) For all equipment that is in degradation failure mode, the corresponding prediction model and parameters are shown in Table 3.
The state probability prediction curve of each (sub-) system is obtained and shown in Figure 10.
Compared to the predicted result above, a system whose devices are all in degradation failure mode has a low probability of staying intact at 1030 h, specifically, only 0.622. However, the state probability predicted above at 1030 h is 0.712, and the working state was maintained with a probability of 0.46 until 1200 h. With a probability threshold of 0.45, the remaining life of the system will be less than 1064 h. The main reason for the change of system state is subsystem C , in which the performance parameter of device A follows the Gamma degradation process (Figure 11). The probability of device A to be at state 2 is reduced to 0.5224 at 1000 h, making the system state significantly degraded; after this prediction point, the performance degradation of device E in subsystem F is deteriorated, and the system fails at a faster rate. On one hand, subsystem C has direct effect on the state of system G . On the other hand, the state of G is indirectly affected through the correlation with subsystem F . Therefore, the degradation rate of a system state after integration is faster.
It can be verified through the above simulation that a BN can achieve the integration of lifetime information from multiple sensors in a multi-level system when the equipment grade product experiences accidental failure or (and) degradation failure, which reasonably predicts system lifetime from the perspective of state probability, and shows good compatibility and integration for a variety of low-level prediction models.

5. Application Case

5.1. BN Lifetime Prediction of an Unmanned Solar-Powered System

An unmanned system is a typical multi-level system with complex structure and function logic. With the air-ground collaborative work requirements, the system achieves its functional requirements and sharing of information resources layer by layer. At the same time, under environmental influences such as temperature, humidity, vibration, and interference the unit with the degradation property and that with the accidental failure together affect the system state and lifetime. In this paper, a solar-powered unmanned aerial vehicle (UAV) system was chosen as an example to carry out lifetime prediction based on a BN.
After the system is simplified, the hardware is mainly composed of an airborne system, a data link, and a ground system. The energy system inside the airborne system must provide power support for the normal operation of the navigation system, flight control system, power system, communication system, etc., and it is subject to the state of the battery management system and the solar panel. Any abnormality can cause the energy system to degrade or even fail to function properly.
The navigation system contains the master and lead. A flight control computer subsystem for integrated task management, a flight control computer subsystem, and a servo action system are necessary. In addition to the airborne terminal, the airborne communication system must go through the data link to achieve the transmission and exchange of information with ground communication equipment. In the ground subsystem, the power supply system is the energy supply system for flight operation, the console, communication equipment, and other equipment. For the console, the requirement is to simultaneously display and control functions. This BN model fully accounts for the relevance of intermediate nodes. In these subsystems, degradation failure exists in the solar cell and UAV body material, which is modeled with Wiener process and Gamma process, respectively, to establish the performance degradation for lifetime prediction. Specifically, the system DAG is shown in Figure 12, and the 61 nodes contain 37 device-level nodes.
At the same time, considering the characteristics of photovoltaic energy, a special flight profile is used: To ensure the full storage of the solar panels, before every sunrise the aircraft must climb from an altitude of approximately 1000 m to a cruising altitude to 8500 m. In order to reduce energy consumption, a fully charged UAV descends to a low-altitude region in the evening to take advantage of lower air resistance and will continue to fly. The ideal charge-discharge process under the flight height profile is shown in Figure 13.
During the process, the ambient temperature will change cyclically with flight height, and the structure with the temperature degradation characteristic will experience significant state change. A simplified temperature stress profile is given in Figure 14.
The lifetime prediction model considering the temperature stress is established by combining the sensor information of the 37 devices, and the state of the cruise mission in the next 35 d is predicted according to the above temperature profile. Figure 15 shows the state probability distribution of the solar-powered unmanned system and its data link, ground subsystem, and airborne subsystem.
The continual decreasing trend reflects the degradation of operation state of the (sub-) system with time and environmental conditions; the UAV system and airborne subsystem show cyclically violent state changes after flying 100 h approximately. This is because the temperature degradation characteristics of the airborne equipment are subject to periodic ambient temperature stress.

5.2. Analysis Based on the Prediction Results

(1) Task risk analysis
In order to make full use of the predicted information to identify the health of the UAV system, a task risk analysis was carried out on the system based on health status. Based on fuzzy theory, the probabilistic predictive values for each time are divided into four healthy states, such as “healthy”, “minor fault”, “catastrophic fault”, and “close to failure” with the trapezoidal fuzzy number. The probability distribution curves of four states are shown as dotted lines in Figure 16. The risk of each state can be weighted to achieve the risk assessment of system at all times (as shown in Figure 16).
It can be seen from Figure 16 that the UAV system can work normally in the first 438 h, but frequent abnormal jumps between “healthy” and “close to failure” occurred between 439 and 496h, and the worst health state was maintained after 497 h, which leads to a continuous increase in systemic risk, especially after 512 h when the task risk has exceeded 0.5. Continuing to perform the task poses a considerable probability of failure and safety risks for personnel and equipment. Additionally, at the time of 439 h when state mutation firstly occurs, the reliability of each device is predicted and shown in Figure 17.
It can be seen from the figure that the probability that the lithium battery and the data link (nodes 16, 20, 24, 28, and 37) can work normally is relatively low when the system state begins to appear abnormally, which may be the fundamental factor causing the deterioration of the system state.
(2) Sensitivity analysis of nodes
Sensitivity can be used to reflect the sensitivity of system output to the amount of input change so that the weak parts of the system can be identified. Due to the conditional independence assumption in BN, nodes at the same level are independent of each other, and thus the sensitivity of each node to the system is defined as:
Θ i , S = R s R X i = Δ R s Δ R X i
In the BN model of the UAV system, the sensitivity of each subsystem is shown in Figure 18.
According to the figure, the system is more sensitive to the state change of node 60 (airborne subsystem). Further analysis of the airborne subsystem shows that the sensitivity of nodes 16, 20, 24, and 28 can reach 14 orders of magnitude at a certain local time (the figures of which are not shown in this paper). In addition, node 37 (data link) is both a subsystem and root node whose sensitivity is at a high level. As a result, the airborne subsystem and data link subsystem will have a significant impact on the system state.
(3) Recommendations for reliability design and maintenance
In reliability design, sensitivity analysis can be used to improve the inherent reliability of critical equipment by design improvement. Meanwhile, close attention should be paid to redundancy design within the allowable range of aircraft load. This will create significant improvements in the reliability, availability, and longevity of UAV systems. From the point of view of maintenance and repair, the maintenance strategy T ( C , A * , P r i s k | t ) related to maintenance cost C ( t ) , instantaneous availability A * ( s ) and work risk P r i s k ( t ) can be formulated according to the given lifetime index. Combined with system mission risk and the average cost of preventive and post-maintenance, a maintenance cost-time model can be established. Maintenance or replacement activities are carried out at higher cost-efficiency-ratio point in accordance with the priority of the equipment status and sensitivity.

6. Conclusions

In this paper we propose a modeling method for lifetime prediction of systems from integrated multi-sensor information at device level based on a BN in order to conduct system-level lifetime prediction for complex, multi-sensor systems. The simulation results show that the method can effectively solve the complex logical association of data and the quantitative description of the uncertainty in the multi-level system, and that it has good integration ability and compatibility with various prediction models of devices. At the same time, the application research of a multi-level solar-powered UAV system with multi-device information is carried out, and the task risk and node sensitivity analysis are conducted according to the prediction results of system lifetime based on the temperature profile. The research shows that the system prediction method from the integrated multi-sensor lifetime based on a BN has the following advantages:
(1)
By making full use of multi-sensor information, data association, and quantitative expression, integration and lifetime prediction can be achieved on multi-level, complex, dynamic, multi-source logic data.
(2)
The diversified life expectancy information of the model output is based on different aspects and complement each other, which provides comprehensive data support and theoretical guidance for design improvements and maintenance from system view and in the entire lifecycle process that considers a trade-off of economics, technology, risk, and effectiveness.
However, simulation and application studies show that the proposed method has a high dependence on the prediction model at the device level. The more accurate lifetime prediction for a complex multi-level system is based on accurate device-level prediction information, which imposes higher requirements on information acquisition, processing, and analysis from multiple sensors. Therefore, there might be limitations in practice by only considering data from sensors. Future work on lifetime prediction methods based on a BN should focus on the study and understanding of prediction accuracy, and appropriately take into account other sources and levels data like experimental data and reasonable expert knowledge and so on to conduct integration. In addition, the prediction dimension can be extended from two states to multiple states, and the intermediate process of various mechanisms can be studied in depth by considering the intermediate state to obtain more abundant lifetime prediction information.

Acknowledgments

This work was supported by the Aero-Science Fund (Grant No. 2015ZD51044).

Author Contributions

Jingbin Wang and Lizhi Wang proposed the idea of the research, conceived and designed the simulation experiments; Xiaohong Wang and Lizhi Wang conceived and provided the application case; Jingbin Wang performed the experiments, analyzed the data and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ge, Q.; Wei, Z.; Cheng, T.; Chen, S.; Wang, X. Flexible Fusion Structure-Based Performance Optimization Learning for Multisensor Target Tracking. Sensors 2017, 17, 1045. [Google Scholar] [CrossRef] [PubMed]
  2. Xu, H.; Liu, J.; Hu, H.; Zhang, Y. Wearable Sensor-Based Human Activity Recognition Method with Multi-Features Extracted from Hilbert-Huang Transform. Sensors 2016, 16, 2048. [Google Scholar] [CrossRef] [PubMed]
  3. Li, Z.; Feng, L.; Yang, A. Fusion Based on Visible Light Positioning and Inertial Navigation Using Extended Kalman Filters. Sensors 2017, 17, 1093. [Google Scholar]
  4. Rodger, J.A. Toward reducing failure risk in an integrated vehicle health maintenance system: A fuzzy multi-sensor data fusion Kalman filter approach for IVHMS. Expert Syst. Appl. 2012, 39, 9821–9836. [Google Scholar] [CrossRef]
  5. Jing, L.; Wang, T.; Zhao, M.; Wang, P. An Adaptive Multi-Sensor Data Fusion Method Based on Deep Convolutional Neural Networks for Fault Diagnosis of Planetary Gearbox. Sensors 2017, 17, 414. [Google Scholar] [CrossRef] [PubMed]
  6. Hsu, Y.L.; Chou, P.H.; Chang, H.C.; Lin, S.L.; Yang, S.C.; Su, H.Y.; Chang, C.C.; Cheng, Y.S.; Kuo, Y.C. Design and Implementation of a Smart Home System Using Multisensor Data Fusion Technology. Sensors 2017, 17, 1631. [Google Scholar] [CrossRef] [PubMed]
  7. Whitmore, G.A.; Schenkelberg, F. Modelling accelerated degradation data using Wiener diffusion with a time scale transformation. Lifetime Data Anal. 1997, 3, 27–45. [Google Scholar] [CrossRef] [PubMed]
  8. Bagdonavicius, V.; Nikulin, M.S. Estimation in degradation models with explanatory variables. Lifetime Data Anal. 2001, 7, 85–103. [Google Scholar] [CrossRef] [PubMed]
  9. Wan, Y.; Huang, H.L.; Das, D.; Pecht, M. Thermal reliability prediction and analysis for high-density electronic systems based on the Markov process. Microelectron. Reliab. 2016, 56, 182–188. [Google Scholar] [CrossRef]
  10. Ali, J.B.; Chebel-Morello, B.; Saidi, L.; Malinowski, S.; Fnaiech, F. Accurate Bearing remaining useful life prediction based on Weibull distribution and artificial neural network. Mech. Syst. Sig. Process. 2015, 56, 150–172. [Google Scholar]
  11. Alghassi, A.; Perinpanayagam, S.; Samie, M. Stochastic RUL Calculation Enhanced With TDNN-Based IGBT Failure Modeling. IEEE Trans. Reliab. 2016, 65, 558–573. [Google Scholar] [CrossRef]
  12. Son, J.; Zhou, S.; Sankavaram, C.; Sankavaram, C.; Du, X.; Zhang, Y. Remaining useful life prediction based on noisy condition monitoring signals using constrained Kalman filter. Reliab. Eng. Syst. Saf. 2016, 152, 38–50. [Google Scholar] [CrossRef]
  13. Li, K.; Wu, J.; Zhang, Q.; Su, L.; Chen, P. New Particle Filter Based on GA for Equipment Remaining Useful Life Prediction. Sensors 2017, 17, 696. [Google Scholar] [CrossRef] [PubMed]
  14. Qian, Y.N.; Yan, R.Q.; Gao, R.X. A multi-time scale approach to remaining useful life prediction in rolling bearing. Mech. Syst. Sig. Process. 2017, 83, 549–567. [Google Scholar] [CrossRef]
  15. Nguyen, K.A.; Do, P.; Grall, A. Multi-level predictive maintenance for multi-component systems. Reliab. Eng. Syst. Saf. 2015, 144, 83–94. [Google Scholar] [CrossRef]
  16. Zheng, X.; Fang, H. An integrated unscented kalman filter and relevance vector regression approach for lithium-ion battery remaining useful life and short-term capacity prediction. Reliab. Eng. Syst. Saf. 2015, 144, 74–82. [Google Scholar] [CrossRef]
  17. Khorasgani, H.; Biswas, G.; Sankararaman, S. Methodologies for system-level remaining useful life prediction. Reliab. Eng. Syst. Saf. 2016, 154, 8–18. [Google Scholar] [CrossRef]
  18. Coolen, F.P.A.; Coolen-Maturi, T. The structure function for system reliability as predictive (imprecise) probability. Reliab. Eng. Syst. Saf. 2016, 154, 180–187. [Google Scholar] [CrossRef] [Green Version]
  19. Jackson, C.; Mosleh, A. Bayesian inference with overlapping data: Reliability estimation of multi-state on-demand continuous life metric systems with uncertain evidence. Reliab. Eng. Syst. Saf. 2016, 145, 124–135. [Google Scholar] [CrossRef]
  20. George-Williams, H.; Patelli, E. A hybrid load flow and event driven simulation approach to multi-state system reliability evaluation. Reliab. Eng. Syst. Saf. 2016, 152, 351–367. [Google Scholar] [CrossRef]
  21. Meineri, E.; Dahlberg, C.J.; Hylander, K. Using Gaussian Bayesian Networks to disentangle direct and indirect associations between landscape physiography, environmental variables and species distribution. Ecol. Modell. 2015, 313, 127–136. [Google Scholar] [CrossRef]
  22. Pearl, J. Fusion, propagation, and structuring in belief networks. Artif. Intell. 1986, 29, 241–288. [Google Scholar] [CrossRef]
  23. Weber, P.; Medina-Oliva, G.; Simon, C.; Iung, B. Overview on Bayesian networks applications for dependability, risk analysis and maintenance areas. Eng. Appl. Artif. Intell. 2012, 25, 671–682. [Google Scholar] [CrossRef]
  24. Tien, I.; Kiureghian, A.D. Algorithms for Bayesian network modeling and reliability assessment of infrastructure systems. Reliab. Eng. Syst. Saf. 2016, 156, 134–147. [Google Scholar] [CrossRef]
  25. Price, L.C.; Peiris, H.V.; Frazer, J.; Easther, R. Designing and testing inflationary models with Bayesian networks. J. Cosmol. Astropart. Phys. 2016, 2016, 049. [Google Scholar] [CrossRef]
  26. Garvey, M.D.; Carnovale, S.; Yeniyurt, S. An analytical framework for supply network risk propagation: A Bayesian network approach. Eur. J. Oper. Res. 2015, 243, 618–627. [Google Scholar] [CrossRef]
  27. Cai, B.; Liu, H.; Xie, M. A real-time fault diagnosis methodology of complex systems using object-oriented Bayesian networks. Mech. Syst. Sig. Process. 2016, 80, 31–44. [Google Scholar] [CrossRef]
  28. Cai, B.; Liu, Y.; Fan, Q.; Zhang, Y.; Liu, Z.; Yu, S.; Ji, R. Multi-source information fusion based fault diagnosis of ground-source heat pump using Bayesian network. Appl. Energy 2014, 114, 1–9. [Google Scholar] [CrossRef]
  29. Khakzad, N. Application of dynamic Bayesian network to risk analysis of domino effects in chemical infrastructures. Reliab. Eng. Syst. Saf. 2015, 138, 263–272. [Google Scholar] [CrossRef]
  30. Tang, D.; Makis, V.; Jafari, L.; Yu, J. Optimal maintenance policy and residual life estimation for a slowly degrading system subject to condition monitoring. Reliab. Eng. Syst. Saf. 2015, 134, 198–207. [Google Scholar] [CrossRef]
  31. Cai, B.; Zhao, Y.; Liu, H.; Xie, M. A Data-Driven Fault Diagnosis Methodology in Three-Phase Inverters for PMSM Drive Systems. IEEE Trans. Power Electron. 2017, 32, 5590–5600. [Google Scholar] [CrossRef]
  32. Yontay, P.; Pan, R. A computational Bayesian approach to dependency assessment in system reliability. Reliab. Eng. Syst. Saf. 2016, 152, 104–114. [Google Scholar] [CrossRef]
  33. Hu, J.; Zhang, L.; Ma, L.; Liang, W. An integrated safety prognosis model for complex system based on dynamic Bayesian network and ant colony algorithm. Expert Syst. Appl. 2011, 38, 1431–1446. [Google Scholar] [CrossRef]
  34. Li, N.; Feng, X.; Jimenez, R. Predicting rock burst hazard with incomplete data using Bayesian networks. Tunn. Undergr. Space Technol. 2017, 61, 61–70. [Google Scholar] [CrossRef]
  35. Kabir, G.; Demissie, G.; Sadiq, R.; Tesfamariam, S. Integrating failure prediction models for water mains: Bayesian belief network based data fusion. Knowl.-Based Syst. 2015, 85, 159–169. [Google Scholar] [CrossRef]
  36. Murphy, K. Bayesian Network Toolbox (BNT). Available online: http://www.cs.ubc.ca/murphyk/Software/BNT/bnt.html (accessed on 11 August 2017).
Figure 1. Seven-node DAG of system G .
Figure 1. Seven-node DAG of system G .
Sensors 17 02123 g001
Figure 2. Algorithm and build thought of JT.
Figure 2. Algorithm and build thought of JT.
Sensors 17 02123 g002
Figure 3. Modeling process of BN prediction at system level based on integrated multi-sensor information.
Figure 3. Modeling process of BN prediction at system level based on integrated multi-sensor information.
Sensors 17 02123 g003
Figure 4. Logical structure in which the parent node is the root node.
Figure 4. Logical structure in which the parent node is the root node.
Sensors 17 02123 g004
Figure 5. Logical structure in which the parent node contains both the intermediate and root nodes.
Figure 5. Logical structure in which the parent node contains both the intermediate and root nodes.
Sensors 17 02123 g005
Figure 6. Logical structure in which the same parent node points to multiple subnodes.
Figure 6. Logical structure in which the same parent node points to multiple subnodes.
Sensors 17 02123 g006
Figure 7. Probability distribution (1) of parent and child nodes at state R for G system group.
Figure 7. Probability distribution (1) of parent and child nodes at state R for G system group.
Sensors 17 02123 g007
Figure 8. Simulation curve of performance degradation of device E . .
Figure 8. Simulation curve of performance degradation of device E . .
Sensors 17 02123 g008
Figure 9. Simulation curve of performance degradation of device D .
Figure 9. Simulation curve of performance degradation of device D .
Sensors 17 02123 g009
Figure 10. Probability distribution (2) of parent and child nodes at state R for G system group.
Figure 10. Probability distribution (2) of parent and child nodes at state R for G system group.
Sensors 17 02123 g010
Figure 11. Simulation curve of performance degradation of equipment A .
Figure 11. Simulation curve of performance degradation of equipment A .
Sensors 17 02123 g011
Figure 12. DAG of a solar-powered unmanned system.
Figure 12. DAG of a solar-powered unmanned system.
Sensors 17 02123 g012
Figure 13. Correspondence between flight height profile of UAV and charge-discharge cycle.
Figure 13. Correspondence between flight height profile of UAV and charge-discharge cycle.
Sensors 17 02123 g013
Figure 14. Simplified temperature profile of solar-powered UAV flight.
Figure 14. Simplified temperature profile of solar-powered UAV flight.
Sensors 17 02123 g014
Figure 15. Probability prediction curve of unmanned system state.
Figure 15. Probability prediction curve of unmanned system state.
Sensors 17 02123 g015
Figure 16. Probability distribution and risk prediction of each state of unmanned system.
Figure 16. Probability distribution and risk prediction of each state of unmanned system.
Sensors 17 02123 g016
Figure 17. Reliability prediction of each piece of equipment at 439 h.
Figure 17. Reliability prediction of each piece of equipment at 439 h.
Sensors 17 02123 g017
Figure 18. System sensitivity to three subsystems.
Figure 18. System sensitivity to three subsystems.
Sensors 17 02123 g018
Table 1. CPT of node G .
Table 1. CPT of node G .
C F G
FTFTFT
101010
100101
011001
010101
Table 2. Prediction model (1) of four pieces of equipment.
Table 2. Prediction model (1) of four pieces of equipment.
NodeDescription of Prediction ModelLifetime Prediction Model
A Exponential distribution R A ( t ) = exp ( t / 3000 )
B Weibull distribution   R B ( t ) = exp [ ( t / 2300 ) ] 1.5
D Wiener degradation process; Drift parameters are selected as the Arrhenius model R D ( t ) = exp [ 2 d ( s ) ( L Y 0 ) σ 2 ] Φ ( L Y 0 + d ( s ) t σ t ) + Φ ( L Y 0 d ( s ) t σ t ) ( E a D = 0.473 , A D = 4.7 × 10 5 ; Y 0 D = 5 , L D = 10 , σ D = 0.01 )
E Degradation under cyclic stress R E ( t ) = Φ ( ( 0.0001 ( t / 10 ) 3 + 200 80 ) / 5 )
Table 3. Prediction model (2) of four pieces of equipment 1.
Table 3. Prediction model (2) of four pieces of equipment 1.
NodeDescription of Prediction ModelLifetime Prediction Model
A Gamma process; Scale parameters are selected as the Arrhenius model R A ( t ) = Γ ( v ( s ) t , L / u ) / Γ ( v ( s ) t ) ( E a A = 0.453 , A A = 3.9 × 10 5 ; Y 0 A = 0 , L A = 12 , u A = 2.3 )
B Wiener process ( E a B = 0.473 , A B = 6.7 × 10 5 ; Y 0 B = 6 , L B = 12 , σ B = 0.032 )
1 Prediction models of D and E are the same as in Table 3.

Share and Cite

MDPI and ACS Style

Wang, J.; Wang, X.; Wang, L. Modeling of BN Lifetime Prediction of a System Based on Integrated Multi-Level Information. Sensors 2017, 17, 2123. https://doi.org/10.3390/s17092123

AMA Style

Wang J, Wang X, Wang L. Modeling of BN Lifetime Prediction of a System Based on Integrated Multi-Level Information. Sensors. 2017; 17(9):2123. https://doi.org/10.3390/s17092123

Chicago/Turabian Style

Wang, Jingbin, Xiaohong Wang, and Lizhi Wang. 2017. "Modeling of BN Lifetime Prediction of a System Based on Integrated Multi-Level Information" Sensors 17, no. 9: 2123. https://doi.org/10.3390/s17092123

APA Style

Wang, J., Wang, X., & Wang, L. (2017). Modeling of BN Lifetime Prediction of a System Based on Integrated Multi-Level Information. Sensors, 17(9), 2123. https://doi.org/10.3390/s17092123

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