Abstract
Uncertainty-based design and optimization is becoming one of the research focuses for high safety of complex engineering systems. Evaluating failure probability is essential and necessary in uncertainty-based design and optimization. Reliability analysis method using original Monte Carlo simulation usually shows good performance in evaluation accuracy. However, in situations for computing small failure probabilities, the calculation efficiency of Monte Carlo simulation is low generally. In this study, an enhanced Monte Carlo simulation method is utilized to solve the above challenge for assessment of the probability of rare failure events in uncertainty-based design and optimization. A mathematic example and a speed reducer design problem are given to illustrate the utilization of the proposed approach.
Keywords
Introduction
Unavoidable uncertainties exist in every aspect of practical engineering. Without the tolerance space for uncertainties, these variations cause the design solutions infeasible.1–5 In complex engineering systems, even a small variation may affect system performance by the transport between linking subsystems. To solve the above problem, uncertainty-based design and optimization (UBDO) is becoming one of research focuses.6–16 To reduce the computation cost in uncertainty analysis in UBDO, Sues et al. 17 utilize surrogate models which are created at the top level instead of the finite element analysis models. Sues and Cesare 18 propose a framework for UBDO that uncertainty analysis can be divided from the optimization process. Uncertainty analysis is performed initially, followed by the execution of optimization. Moreover, a current multi-stage conduction strategy is utilized to introduce existing uncertainty measure approaches into UBDO. In the study of Meng et al., 12 the possibility and probability analysis is introduced into UBDO, proposing the UBDO with mixed variables in order to deal with the multi-sources uncertainties at the same time. Zhang et al. 19 utilize the interval design information of variables to consider the lack of knowledge and propose the reliability evaluation approach to establish the reliability of complex systems.
Although many researches are finished in UBDO, there are still some problems requiring more investments. In many situations, high-quality performances, for example, long life, are needed. To obtain the accurate uncertainty estimation of safety, Monte Carlo simulation (MCS) is usually introduced in UBDO. This method is robust to the dimension and type of inputs.20,21 It also has the ability to give accurate evaluations of the probability if sufficient simulation samples can be given. MCS, however, usually cannot be utilized in calculating small probabilities directly. Higher reliability of safety means smaller probabilistic of rare failure events. In this situation, the information of rare simulation samples about failure is necessary for evaluating reliability, which requires a large number of simulation samples until failure events happen. The needs of the large number of simulation samples will use longer evaluation time of reliability assessment.
The UBDO problem considering the above challenges is dealt here. In this article, failure probability is denoted by a product of conditional probabilities of a series of involved failure events using the subset simulation strategy.
22
Then, in the conditional probability spaces, the original problem is changed through a sequence of simulations of more frequent events. A sequence of next samples in every failure event
UBDO and commonly used uncertainty analysis methods
Here, a UBDO problem can be denoted as
where
During the optimization process, system reliability can be guaranteed by probabilistic constraints. Theoretically, the reliability can be denoted using the cumulative distribution function of
In practical engineering, different approximation or simulation approaches are widely utilized for reliability analysis in equation (2). Sampling-based methods usually are robust to utilize. Generally, accurate reliability estimation can be obtained when a large number of simulations are introduced. However, this method is not suitable for many practical challenges where computationally expensive
The strategy of an enhanced MCS approach
The division of intermediate events
The original strategy of sequential Monte Carlo simulation (SMCS) is to change an initial probability evaluation problem into a series of conditional probability evaluation problems. Use
The enhanced metropolis algorithm and its application in SMCS
The enhanced metropolis algorithm is a powerful means which usually is utilized into performing uncertainty estimation.
27
In other words, random samples which have a certain distribution can be generated by this method. The significance of the metropolis method in this article is that if a simulation point follows a conditional distribution
For every
Step 1. Set
Given the value of
Step 2. Accept or reject
The procedure of SMCS
Based on the enhanced metropolis method, the strategy of SMCS can be given as follows:
Step 1. Obtain the value of
where
Step 2. Calculate each conditional probability by SMCS based on the metropolis method. In SMCS, the beginning points of the
where the conditional PDF of
Step 3. The value of
UBDO using SMCS in SORA (UBDO-SMCS-SORA)
The sequential optimization and reliability assessment (SORA) approach is originally developed to enhance the efficiency in reliability-based design 28 and introduced into UBDO.29–33 The basic idea of SORA is to decouple a UBDO problem into reliability analysis problem and deterministic design problem, respectively. Here, the idea of SORA is introduced and the strategy of UBDO-SMCS-SORA is to be discussed in detail.
Step 1. Solve the deterministic design problem in equation (9) and obtain the value of
Step 2. Perform MCS directly based on
Step 3. Construct a vector
Step 4. Solve the deterministic design problem with the modified constraints as shown in equation (10) and obtain
Step 5. Given
The probabilistic constrain
Step 6. Once the value of
which

The strategy of SORA.

The flowchart of UBDO-SMCS-SORA.
The speed reducer example
Here, a UBDO problem is introduced to show the effectiveness of the given approach. A structure design problem of a speed reducer is considered here, which is shown in Figure 3.34,35 The design solutions obtained by MCS-based design and optimization (DO-MCS) are utilized as reference value.

Design variables of the speed reducer design.
Here, aleatory uncertainties are associated with design variables and can be described by normal distribution.
The UBDO problem of speed reducer design is as follows
where
Here, the target failure probabilistic
Design solutions of the example.
DO-MCS: Monte Carlo simulation–based design and optimization; UBDO-SMCS: uncertainty-based design and optimization–sequential Monte Carlo simulation.
Conclusion
In this study, an enhanced MCS–based UBDO method is given to solve the efficiency problem for the assessment of the probability of rare failure events. To reduce the computation burden and enhance the reliability analysis efficiency, subset simulation is utilized and integrated within SORA. Subset simulation resolves the efficiency problem by changing original problem into the calculation of some larger conditional probabilities. The modified metropolis algorithm is introduced to obtain simulation samples in SMCS during optimization process. A structure design problem is given to illustrate the effectiveness of UBDO-SMCS-SORA. Compared with the DO-MCS, the proposed approach can provide solutions as accurate as the DO-MCS but using much less simulation samples.
Footnotes
Academic Editor: Dong Wang
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is partially supported by the National Natural Science Foundation of China (grant no. 51605047); the China Postdoctoral Science Foundation (grant no. 2016M602687); the Natural Science Foundation of Guangdong Province, China (Grant No. 611228787036); and the Open Research Fund of Key Laboratory of szjj2017-100 and szjj2013-03, Xihua University.
