Abstract

This research aims to evaluate the calculation accuracy and efficiency of the artificial neural network-based important sampling method (ANN-IS) on reliability of structures such as drum brakes. The finite element analysis (FEA) result is used to establish the ANN sample in ANN-based reliability analysis methods. Because the process of FEA is time-consuming, the ANN sample size has a very important influence on the calculation efficiency. Two types of ANNs used in this study are the radial basis function neural network (RBF) and back propagation neural network (BP). RBF-IS and BP-IS methods are used to conduct reliability analysis on training samples of three different sizes, and the results are compared with several reliability analysis methods based on ANNs. The results show that the probability of failure of the RBF-IS method is closer to that of the Monte-Carlo simulation method (MCS) than those of other methods (including BP-IS). In addition, the RBF-IS method has better calculation efficiency than the other methods considered in this study. This research demonstrates that the RBF-IS method is well suited to structure reliability problems.

1. Introduction

Since brakes serve the essential purpose of slowing down and stopping vehicles, the reliability of brakes significantly impacts vehicle safety. Due to their advantages of large braking torque and low cost, drum brakes are widely used in automobiles [13]. Resonance occurs when the driving frequency is close to the natural frequency of the brake and can cause fatigue damage to the brakes and other components, potentially leading to failure of the braking system. Therefore, research on vibration reliability analysis of brake component is necessary.

Over the years, many scholars have contributed to significant research on structure reliability analysis. The most common methods used for reliability analysis include the first-order reliability method (FORM), second-order reliability method (SORM), and MCS. FORM [46] and SORM [68] of the approximation method are widely used for the reliability analysis due to the simplicity and efficiency of these methods. In addition, because the MCS method is very simple and easy to program, it is widely used for the reliability analysis of the engineering structures. However, this method requires a significant computation effort for low probability of failure (POF) problem [9]. In order to address this issue, several methods have been proposed to reduce computation effort, such as the IS [10, 11], directional importance sampling simulation method [12], and subset simulation method [13, 14].

Since the limit state function (LSF) of most engineering structures is generally either nonlinear or implicit, the FEA method is required to compute structure responses, such as deformations, frequencies, and stresses. FEA is often combined with reliability analysis methods to conduct the engineering structure reliability analysis. While it is easy to combine FEA with the MCS method to conduct reliability analysis, this method is computationally expensive. In recent years, metamodeling techniques have been developed to overcome this issue, such as the model tree (MT), evolutionary polynomial regression (EPR), multivariate adaptive regression spline (MARS), gene expression programming (GEP) [15], response surface method (RSM) [1618], support vector machine [19, 20], kriging surrogate model [2124], and ANN [2532]. Metamodeling techniques are adopted to establish the approximate models, which can replace the original implicit LSF. The ANN is widely used to solve reliability problems related to practical engineering structure. In this setup, FEA is used to obtain training sample of the ANN, the random variables of the structure are used for ANN input, the structure response is used for ANN output, and explicit form of the function relationship between the structure response and random variables can be established. The BP-based MCS method is used to predict the structure reliability, and the result is compared with three reliability analysis methods including the traditional MCS, polynomial-FOSM, and BP-based advanced first-order second moment method (AFOSM) [27]. The combination of an adaptive RBF metamodeling technique and a FORM proposes a new reliability analysis method for the practical tunnel engineering problems [28]. The uniform design method is used to improve the disadvantage of the conventional ANN-based RSM method [29]. The development and use of ANN-based reliability analysis methods is surveyed to propose the improvement techniques of the ANN training set [30]. Two ANNs (including a multilayer perceptron network and RBF) based several reliability analysis methods is used to conduct a reliability analysis of the laminated composite structures [31]. In the ANN-based reliability analysis method, the most of the calculation time is spent conducting the FEA, and the ANN sample data size has a significant influence on the calculation efficiency of the ANN-based reliability analysis method.

This study evaluates the calculation accuracy and efficiency of the ANN-IS method on reliability of structures. Using two numeric verification problems and a drum brake vibration reliability problem as examples, the RBF-IS method and BP-IS method are used to conduct the reliability analysis on the training samples equal to 100, 200, and 300 units. The calculation results are compared with several methods including MCS, AFOSM, BP-MCS, BP-AFOSM, RBF-MCS, and RBF-AFOSM.

