Collocation-based stochastic finite element analysis for random field problems

A stochastic response surface method (SRSM) which has been previously proposed for problems dealing only with random variables is extended in this paper for problems in which physical properties exhibit spatial random variation and may be modeled as random fields. The formalism of the extended SRSM is similar to the spectral stochastic finite element method (SSFEM) in the sense that both of them utilize Karhunen–Loeve (K–L) expansion to represent the input, and polynomial chaos expansion to represent the output. However, the coefficients in the polynomial chaos expansion are calculated using a probabilistic collocation approach in SRSM. This strategy helps us to decouple the finite element and stochastic computations, and the finite element code can be treated as a black box, as in the case of a commercial code. The collocation-based SRSM approach is compared in this paper with an existing analytical SSFEM approach, which uses a Galerkin-based weighted residual formulation, and with a black-box SSFEM approach, which uses Latin Hypercube sampling for the design of experiments. Numerical examples are used to illustrate the features of the extended SRSM and to compare its efficiency and accuracy with the existing analytical and black-box versions of SSFEM.