“…However, even if the performance function is explicit, the nonlinearity of the function and the non-Gaussian distribution of random variables always make it impossible to obtain the exact solution of this multi-fold integral [1], [2]. Various methods have been used to analyze the failure probability of structures subjected to random material parameters, such as the first-and second-order reliability methods (FORM/SORM) [3], [4], the simulation methods [5]- [7], the traditional perturbation stochastic finite element methods (PSFEM) [8], [9], the response surface method [10], [11], and so on. However, more or less problems exist in the above mentioned methods.…”