This study is structured as follows: Section 2 presents a vibration reliability analysis based on the ANN-IS method, Section 3 provides the numerical analyses and results, and Section 4 presents the discussion and conclusions.

2. Vibration Reliability Analysis Based on the ANN-IS Method

2.1. Basis Theory of the Vibration Reliability Analysis

Resonance occurs when the driving frequency is close to the natural frequency of a structure, such that the structure enters a quasifailure state [33, 34]. According to reliability theory, the state function of the structure with vibration is defined aswhere f is the natural frequency of the structure, and is the driving frequency.

Based on the relationship between the structure natural frequency f and the driving frequency , the state function of the quasifailure can be expressed aswhere is the specified range, in which the value is generally between 0.1 and 0.15 times the natural frequency of structure.

The mean value and variance of the state function can be expressed as

The quasifailure probability of the structure is expressed as

If the natural frequency and the driving frequency are independent normally distributed, then the probability of failure can be expressed as

2.2. Importance Sampling Method

The importance sampling method is a commonly improved MCS method. This method is widely used because of its high sampling efficiency and small variance. The basic theory of this method is that the importance sampling density function is adopted to replace the original probability density function. The POF base on the IS method is computed aswhere is the POF of the reliability performance function, is the importance sampling density function, is the joint probability density function of random variables, is the m-dimension variable space, and is the indicator function of failure domain.

The estimated value of POF is computed as

The mean value, variance, and coefficient of variation (COV) of the estimated failure probability are computed as

2.3. Artificial Neural Network

When the reliability analysis method is used to for a complex structure, the responses of the structure need to be calculated by sophisticated numerical methods, such as FEA. However, this method is not suitable for reliability analysis of a complex structure because of its long computation time. In order to reduce the computation time, the ANN is commonly used to establish the LSF.

2.3.1. BP Neural Network

BP neural network is a widely used type of the ANN. The structure of a BP neural network is shown in Figure 1.

In this study, the BP model includes a hidden layer, an input layer, and an output layer, where the design parameters (material properties and structural dimensions) of drum brake are used for input data, and the natural frequency of the drum brake is used for output data. According to references [27, 35, 36], the number of neurons in the hidden layer is calculated aswhere is the neuron number of the hidden layer, n is the neuron number of the input layer, m is the neuron number of the output layer, and a is an adjusting constant; according to references [27, 35, 36], the range of a is the integer between 1 and 10.

The output model of the hidden layer neuron is computed aswhere IW is the weight matrix between the input layer and the hidden layer, b is the threshold matrix, X is the vector of input variables, and f () is the neuron transfer function between the input layer and hidden layer. Logsig and tansig functions are widely used in the transfer function.

The output model of the output layer neuron is computed aswhere LW is the weight matrix between the hidden layer and the output layer, b is the threshold matrix, and f is the neuron transfer function between the hidden layer and output layer. Purelin function is widely using in the transfer function.

If the difference between the model output value and the expected value exceeds the expected error value, gradient search technique is used to adjust the weight and threshold values, and the output value is recalculated. This process is repeated until the difference between the model output value and expected value is found to be less than the expected error value.

2.3.2. RBF Neural Network

RBF is a two-layer forward neural network. The RBF model includes the input layer, output layer, and one hidden layer. The RBF neural network structure model is shown in Figure 2.

RBF typically uses the Gaussian function. The activation function is computed aswhere is the Euclidean norm, is the Gaussian function center, is the Gaussian function variance, and is the hidden layer neuron number.

The output function of the RBF model is computed aswhere is the input parameters vector, is the output layer neuron number, and is the weight matrix between the hidden layer and the output layer.

The variance of basis function is computed aswhere is the expected output value, and is the sample number. RBF has a faster training speed.

2.4. Vibration Reliability Analysis Method of the Drum Brake Based on ANN-IS

