Anisotropic Strength Identification of Sandy Cobble Soil Using Multiscale Limit Analysis

Sandy cobble soil is a composite made of soil matrix and cobbles, and the estimation of its shear strength always requires expensive large-scale experiments. The strength of the sandy cobble soil exhibits macroscopic anisotropy with respect to the direction of the major stress due to the observed dominant distribution of the cobble dip angle. In the present paper, a numerical homogenization procedure for anisotropic strength identification of the sandy cobble soils is established, which can take into account the influencing factors of the size, shape, and inclination of the cobbles and the mesoscopic strength of the soil–rock interface. To consider the condition of plain strain, the particle size distribution of the cross section of the stratum is derived based on the fractal theory and the transformation method of Walraven. The mesostructure of the sandy cobble soils is randomly produced using ellipses to model the cross section of the cobbles. An iterative procedure is utilized to represent the major stress orientation-dependent macroscopic strengths. The results are validated against the data from indoor experiments and global mesoscopic computations. It is shown that the macroscopic strength of the sandy cobble mixtures can be accurately determined and the iterative multiscale limit analysis method is reliable and efficient. Parameter analysis is finally conducted to discuss the effect of the mesoscopic properties on the macroscopic strength.