SuSIE: Subset Simulation with Intrepid Exploration

2026-07-24

AI Summary

This study addresses the limitations of conventional subset simulation in accurately estimating failure probabilities when dealing with multiple failure regions, discontinuous, or highly nonlinear performance functions—challenges largely stemming from its reliance on the random-walk Metropolis sampler. To overcome these issues, this work proposes the first integration of the Intrepid MCMC sampler into the subset simulation framework, effectively mitigating sampling inefficiencies in multimodal, discontinuous, and high-dimensional failure domains. The proposed approach significantly enhances both the accuracy and robustness of failure probability estimation for complex reliability problems. Its superior performance is demonstrated across a range of benchmark examples spanning dimensions from 2 to 1003, confirming its effectiveness and scalability in practical applications.

Abstract

This work explores the challenges associated with the subset simulation framework - a well-established algorithm for estimating structural reliability - in settings involving multiple, possibly disconnected regions of failure, or involving discontinuous or sharply changing performance functions. We demonstrate that these drawbacks of subset simulation stem from the limitations of the Markov chain Monte Carlo (MCMC) sampler employed within the method, not from the subset simulation framework itself. Traditional random-walk Metropolis algorithms, which are known to struggle severely with sampling from multimodal distributions, are conventionally applied within subset simulation, which leads to inaccurate failure probability estimates in such cases. In this work, we instead utilize a modified version of the Intrepid MCMC sampler, which has recently been shown to be more effective than vanilla random-walk Metropolis algorithms in sampling from multimodal probability distributions. The subset simulation method with the proposed Intrepid sampler is demonstrated to address these complicating features. Several illustrative examples are considered, ranging from 2 to 1003 dimensions and exhibiting multiple failure regions or highly nonlinear performance functions, including both analytical problems and structural engineering applications.

Problem

• subset simulation
• structural reliability
• multimodal distributions
• failure probability estimation
• discontinuous performance functions

Innovation

• Subset Simulation
• Intrepid MCMC
• Multimodal Sampling
• Structural Reliability
• Failure Probability Estimation