The ANN-IS method is used to conduct the drum brake vibration reliability analysis based on the theory in Section 2.1 using ANN samples of different sizes. The steps of the method are given (for a detailed flowchart of the ANN-IS method, Figure 3.) Pseudocode of the proposed method can be viewed from the Supplementary Data 1.(1)Define the ANN sample data size . In this study, (2)The design of the experiment method is used to generate ANN sample points, and FEA is adopted to obtain the natural frequency of the drum brake at the training sample points.(3)According to the relation criterion in which the difference between the natural frequency and driving frequency of the drum brake could not exceed the specific value, the ANN is used to establish the LSF. In this study, the specific value is defined 0.1 times of the natural frequency.(4)The optimization algorithm such as AFOSM is used to evaluate the most probable point of the LSF(5)The important sampling density function is constructed with the most probable point as the sampling center(6)Define the sample size . Start iteration , according to the important sampling density function, generate sample point of the input variables (7)Trained ANN is used to calculate the LSF value at the sample point , and indicator function value calculates using equation (7). The Num1 and Num2 are calculated.(8)If the iteration number , go to step 9. If not, iteration number , go to step 6.(9)The estimated value of the POF and COV of the POF is calculated using equation (8) and equation (9).

3. Numerical Analyses and Results

To evaluate the calculation accuracy and efficiency of the proposed reliability analysis method, three examples are used to conduct the reliability analysis including the drum brake problem and two verification problems from the reference literature. In this study, BP and RBF neural networks are used to establish the LSF, and the same sample data are used for two neural networks. ANN sample data size is defined as , and the ANN sample data are divided into training data and testing data of the ANN [3537] where 90% of ANN sample data are used for training data, and 10% of ANN sample data is used for testing data. In order to evaluate the error of ANNs more accurately, the same testing sample data are used to compare the mean square errors (MSEs) of three ANN sample data , in which the testing sample data have 10 sample points.

3.1. Case 1: First Verification Problem

The LSF with two random variables is as follows [17]:where is the random variables of standard normal distribution.

ANN sample data sizes are defined such as ANN sample points are generated to calculate the original LSF values. The trained ANNs are used to establish the LSF under three training samples of difference sizes, respectively. The same training sample is used for both RBF and BP. Figure 4 shows the comparison between the original LSF values and the predicted LSF values by ANNs. The maximum error point of the testing sample is the third sample point.

The ANN sample sizes are 100, 200, and 300, the mean square errors (MSEs) of RBF are 4.03 E − 7, 9.81 E − 10, and 1.82 E − 10; the MSEs of BP are 6.99 E − 6, 1.37 E − 6, and 4.44 E − 7, respectively. Figure 5 shows the comparison of the MSE between RBF and BP.

The result shows the MSE of RBF is smaller than that of BP. In addition, the MSE of RBF is already very small for an ANN sample size of 200. Define the number of samples, in this study, . Two ANN-IS methods and several reliability analysis methods including AFOSM, IS, BP-MCS, BP-AFOSM, RBF-MCS, and RBF-AFOSM are used to calculate the reliability indices on training sample of three difference sizes. The calculated results are compared with the MCS method. The result is given in Tables 13 and Figure 6.

The POF of the MCS, AFOSM, and IS is 3.58 E − 04, 4.05 E − 04, and 3.58 E − 04, respectively. The POF of the RBF-IS method is 3.59 E − 04, 3.59 E − 04, and 3.58 E − 04, and the relative errors are 0.3739%, 0.1528%, and 0.0058%, respectively. The POF of the BP-IS method is 3.61 E − 04, 3.52 E − 04, and 3.56 E − 04, and relative errors are 0.8292%, 1.5901%, and 0.4791%, respectively. The results show that the POF of the RBF-IS method is closer to that of MCS than that of the BP-IS method. In addition, the used of number of samples for RBF-IS is smaller than the ANN-based MCS method, so that the RBF-IS method has better calculation efficiency than the ANN-based MCS method.

3.2. Case 2: Second Verification Problem

The LSF with random variables is as follows [8]; the mean and standard deviation of the random variables is given as Table 4.where is the random variable of the normal distribution.

For Case 2, we use the same ANN sample data sizes used in Case 1. The trained RBF and BP are used to establish the LSF on the training samples of three difference sizes, respectively. Figure 7 shows the comparison between the original LSF values and the predicted LSF values by ANNs.

