Sparse polynomial chaos expansion for high-dimensional nonlinear damage mechanics

Finite Element Simulations in solid mechanics are nowadays common practice in engineering. However, considering uncertainties based on this powerful method remains a challenging task, especially when nonlinearities and high stochastic dimensions have to be taken into account. Although Monte Carlo Simulation (MCS) is a robust method, the computational burden is high, especially when a nonlinear finite element analysis has to be performed behind each sample. To overcome this burden, several “model-order reduction” techniques have been discussed in the literature. Often, these studies are limited to quite smooth responses (linear or smooth nonlinear models and moderate stochastic dimensions). This paper presents systematic studies of the promising Sparse Polynomial Chaos Expansion (SPCE) method to investigate the capabilities and limitations of this approach using MCS as a benchmark. A nonlinear damage mechanics problem serves as a reference, which involves random fields of material properties. By this, a clear limitation of SPCE with respect to the stochastic dimensionality could be shown, where, as expected, the advantage over MCS disappears. As part of these investigations, options to optimize SPCE have been studied, such as different error measures and optimization algorithms. Furthermore, we have found that combining SPCEs with sensitivity analysis to reduce the stochastic dimension improves accuracy in many cases at low computational cost.