Bearing Capacity of Spatially Random Cohesive Soil Using Numerical Limit Analyses

This paper describes a probabilistic study of the two dimensional bearing capacity of a vertically loaded strip footing on spatially random, cohesive soil using Numerical Limit Analyses (NLA‐CD). The analyses uses a Cholesky Decomposition (CD) technique with mid‐point discretization to represent the spatial variation in undrained shear strength within finite element meshes for both upper and lower bound analyses, and assumes an isotropic correlation length. Monte Carlo simulations are then used to interpret the bearing capacity for selected ranges of the coefficient of variation in undrained shear strength and the ratio of correlation length to footing width. The results are compared directly with data from a very similar study by Griffiths et al. in which bearing capacity realizations were computed using a method of Local Average Subdivision (LAS) in a conventional displacement‐based Finite Element Method (FEM‐LAS). These comparisons show the same qualitative features, but suggest that the published FEM calculations tend to overestimate the probability of failure at large correlation lengths. The NLA method offers a more convenient and computationally efficient approach for evaluating effects of variability in soil strength properties in geotechnical stability calculations