Based　on　the　incremental　method　of　elastic-plastic　mechanics　and　bilinear　projection　operators,　by　combining　the　strength　reduction　method　with　the　phi　-　nu　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.