A new adaptive importance sampling scheme for reliability calculations

Abstract

An adaptive importance sampling methodology is proposed to compute the multidimensional integrals encountered in reliability analysis. It is based on a Markov simulation algorithm due to Metropolis et al. (Metropolis, Rosenbluth, Rosenbluth and Teller, Equations of state calculatons by fast computing machines. Journal of Chemical Physics, 1953;21(6): 1087-1092). In the proposed methodology, samples are simulated as the states of a Markov chain and are distributed asymptotically according to the optimal importance sampling density. A kernel sampling density is then constructed from these samples which is used as the sampling density in an importance sampling simulation. The Markov chain samples populate the region of higher probability density in the failure region and so the kernel sampling density approximates the optimal importance sampling density for a large variety of shapes of the failure region. This adaptive feature is insensitive to the probability level to be estimated. A variety of numerical examples demonstrates the accuracy, efficiency and robustness of the methodology.