A Decomposed Karhunen–Loève Expansion Scheme for the Discretization of Multidimensional Random Fields in Geotechnical Variability Analysis

Bin Zhu,Tetsuya Hiraishi

DOI: https://doi.org/10.1007/s00477-023-02625-8

IF: 3.821

2024-01-01

Stochastic Environmental Research and Risk Assessment

Abstract:The efficient discretization of a multidimensional random field with high definition and large geometric size remains a significant challenge. Compared with the simulation of one-dimensional or two-dimensional random fields, the generation of three-dimensional (3-D) random fields using the traditional Karhunen–Loève (K–L) expansion method tends to require a relatively long computational time and a larger amount of physical memory. In the present study, a decomposed K–L expansion scheme is proposed, which is applicable when a separable autocorrelation function (ACF) is used. The separability in the ACF is a precondition for the implementation of the proposed method. The proposed method decomposes the discretization of a 3-D random field into that of separate one-dimensional random fields and computes the eigenpairs in each dimension respectively. The accuracy and efficiency of the proposed method are demonstrated, and the numerical solutions of the eigenpairs used for the random field discretization are validated by comparing them with the theoretical solutions. Comparisons with the traditional K–L expansion method showed that the proposed scheme significantly reduced the computational time and memory requirements. This makes it potentially useful for the discretization of a multidimensional random field with a large geometric size and high definition, which can be used in stochastic finite element analysis. The proposed random field discretization method was applied to analyze a 3-D saturated clay slope under undrained conditions. The results of the stochastic finite element analysis suggest that the 3-D model is rather advantageous for estimating slope stability when the spatial variability of soil properties is considered. The 3-D model can capture more spatial details in soil distribution and its effects on failure mechanisms and safety factors from a theoretical point of view. Therefore, the proposed decomposed K–L expansion method for 3-D random field discretization has strong application potential in geotechnical problems.