Abstract
Based on the incremental method of elastic-plastic mechanics and bilinear projection operators, by combining the strength reduction method with the φ-v inequality, this paper proposes a virtual element method strength reduction technique for slope stability analysis. The deformations of a homogeneous slope and heterogeneous slope are solved under different strength reduction factors, and the mesh dependency problem of the method is discussed. Numerical examples verify the correctness and effectiveness of the proposed method. The results demonstrate that due to the differences in the physical and mechanical properties of soil and rocks, stress becomes concentrated in the contact zone between soil and rocks, and thus, the plastic zone surrounds the rocks. Therefore, it is difficult to form regular connections in a plastic zone, such as a soil slope. This method can be used to analyse the stability of a stony soil slope and to study the effects of particle size, rock content, rock density, and rock spatial distribution on the mechanical behaviour of stony soil slopes.