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.

Performance of Neumann Expansion Preconditioners for Iterative Methods with Geotechnical Elastoplastic Applications

Theoretical Elastoplastic Analysis for Foundations with Geosynthetic-Encased Columns

Robust partitioned block preconditioners for large-scale geotechnical applications with soil-structure interactions

ELASTIC-PLASTIC ANALYSIS OF 3D SOLIDS USING THE FINITE PARTICLE METHOD

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

GPU-accelerated iterative solutions for finite element analysis of soil–structure interaction problems

A robust and efficient iterative method for hyper-elastodynamics with nested block preconditioning

Elasto-plastic solutions for expanding cavities in strain-hardening and/or softening soils

Drained elasto-plastic solution for cylindrical cavity expansion in refuse soil

Study on 1-D Nonlinear Consolidation Behavior for Ningbo Soft Clay with Threshold Gradient

On the Preconditioners for Fast Multipole Boundary Element Methods for 2D Multi-Domain Elastostatics

Coupled thermomechanical analysis for cavity expansion under non-isothermal conditions

Generalized Neumann Expansion and Its Application in Stochastic Finite Element Methods

Symmetrization Method for the Generalized Plasticity Model with Nonassociated Plastic Flow Rule

Parsimonious Universal Function Approximator for Elastic and Elasto-Plastic Cavity Expansion Problems

2.5D finite element algorithm based on Drucker-Prager yield criterion

A robust stress update algorithm for elastoplastic models without analytical derivation of the consistent tangent operator and loading/unloading estimation

An adaptive model order reduction method for nonlinear seismic analysis of civil structures based on the elastic–plastic states

Fast Auxiliary Space Preconditioners for Linear Elasticity in Mixed Form

Finite Element Simulations of an Elasto-Viscoplastic Model for Clay

Stability and Deformation Analysis for Geotechnical Problems with Nonassociated Plasticity Based on Second-Order Cone Programming