The MSEs of RBF are 1.72 E − 10, 2.55 E − 11, and 1.49 E − 11; the MSEs of BP are 5.88 E − 10, 1.59 E − 10, and 2.52 E − 11, respectively. Figure 8 shows the comparison of the MSE between RBF and BP.

The MSEs of RBF are the smaller than that of BP. Two ANN-IS methods and several reliability analysis methods including AFOSM, IS, BP-MCS, BP-AFOSM, RBF-MCS, and RBF-AFOSM are used to calculate the reliability indices on training samples of three difference sizes. The calculated results are compared with the MCS method. The result is given in Tables 57 and Figure 9.

The POF of the MCS, AFOSM, and IS is 3.383 E − 04, 3.211 E − 04, and 3.368 E − 04, respectively. The number of samples for MCS method is 1.0 E + 7 and that of ANN-based MCS is 1.0 E + 6. The POF of the RBF-IS method is 3.378 E − 04, 3.382 E − 04, and 3.383 E − 04, and relative errors are 0.1546%, 0.0279%, and 0.0136%, respectively. The POF of the BP-IS method is 3.356 E − 04, 3.391 E − 04, and 3.381 E − 04, and relative errors are 0.8020%, 0.3194%, and 0.0645%, respectively. The results of two verification problems show that the POF of the RBF-IS method is closer to that of the MCS method than that of the BP-IS method. In addition, the number of samples used for the RBF-IS method is smaller than that used for the ANN-based MCS method, demonstrating that the RBF-IS method offers better calculation efficiency than does the ANN-based MCS method. Through two verification problems, the calculation accuracy and efficiency of the proposed method are verified. Therefore, the proposed method is applied to the reliability analysis of the drum brake.

3.3. Application Problem: Vibration Reliability Analysis of the Drum Brake

Resonance will occur when the driving frequency is close to the natural frequency of the brake, and the resonance will cause fatigue damage to the brakes and other components and that may lead to failure of the braking system. In this study, the vibration reliability problem of the drum brake including random parameters is researched. The material properties and structural dimensions of drum brakes are random due to the influence of technical conditions, manufacturing and installation errors, manufacturing environment, and material characteristics. Assuming that the design parameters all obey the normal distribution, the coefficients of variation of the material properties and structural dimensions are determined to be 0.05 and 0.005. In this study, the density (Dd), the elastic modulus (Ed), and Poisson’s ratio (Pd) of the brake drum; the density (Ds), the elastic modulus (Es), and Poisson’s ratio (Ps) of the friction plate; the density (Dp), the elastic modulus (Ep), and Poisson’s ratio (Pp) of the brake shoe; width of the brake drum (Wd), the brake shoe (Ws), and the friction plate (Wp); and the initial angle (θ1) and envelop angle (θ2) of the friction plate are taken as design parameters.

In order to reduce the calculation workload, sensitivity analysis of the design parameters is conducted [15, 35]. Then the “Six sigma analysis” module in the Workbench software is used to conduct the sensitivity analysis on the design parameters of the drum brake. Design of the experiment method is used to generate 50 sample points, and then, modal analysis is conducted to obtain the influence degree of the design parameters on the natural frequency of drum brakes. Figure 10 shows the result of the sensitivity analysis.

Based on the sensitivity analysis results, in order to reduce the calculation workload, the parameters that have little influence on the natural frequency of the brake were screened. Finally, nine parameters were selected as design parameters. The mean and standard deviation of the design parameters is given in Table 8.

According to the theory of Section 2.1, the LSF of the vibration reliability on the drum brake is aswhere is the design parameters of the drum brake, and the design parameters are random; is the natural frequency of the drum brake; is the driving frequency of the drum brake and is the 51 Hz in this study; and the specified range is 0.1 times of the natural frequency.

Figure 11 shows the FEA result at the mean value of the random parameters.

The ANN sample data sizes are defined such as , and the design of the experiment method is used to generate training sample points, respectively. FEA is adopted to obtain the natural frequency of drum brake at the training sample points. The ANN sample data can be viewed from the Supplemental Data 2. Trained RBF and BP are used to establish the LSF, respectively. Figure 12 shows the comparison between the FEA values and the predicted LSF values by ANNs.

