[论文解读] The Geometry of Limit State Function Graphs and Subset Simulation

Karl Breitung|arXiv (Cornell University)|May 12, 2017

Probabilistic and Robust Engineering Design参考文献 13被引用 3

一句话总结

本文研究了极限状态函数（LSF）的几何结构，揭示了当底层假设（如单峰失效域或单一主导设计点）被违反时，子集模拟（SuS）可能产生不准确的结果。研究表明，SuS可能无法捕捉复杂失效域的拓扑结构，因此建议在模拟前结合FORM/SORM以检测结构特征，从而提高小概率估计的可靠性。

ABSTRACT

In the last fifteen the subset sampling method has often been used in reliability problems as a tool for calculating small probabilities. This method is extrapolating from an initial Monte Carlo estimate for the probability content of a failure domain found by a suitable higher level of the original limit state function. Then iteratively conditional probabilities are estimated for failures domains decreasing to the original failure domain. But there are assumptions not immediately obvious about the structure of the failure domains which must be fulfilled that the method works properly. Here examples are studied that show that at least in some cases if these premises are not fulfilled, inaccurate results may be obtained. For the further development of the subset sampling method it is certainly desirable to find approaches where it is possible to check that these implicit assumptions are not violated. Also it would be probably important to develop further improvements of the concept to get rid of these limitations.

研究动机与目标

• 识别子集模拟（SuS）中隐含的几何假设，这些假设在被违反时会导致失效概率估计不准确。
• 通过反例证明，即使样本量很大，当存在多个或非凸失效域时，SuS也可能失效。
• 主张SuS作为独立方法缺乏检测极限状态函数中结构性缺陷的诊断能力。
• 建议在模拟前通过FORM/SORM分析检测设计点和曲率信息，以补充SuS。
• 倡导从单纯概率计算的结构可靠性范式，转向对失效域结构的识别。

提出的方法

• 使用分析和数值示例，测试SuS在具有复杂几何结构的极限状态函数上的表现，包括多个设计点和非凸失效域。
• 应用标准子集模拟算法，通过迭代条件概率：P(F) = ∏P(F_{k+1}|F_k)，其中F_k = {g(u) < a_k}且a_k → 0。
• 将SuS结果与FORM/SORM近似值进行比较，以识别由未建模的几何复杂性引起的差异。
• 采用马尔可夫链蒙特卡罗（MCMC）采样方法，在每一子集层级生成条件样本。
• 分析LSF g(u)的互补累积分布函数（CCDF），以评估SuS是否捕获了所有显著的失效贡献。
• 建议将FORM/SORM作为预处理步骤，以检测多个beta点和曲率，从而验证SuS的适用性。

实验结果

研究问题

• RQ1 在何种几何条件下，子集模拟无法准确估计小失效概率？
• RQ2 多个或非凸失效域如何影响子集模拟的收敛性和准确性？
• RQ3 当极限状态函数具有复杂拓扑特征且未被该方法捕捉时，子集模拟的可信度在多大程度上受到影响？
• RQ4 FORM/SORM分析能否作为诊断工具，在应用子集模拟前检测失效域中的结构性缺陷？
• RQ5 是否需要从以概率为中心的可靠性分析，转向对极限状态曲面结构的感知分析？

主要发现

• 当失效域具有多个设计点或非凸几何结构时，即使样本量很大，子集模拟也可能产生误导性结果。
• 该方法假设存在单一主导失效模式和单峰结构，但在高维实际问题中可能不成立。
• SuS缺乏内部诊断能力，无法检测是否已捕获所有显著的失效贡献，从而导致未被察觉的错误。
• 当极限状态函数为黑箱时，SuS无法检测如多个beta点或曲率等结构性特征，而FORM/SORM可以。
• 本文结论认为，SuS不应作为独立方法使用；而应在模拟前通过FORM/SORM验证失效域的几何结构。
• 若无此类验证，SuS可能被应用于其根本不适用的问题，最终仍需完整蒙特卡罗计算才能获得准确结果。