Research Article  Open Access
A KrigingBased Active Learning Algorithm for Mechanical Reliability Analysis with TimeConsuming and Nonlinear Response
Abstract
When the reliability analysis of the mechanical products with high nonlinearity and timeconsuming response is carried out, there will be the problems of low precision and huge computation using the traditional reliability methods. To solve these issues, the active learning reliability methods have been paid much attention in recent years. It is the key to choose an efficient learning function (such as U, EFF, and ERF). The aim of this study is to further decrease the computation and improve the accuracy of the reliability analysis. Inspired from these learning functions, a new pointselected learning function (called HPF) is proposed to update DOE, and a new point is sequentially added step by step to the DOE. The proposed learning function can consider the features like the sampling density, the probability to be wrongly predicted, and the local and global uncertainty close to the limit state. Based on the stochastic property of the Kriging model, the analytic expression of HPF is deduced by averaging a hybrid indicator throughout the real space. The efficiency of the proposed method is validated by two explicit examples. Finally, the proposed method is applied to the mechanical reliability analysis (involving timeconsuming and nonlinear response). By comparing with traditional mechanical reliability methods, the results show that the proposed method can solve the problems of large computation and low precision.
1. Introduction
In the mechanical reliability analysis (MRA), when the joint probability density function f_{X}(x) of the random variable vector X = [X_{1}, X_{2}, …, X_{n}]^{T} is known, the failure probability can be expressed aswhere is the performance function of the mechanical system. If > 0, the system is in the safety state; otherwise, it is in the failure state. and the failure domain F and the safety domain S of input random variables are and , respectively.
However, the performance function of the MRA is obtained by solving the numerical model (such as finite element method and finite different method) which is usually timeconsuming and highly nonlinear. For this situation, the moment reliability methods [1] are not applicable because the calculated reliability’s accuracy is very low due to ignorance of the Taylor expansion’s highorder items. While the random simulation methods [2, 3] (such as Monte Carlo simulation (MCS), important sampling (IS), subset simulation (SS), line sampling (LS), and Latin hypercube sampling) are unaffordable because their calculation computation is too huge for calling a large of numerical models. It also takes too much time to calculate the reliability by the workstation computer. Therefore, in MRA, it is the core problem to minimize the number of calls to the numerical calculation under the premise of ensuring the calculated accuracy.
To solve the problem of large computation in MRA, several surrogate models (including polynomial response surfaces [1], sparse polynomials [4], neural networks [5], support vector machines [6], and Kriging [7]) are usually replaced with the real performance function. Since the speed of calculating the surrogate function is much faster than the real performance function, the surrogate model techniques are widely applied to MRA. However, there are two weaknesses in these surrogate model techniques: firstly, the sample points of the design of experiment (DOE) must be arranged in advance before the surrogate model of the performance function is constructed. This fixed DOE’s scheme will cause excessive sampling of areas where design variables are not important and then results in unnecessary computational waste. Secondly, the traditional surrogatebased reliability methods only guarantee that the fitting performance function is as accurate as possible, while the accuracy of the reliability analysis cannot be guaranteed.
Kriging, as a surrogate model, shows better fitting accuracy than other surrogate models [4, 8, 9]. Kriging not only is a interpolation model but also provides the mean and variance at an unknown point. Because of the statistical information of Kriging, the Krigingbased active learning reliability method (KARM) has received much attention in the recent years. The basic idea of KARM is to sequentially refresh the DOE, and the constructed surrogate model gradually reaches to the real performance function. Inspired by the famous efficient global optimization (EGO) method [10], Bichon et al. [11, 12] proposed an efficient global reliability analysis (EGRA) method to search the DOE’s points in the vicinity of the limit state. Echard et al. proposed an active learning reliability method (AKMCS/AKIS) [13, 14] by combing MCS/IS and Krigng. Yang et al. [15] proposed a reliability mixed algorithm (ALKHRA) based on probability and convex set hybrid reliability analysis. Wen et al. [16] proposed a dynamic sampling area strategy (SKRA) to avoid selecting samples in areas with low probability density, to improve the efficiency of selecting the best point. Among these methods, AKMCS [13] is one of the most typical KARMs, whose learning function U is an innovative pointselected learning function that measures the probability of classifying a point into the wrong domain. Based on AKMCS, there are many scholars who have extended the study of U in order to solve some special conditions (such as small failure probability, multiple dimension, multiple failure regions, and systems) of reliability analysis. For example, Echard et al. [14] combined the AKMCS algorithm idea with IS (AKIS) to make it suitable for solving the reliability problem with small failure probability. Tong et al. [17] combined the AKMCS idea with the SSIS method for the small failure probability problem, to improve the robustness of multidimensional random variables. Cadini et al. [18] improved the AKIS to solve the reliability problem of multiple areas and small failure probability. Fauriat and Gayton [19] presented the AKSYS method which is adaptation to system problems.
Many studies [20–23] indicated that it is the key problem to choose a good learning function for guiding the DOE’s strategy of KARM. Some efficient learning functions have been presented in recent years. Yang et al. [20] defines the expected risk function (ERF), which provides an indication of how large the predicted performance function value is if the sign of the predicted performance function is wrongly predicted. Lv et al. [21] proposed a new learning function based on information entropy theory and applied it to the reliability algorithm in combination with the LS method. Sun et al. [22] proposed a learning function that considers the degree of improvement in prediction accuracy and probability density based on a dynamic point selection update strategy and established a relatively highprecision Kriging prediction model relatively quickly. Wang and Sun [23] roughly determined the region where = 0 according to Kriging prediction statistics and gradually increased the number of sample points in the region to maximize a certain accuracy index of the Kriging model.
The aim of this paper is to further decrease the computation and improve the accuracy of reliability analysis. A Krigingbased active learning reliability method is presented, in which an efficient pointselected learning function (called HPF) is proposed to update DOE. The uncertainty of HPF could measure the features like the sampling density, the probability to be wrongly predicted, and the local and global uncertainties close to the limit state. The new pointselected learning function can help select the next point effectively, and only a small number of real performance function evaluations are required to build an accurate Kriging performance function. Finally, the gear reliability analysis is taken as the object of this study, whose performance function response’s behavior is nonlinear and timeconsuming. It is a difficult task that gear reliability is evaluated accurately and efficiently by traditional mechanical reliability methods (such as crude MCS, FORM, and RSM). The proposed method is applied to the reliability analysis of gear. This study will provide reference for the reliability design of the gear vibration.
This article is organized as follows: Section 2 introduces the structural reliability method based on the Kriging surrogate model and Monte Carlo simulation. Section 3.1 presents a new hybrid pointselected learning function and then describes the analytic expression of HPF. Section 3.2 describes the procedures of the proposed Krigingbased active reliability method. Section 4 illustrates the efficiency and correctness of the proposed method through two academic examples. In Section 5, the proposed method is applied to analyze the reliability analysis of gear. Section 6 is the conclusion.
2. Kriging Surrogate Model for Structural Reliability Analysis
Following [4], the Kriging response can be expressed aswhere β is the vector of regression coefficient; k(x) is the vector of polynomial, where constant polynomials are selected as the basic functions; and z(x) is the random process, and the covariance of z(x) iswhere σ^{2} is the variance of the process and R(θ, w, x) is the correlation function between z(w) and z(x), where a Gaussian correlation function is selected as
Under the consideration of linear predictor, the predictor at the unobserved point x is a normal distribution variable before calculating the true response. And the predictor follows a Gaussian distribution:where the Kriging predictor value and standard variance can be expressed aswhere is the regression parameters of generalized least squares solution; R is the matrix of stochasticprocess correlations between design sites; and R = [R · ij] m × m, where R · ij = R(θ, s i, θ j), i, j = 1, ⋅⋅⋅, m; ; ; .
The optimal coefficients corresponding to maximum likelihood estimation should solvewhere R is the determinant of R.
When is found, a Kriging prediction model can be established.
Then, the Kriging surrogate expression will be substituted for the performance function to calculate the failure probability, and it is used to calculate the failure probability:
Because is explicit expression, the computer can complete the calculation of the performance function at millions of sample points in a short time. Therefore, Monte Carlo simulation, as the most common and robust method, is applied to reliability analysis in this paper. The failure probability can be expressed aswhere x = [x_{1}, x_{2}, ⋅⋅⋅, x_{n}]^{T} is the vector of input random variables with ndimension in the domain; f_{X}(x) is the joint probability density function of x; n_{MCS} is the number of random samples generated by MCS; x_{i}, (i = 1, 2, … , n_{MCS}) is the ith sample in the Monte Carlo population; is a indicator function, and it can be defined as
The coefficient of variation of is calculated as
3. KrigingBased Active Learning Reliability Method
3.1. PointSelected Learning Function
The idea of the Krigingbased active learning reliability method (KARM) is to perform a relatively small DOE and sequentially add a new point step by step to the DOE, and then the constructed Kriging surrogate model gradually approaches to the real performance function. So the surrogate model is to be replaced with the real performance function when reliability analysis is conducted. Studies show that the key of the Krigingbased active learning reliability method is to choose an efficient pointselected learning function, where some widely used learning functions consist of such as EFF [11], U [13], and ERF [15]. Inspired by the own advantages of these efficient learning functions, combination of some efficient features is proposed and deduced to further improve the efficiency and accuracy of the KARM. The proposed learning function (called “HPF”) could consider the features like the sampling density, the probability to be wrongly predicted, and the local and global uncertainty close to the limit state. And, the derivation of HPF is elaborated below.
Firstly, based on equation (5) of Kriging theory in Section 2, is assumed to a random variable and under the Kriging framework. As shown in Figure 1, is the probability density and and are the Kriging predictor value and standard variance, to equation (1), the failure probability p_{f} is only related to the positive or negative of the state function instead of the absolute value of in the X domain. The accuracy of the Krigingbased reliability analysis can be interpreted as the accuracy of = 0 relative to = 0. As shown in Figure 2, and = 0 divide the X space into four parts: the signs of and in the domain D_{1} and D_{3} are the same and the signs of and in the domain D_{2} and D_{4} are converse. The probability that the sign of is opposite to the sign of can be expressed as
In the Kriging framework, the predicted sign at the point x is a random variable , and obeys the distribution of 01 with the parameter . According to [24], the accuracy of can be measured as
It is shown from equation (14) that the larger the value of is, the larger is the accuracy of probability of predicting to be wrong, where is the learning function of AKMCS [13] and indicates the sampling points’ density; so will be regarded as one part of the proposed indication function.
Secondly, as shown in Figure 3(a), the larger the value of [15] is, the higher the uncertainty of the corresponding sample point in the global area is. According to [11], the point with the larger value of is to be the higher uncertainty in local region , which is shown in Figure 3(b). Therefore, could balance the local and global uncertainty of , and then it measures the uncertainty more reasonably.
(a)
(b)
However, is a random variable, so the expectation of is used to measure this uncertainity. Just like EFF [11] and ERF [15], we average throughout the real space. So is expressed as
This integral is expressed analytically as (the derivation process can be seen in Appendix)where is the sign function and
Above all, a hybrid pointselected function (HPF) is defined aswhere and are the Kriging predictor value and standard variance seen in equations (6) and (7), respectively; Φ(⋅) and φ(⋅) are the cumulative distribution function and the probability density function of the standard normal distribution, respectively; and k = 3 is selected in this paper; is the sign function, and Note that equation (17) is hereafter referred to as ‘HPF’ in this study because HPF is the abbreviation of the hybrid pointselected function.
3.2. Stopping Criterion
The stopping criterion directly determines the number of calls to the real performance function (referred to as ‘n_{call}’). If n_{call} is too small, the accuracy of the constructed surrogate model is too low; while if n_{call} is too large, the computation is too large. Therefore, it is the key point to choose a precise and efficient criterion in the active learning reliability methods.
In AKMCS [13], the stopping criterion can be expressed aswhere U(x) is the learning function of AKMCS [13], , S is the sampling points of MCS, and .
The stopping criterion only considers the sampling point with the worst predicting accuracy. It does not consider the prediction accuracy of all sampling points of MCS. This stopping criterion is too rough and not rigorous. In order to more accurately predict the sign of all sampling points, a stopping criterion considering global samples is presented in the literature [25], namely,where the event C denotes all the points of the MCS population are correctly classified, and .
The stopping criterion means that the probability that all the points are correctly classified is greater than 0.99. It requires that the prediction accuracy of every point of the MCS population is more accurate. For example, if the prediction accuracy of every point of the 10^{5} MCS population is assumed as P(C_{i}) = 0.99999, the probability that all points are correctly classified should be . Or there is only one sampling point with P(C_{i}) < 0.99; the active learning process will not stop. Therefore, this stopping criterion is so strict that the accuracy of the calculated failure probability is very high, but it also requires more iterations. Especially when the real performance function is timeconsuming, the calculation of the active learning reliability method is too large.
In order to avoid above problems of the classical stopping criteria, an improved stopping criterion is adopted in this paper. This criterion is based on the fact that the signs of all sampling points are independent, and it could consider the global sampling points as well. This stopping criterion is not affected by a few local points with the low prediction accuracy (such as P(C_{i}) < 0.99). Therefore, it can get a relatively acceptable results faster than others. The improved stopping condition is presented as (the derivation process and more details can be found in the literature [17])where P_{threshold} is the threshold probability of classifying the right sign on S; if P_{threshold} = 0.998, it means that the probability of predicting a right sign on population S is 0.998 at least; (k = 1, 2, ..., n_{c}) are uniformly distributed in [P_{threshold},1), and the larger the n_{c} is, the more accurately the stopping condition (20) will judge the state of population (n_{c} equals to 100 in this paper); , where is the vector corresponding to , e.g., .
To better express how to use the stopping criterion, the judgment process of the stopping criterion is illustrated as follows. After each updating the Kriging prediction model, the following process is conducted.
Firstly, the probabilities of making a right sign at all sampling points are calculated by the current Kriging prediction model, and , P^{i} = Φ(U(x^{i}))(i = 1, 2, ..., n_{MCS}).
Secondly, (k = 1, 2, ..., n_{c}) are uniformly selected from the interval [P_{threshold}, 1), and they are , , respectively.
Thirdly, are calculated. For example, in case of k = 1, the number of sample points whose P^{i} is greater than , where is the vector corresponding to , e.g., P^{i} = Φ(x^{i}) (i = 1, 2, ..., n_{MCS}). In the similar way, when k = 2, 3, ..., n_{c}, , , ... are calculated, respectively.
Finally, on the basis of (k = 1, 2, ..., n_{c}) and , respectively, judge whether , , ..., If any of is meet, then the iteration would stop.
3.3. Procedures in the Present Study
The procedure (Figure 4) of the proposed active learning reliability method can be described as follows: Step 1. Generate a Monte Carlo population in the design space. At this stage, none of them is evaluated on the real performance function , and they are regarded as the candidate points to be evaluated on the real performance if KARM requires it. Step 2. Define the initial DOE with small size (, in which n is the number of random variables over the bounds ± five standard deviations, and calculate corresponding to with the real performance function . Step 3. Construct an Kriging predicted function model from the and . Step 4. Identify the best point by the proposed HPF which is shown in equation (19). The HPF criterion of finding the best point is expressed as Step 5. Predict the performance function value , and estimate the failure probability and its coefficient of variation according to equations (10) and (12), respectively. Step 6. If the stopping condition (equation (20)) is satisfied, stop iteration and go to Step 7; otherwise, N = N + 1, and update the previous DOE with the best point , and go to Step 3. Step 7. If , stop iteration and go to Step 8; otherwise, the population S enlarges and go to Step 3. Step 8. Output the results.
4. Validation
4.1. Case 1 (Modified Rastrigin Function)
This explicit example taken from [13, 18] is used to validate the efficiency and correctness of the proposed method. And, its performance function is a modified Rastrigin function, which is highly nonlinear involving nonconvex and nonconnex domains of failure. And, the performance function can be expressed aswhere is the vector of input random variables and and are standard normal distributed random variables, and they are independent.
This example is performed by the active learning reliability method (which can be seen in Section 3.3). Firstly, 10 points are sampled by the Latin hypercube sampling method as the initial DOE, and 10^{5} points (called population S) generated by MCS are regarded as the candidate points. The failure probability estimated by crude MCS, which is calling the real performance function, is regarded as the “true” solution (called it in this paper). Figure 5 shows the iterative process of the evaluated . As the number of DOE’s points increases, the estimated failure probability converges to the true one . It takes about 281 points to obtain the failure probability with relative high accuracy. Figure 6 shows the initial DOE’s points, iterative DOE’s points, and the Kriging predicted function after the iteration of the proposed method run 281 times. From Figure 6, the Kriging predicted performance function is almost consistent with in the main domain, and most of the sample points’ signs are right (Figure 7).
In order to demonstrate the efficiency of HPF, the active reliability method is tested with the different learning functions (such as EFF, ERF, U, and HPF). To avoid the effect of other random factors, EFF, ERF, U, and HPF are performed on the same Monte Carlo candidate points and initial DOE. Figure 8 shows the iterative graphs of corresponding to by the different pointselected learning functions. By comparison, the convergence speed of HPF is faster than that of other learning functions, and HPF requires the least number of calls to the real performance function than other pointselected learning functions at the same precision level.
To demonstrate the influences of the different stopping criteria on the reliability results, , and the improved criterion equation (20) (P_{threshold} = 0.998) are, respectively, applied in the active learning reliability methods. The results of different methods are listed in Table 1, where is the estimated failure probability; δ is the coefficient of variation for the estimated failure probability; is the relative percentage error in comparison with the reference failure probability and it can be calculated by ; and is the number of calls to the real performance function.

It can be seen from Table 1 that the failure probabilities obtained by three iteration termination criteria are all quite accurate. The result of the global stopping criterion P(C) ≥ 0.99 is the most accurate one, and its relative error is almost 0. This is because P(C) ≥ 0.99 is the most strict one. However, the improved stopping criterion requires the least number of calls to the real performance function among three criteria, whose relative errors are all less than 0.0028 (). Therefore, the improved stopping condition could get a relatively accurate result faster than others.
4.2. Case 2 (Multidimensional Performance Function)
The second example is a multidimensional performance function, which is a classic illustration in [13, 25]. The explicit expression of the performance function is given aswhere . The means and standard deviations of the six random variables are listed in Table 2.

In order to demonstrate the efficiency of HPF, the active reliability method is tested with the different learning functions (such as EFF, ERF, U, and HPF). EFF, ERF, U, and HPF are performed on the same Monte Carlo candidate points and initial DOE. Figure 9 shows the iterative graphs of corresponding to by the different pointselected learning functions. By comparison, the convergence speed of HPF is faster than that of other learning functions, and HPF requires the least number of calls to the real performance function than other pointselected learning functions at the same precision level.
To demonstrate the influences of the different stopping criteria on the reliability results, , and the improved criterion equation (20) (P_{threshold} = 0.998) are, respectively, applied in active learning reliability methods. The results of different methods are listed in Table 3, where is the estimated failure probability; δ is the coefficient of variation for the estimated failure probability; is the relative percentage error in comparison with the reference failure probability and it can be calculated by ; and is the number of calls to the real performance function.

It can be seen from Table 3 that the failure probabilities obtained by three iteration termination criteria are all quite accurate. The result of the global stopping criterion P(C) ≥ 0.99 is the most accurate one, and its relative error is almost 0. This is because P(C) ≥ 0.99 is the most strict one. However, the improved stopping criterion requires the least number of calls to the real performance function among three criteria, whose relative errors are all less than 0.00104 (). Therefore, the improved stopping condition could get a relatively accurate result faster than others.
5. Application of Mechanical Reliability Analysis
Gear is one of the most important parts of the mechanical equipment. When the gears are engaged, the dynamic transmission error of the gear pair will be generated due to the random errors (such as manufacturing error and assembly error). Too large dynamic transmission error will cause the mechanical system to produce strong impact, vibration, and noise [26], so it is regarded as a failure mode of the gear transmission. In this section, the nonlinear vibration system (Figure 10) of a spur gear pair is set as the object of this study, and the proposed reliability method is applied to the reliability analysis of the gear vibration.
5.1. Response of the Gear
The differential equation of the gear vibration system is expressed as (the detailed procedure of derivation is seen in [27])where the parameters of equation (24) are expressed as follows: , , , , , , , , , , , , , , , and , where and are the masses of the gear and pinion, respectively, and is the equivalent mass of the gear pair; , , and are the stiffnesses corresponding to , , and ; is the gear meshing force; is dimensionless time, in which t is the time; and is the inherent frequency, where is the average value of the timevarying meshing stiffness between the gear and pinion. The static transmission error is converted to a sine function with the Fourier series as . In the same way, , , and is, respectively, expressed as , , .
The fourthorder Runge–Kutta numerical integration method edited by Matlab sofware is used to solve the vibration differential equation, and then the response of gear pair is obtained. It is noted that the process of computing the equation (24) is timeconsuming and it requires about 7.260 seconds to calculate the implicit performance function. Furthermore, the response of this gear system is highly nonlinear.
5.2. Limit State Function of Gear Vibration
The meshing stiffness , meshing damping , tooth side clearance , and rotational speed of the gear vibration system are considered as random factors. According to [26], , , , and obey the Gauss distribution, and the means and standard deviations of all random variables are listed in Table 4.

If the fluctuation of the dynamic transmission error y(t) exceeds the threshold value, the gear vibration system is in the failure state. Based on this failure mode, the performance function of the gear vibration system is expressed aswhere is the fluctuation value of the dynamic transmission error, i.e., , in which and , respectively, represent the maximum and minimum values of the dynamic transmission error during the engagement of the gear pair, namely, and and is the threshold value, which is selected as 0.055 mm in this paper. When , , so the gear is in the safe state. When ,, and then the gear is in the failure state.
It is known from Section 5.1 that y(t) is an implicit response and needs to be obtained by solving equation (25) with the fourthorder Runge–Kutta numerical integration method, so the performance function is also an implicit response as well, and each calculation of the gear performance function requires about 7.260 seconds (i3 cpu). Furthermore, the response of the gear performance function is highly nonlinear. For example, Figure 10 is the bifurcation diagram of the gear vibration system’s response corresponding to the tooth side clearance b_{m}. It is shown from Figure 11 that as the tooth side clearance b_{m} increases; firstly, the system’s motion changes from single period (b_{m} ∈ [0.01, 0.032]) to five times period (b_{m} ∈ [0.32, 0.061]), and then it becomes the quasiperiodic motion (b_{m} ∈ [0.061, 0.068]); finally, it is to be the chaotic motion state. Namely, the motion behavior of this gear system is highly nonlinear. And, the performance function response of the gear vibration is very sensitive to the slight changes of the system’s parameters.
Therefore, this is the type problem of mechanical reliability analysis which is timeconsuming and nonlinear. Traditional mechanical reliability methods (such as MCS, FORM, and RSM) will become nonapplied because of the problems of large computation and low precision.
5.3. Results and Discussion
To validate the efficiency and correctness of the proposed method, the gear vibration reliability analysis is also performed by the traditional mechanical reliability methods (such as crude MCS, FORM [1], and RSM + MCS [5]) and the active learning reliability methods (such as HPF + , HPF + , and HPF + improved condition (P_{threshold} = 0.99)), where the result of the crude MCS is regarded as the true one; FORM is the first order secondmoment reliability method; RSM + MCS means that the polynomial response surface based on Latin DOE is established, and then the constructed response surface is used to calculate the failure probability using MCS.
Table 5 shows the results of different methods. It can be found from Table 5 that FORM, RSM + MCS, HPF +, , and HPF + improved condition cost the cuputime about 473.1 s, 4150 s, 2412 s, 3455 s, and 1827 s, respectively. The computational efficiency of FORM, RSM + MCS, and the active learning reliability methods (, , and HPF + improved condition) can be acceptable in the engineering. But the results of RSM + MCS and FORM are totally different from those of the crude MCS, and their relative errors are 0.9693 and 0.8831, respectively, due to the highly nonlinear performance function. Table 5 also shows that the results of the active learning reliability methods are all consistent with that of crude MCS, and the number of calls to the real performance function require 303(15 + 288), 416(15 + 401), and 240(15 + 225), respectively. And, they cost only 2412 s, 3455 s, and 1827 s cuputime. In comparison with 7.264 × 10^{5} s of the crude MCS, the proposed active learning reliability method can improve the efficiency and accuracy dramatically. Comparing with the two stopping criteria and , the improved stopping criterion requires the least number of calls to the real performance function among three criteria, and its accuracy is enough in the engineering (). The improved stopping condition could get a relatively accurate result faster than others. And, it is more suitable for reliability analysis with a timeconsuming performance function.

6. Conclusions
The aim of this study is to solve the problem of large computation and low precision for the mechanical reliability analysis which is timeconsuming and nonlinear response. A Krigingbased active learning reliability method is presented in this paper. The efficiency of the Krigingbased active learning reliability method mainly depends on the proposed hybrid pointselected learning function (HPF) which is inspired by EFF, ERF, and U. The analytic expression of HPF is deduced by averaging the hybrid indicator throughout the real space. The uncertainty of HPF could measure the features like the sampling density, the probability to be wrongly predicted, and the local and global uncertainty close to the limit state. The new pointselected criterion can help select the next point effectively, and only a small number of true function evaluations are required to build an accurate Kriging performance function. The efficiency and correctness of HPF are validated by two explicit examples. Finally, the proposed method is applied to reliability analysis of the gear vibration (involving timeconsuming and highly nonlinear response). By comparing with traditional mechanical reliability methods, the results show that the proposed method can solve the problems of large computation and low precision.
For many engineering problems, the response of the system is often calculated by the numerical method (timeconsuming) and its behavior is usually highly nonlinear. The proposed active reliability method is independent of the gear reliability analysis and treats the system response as a black box. Thus, the proposed method can be also used to solve other engineering problems to perform reliability analysis.
Appendix
Derivation Process of Equation (16)
When , can be expressed analytically aswhere φ(⋅) is the probability density function of the standard normal distribution.
In the same way, when , the expectation of can be expressed analytically as
Equations (A.1) and (A.2) can be written into a uniform expression aswhere is the sign function and
Derivation completed.
Data Availability
The data of the active learning algorithm are available from the corresponding author upon request. However, the application data of this research cannot be valid because the application example in this paper is based on the background of the national aerospace field, which involves some confidential data.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was financially supported by the Young Doctor Scientific Research Foundation of College (Grant no. 19YB27), State Key Laboratory of Robotics (Grant no. 2017Z18), and Liaoning Provincial Natural Science Foundation (Grant no. 2018010334301). Their financial supports are gratefully acknowledged.
References
 M. D. McKay, R. J. Beckman, and W. J. Conover, “Comparison of three methods for selecting values of input variables in the analysis of output from a computer code,” Technometrics, vol. 21, no. 2, pp. 239–245, 1979. View at: Publisher Site  Google Scholar
 R. E. Caflisch, “Monte Carlo and quasiMonte Carlo methods,” Acta Numerica, vol. 7, pp. 1–49, 1998. View at: Publisher Site  Google Scholar
 S. K. Au and J. L. Beck, “Important sampling in high dimensions,” Structural Safety, vol. 25, no. 2, pp. 139–163, 2003. View at: Publisher Site  Google Scholar
 I. Kaymaz, “Application of Kriging method to structural reliability problems,” Structural Safety, vol. 27, no. 2, pp. 133–151, 2005. View at: Publisher Site  Google Scholar
 V. J. Romero, L. P. Swiler, and A. A. Giunta, “Construction of response surfaces based on progressivelatticesampling experimental designs with application to uncertainty propagation,” Structural Safety, vol. 26, no. 2, pp. 201–219, 2004. View at: Publisher Site  Google Scholar
 J. Cheng, “Hybrid genetic algorithms for structural reliability analysis,” Computers & Structures, vol. 85, no. 1920, pp. 1524–1533, 2007. View at: Publisher Site  Google Scholar
 T. H. Lee and J. J. Jung, “A sampling technique enhancing accuracy and efficiency of metamodelbased RBDO: constraint boundary sampling,” Computers & Structures, vol. 86, no. 1314, pp. 1463–1476, 2008. View at: Publisher Site  Google Scholar
 B. Gaspar, A. P. Teixeira, and C. G. Soares, “Assessment of the efficiency of Kriging surrogate models for structural reliability analysis,” Probabilistic Engineering Mechanics, vol. 37, pp. 24–34, 2014. View at: Publisher Site  Google Scholar
 X. Shi, Â. Palos Teixeira, J. Zhang, and C. Guedes Soares, “Kriging response surface reliability analysis of a shipstiffened plate with initial imperfections,” Structure and Infrastructure Engineering, vol. 11, no. 11, pp. 1450–1465, 2015. View at: Publisher Site  Google Scholar
 D. R. Jones, M. Schonlau, and W. J. Welch, “Efficient global optimization of expensive blackbox functions,” Journal of Global Optimization, vol. 13, no. 4, pp. 455–492, 1998. View at: Publisher Site  Google Scholar
 B. J. Bichon, M. S. Eldred, L. P. Swiler, S. Mahadevan, and J. M. McFarland, “Efficient global reliability analysis for nonlinear implicit performance functions,” AIAA Journal, vol. 46, no. 10, pp. 2459–2468, 2008. View at: Publisher Site  Google Scholar
 B. J. Bichon, J. M. McFarland, and S. Mahadevan, “Efficient surrogate models for reliability analysis of systems with multiple failure modes,” Reliability Engineering & System Safety, vol. 96, no. 10, pp. 1386–1395, 2011. View at: Publisher Site  Google Scholar
 B. Echard, N. Gayton, and M. Lemaire, “AKMCS: an active learning reliability method combining Kriging and Monte Carlo Simulation,” Structural Safety, vol. 33, no. 2, pp. 145–154, 2011. View at: Publisher Site  Google Scholar
 B. Echard, N. Gayton, M. Lemaire, and N. Relun, “A combined importance sampling and Kriging reliability method for small failure probabilities with timedemanding numerical models,” Reliability Engineering & System Safety, vol. 111, pp. 232–240, 2013. View at: Publisher Site  Google Scholar
 X. Yang, Y. Liu, Y. Gao, Y. Zhang, and Z. Gao, “An active learning Kriging model for hybrid reliability analysis with both random and interval variables,” Structural and Multidisciplinary Optimization, vol. 51, no. 5, pp. 1003–1016, 2015. View at: Publisher Site  Google Scholar
 Z. Wen, H. Pei, H. Liu, and Z. Yue, “A Sequential Kriging reliability analysis method with characteristics of adaptive sampling regions and parallelizability,” Reliability Engineering & System Safety, vol. 153, pp. 170–179, 2016. View at: Publisher Site  Google Scholar
 C. Tong, Z. Sun, Q. Zhao, Q. Wang, and S. Wang, “A hybrid algorithm for reliability analysis combining Kriging and subset simulation importance sampling,” Journal of Mechanical Science and Technology, vol. 29, no. 8, pp. 3183–3193, 2015. View at: Publisher Site  Google Scholar
 F. Cadini, F. Santos, and E. Zio, “An improved adaptive Krigingbased importance technique for sampling multiple failure regions of low probability,” Reliability Engineering & System Safety, vol. 131, pp. 109–117, 2014. View at: Publisher Site  Google Scholar
 W. Fauriat and N. Gayton, “AKSYS: an adaptation of the AKMCS method for system reliability,” Reliability Engineering & System Safety, vol. 123, pp. 137–144, 2014. View at: Publisher Site  Google Scholar
 X. Yang, Y. Liu, Y. Zhang, and Z. Yue, “Probability and convex set hybrid reliability analysis based on active learning Kriging model,” Applied Mathematical Modelling, vol. 39, no. 14, pp. 3954–3971, 2015. View at: Publisher Site  Google Scholar
 Z. Lv, Z. Lu, and P. Wang, “A new learning function for Kriging and its applications to solve reliability problems in engineering,” Computers & Mathematics with Applications, vol. 70, no. 5, pp. 1182–1197, 2015. View at: Publisher Site  Google Scholar
 Z. Sun, J. Wang, R. Li, and C. Tong, “LIF: a new Kriging based learning function and its application to structural reliability analysis,” Reliability Engineering & System Safety, vol. 157, pp. 152–165, 2017. View at: Publisher Site  Google Scholar
 J. Wang and Z. Sun, “The stepwise accuracyimprovement strategy based on the Kriging model for structural reliability analysis,” Structural and Multidisciplinary Optimization, vol. 58, no. 2, pp. 595–612, 2018. View at: Publisher Site  Google Scholar
 J. Wang, Z. Sun, Q. Yang, and R. Li, “Two accuracy measures of the Kriging model for structural reliability analysis,” Reliability Engineering & System Safety, vol. 167, pp. 494–505, 2017. View at: Publisher Site  Google Scholar
 N. Lelièvre, P. Beaurepaire, C. Mattrand, and N. Gayton, “AKMCSi: a Krigingbased method to deal with small failure probabilities and timeconsuming models,” Structural Safety, vol. 73, pp. 1–11, 2018. View at: Publisher Site  Google Scholar
 Z. L. Sun and L. Y. Chen, Practical Mechanical Reliability Design Theory and Method, Science Press, Shanghai, China, 2003.
 R. Li and J. Wang, Dynamics of Gear System, Science Press, Shanghai, China, 1997.
Copyright
Copyright © 2019 Cao Tong 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.