The MSEs of RBF are 2.76 E − 5, 3.95 E − 6, and 1.79 E − 6; the MSEs of BP are 1.05 E − 4, 2.28 E − 5, and 4.19 E − 6, respectively. Figure 13 shows the comparison the MSE between RBF and BP.

Two ANN-IS methods are used to calculate the reliability indices on training samples of three difference sizes and are compared with ANN-based several reliability analysis methods such as BP-MCS, BP-AFOSM, RBF-MCS, and RBF-AFOSM. The calculated results are given in Tables 911 and Figure 14.

The POF of the RBF-IS method is 1.184 E − 04, 1.185 E − 04, and 1.184 E − 04 and that of the BP-IS method is 1.187 E − 04, 1.185 E − 04, and 1.184 E − 04, respectively. The results show that the COV (POF) values of different methods are less than 3 E − 2, and the calculated results from different reliability analysis methods are basically consistent. The POF of the RBF-IS method relative to change ANN sample data size is smaller than those of other methods. In addition, the ANN-IS method has better calculation efficiency than the ANN-based MCS method from number of samples. The number of samples of the ANN-based MCS method is 1 E + 7 and that of the ANN-IS method is 1 E + 4.

4. Discussion and Conclusions

This research evaluated the calculation accuracy and efficiency of the ANN-IS method on reliability of structures such as drum brakes. By analyzing the results of two numeric verification problems and the drum brake vibration reliability problem, we arrive at the following conclusions.

First, the result of the RBF-IS method is closer to that of the MCS method than is the BP-IS method result. After reviewing the results of two numeric verification problems, the results calculated from the MCS method are found to be consistent with reference MCS results. The relative errors of the RBF-IS method for the two numeric verification problems are 0.0058% and 0.0136%, respectively, which are smaller than the results from other methods for the same ANN sample size.

Second, the RBF-IS method has a better calculation efficiency than the other methods used in this study. Through numeric analysis, find that use of the RBF-IS method results in smaller changes to the POF as the ANN sample size changes than when other methods are used. Additionally, the POF of the RBF-IS method at an ANN sample size of 200 is found to be closer to that of the MCS method than those of other methods (including BP-IS) at an ANN sample data size of 300. For reliability of structures such as drum brake, the time required for the FEA increases with both the ANN sample size and the complexity of the structural model. Since the RBF-IS method can use a smaller ANN sample size than other methods while still main training good calculation accuracy, the required calculation time of the RBF-IS method is smaller than those of other methods (including BP-IS). Finally, the number of samples used by the RBF-IS method is 1 E + 4 and those used by the ANN-based MCS method is 1 E + 7. This indicates that the RBF-IS method offers a better calculation efficiency than the ANN-based MCS method.

The disadvantage of the RBF-IS method is that an optimization algorithm such as AFOSM is required to calculate the most probable point of the LSF, which requires the explicit form of the LSF to be known in advance.

Based on this research, find that the RBF-IS method is well suited to reliability problems such as the drum brake vibration reliability problem.

Data Availability

The datasets used to support the findings of this research are included within the Supplementary Materials.

Additional Points

In this study, according to the theory of the vibration reliability analysis, the ANN-IS method is used to conduct the vibration reliability analysis on the drum brake. In addition, the calculation results of two ANN-IS methods (RBF-IS and BP-IS) are compared for structure reliability problems. This research demonstrates that the RBF-IS method is well suited to structure reliability problems such as the drum brake vibration reliability problem.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Z. Yang and U. S. Pak conceptualized the study; U. S. Pak and C. U. Kwon performed finite element analysis; U. S. Pak performed reliability analysis; Z. Yang, U. S. Pak, and C. U. Kwon reviewed and edited the article. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

The authors are grateful for the support from Chinese National Natural Science Foundation (U1710119 and U1708254) and Fundamental Research Funds for the Central Universities (N2003022).

Supplementary Materials

Supplementary Data 1. Pseudocode of the proposed method. Supplementary Data 2. ANN sample data and testing data of the drum brake. Table 1. FEA results of drum brake when ANN sample data size is 100. Table 2. FEA results of drum brake when ANN sample data size is 200. Table 3. FEA results of drum brake when ANN sample data size is 300. Table 4. FEA results of ANN testing sample data. (Supplementary Materials)