The　assumed　enhanced　strain　(AES)　method　has　been　successfully　used　to　deal　with　the　strain　localization　or　shear　band　problems.　With　AES　treating　the　localization　deformation　as　strong　discontinuity　in　displacement　field,　the　structural　phenomenon　and　the　associated　strength　loss　manifested　in　strain　localization　could　be　reflected　in　a　way　independent　of　mesh　density　and　alignment.　On　the　other　hand,　the　numerical　manifold　method　(NMM)　excels　at　dealing　with　non-homogeneity　and　discontinuity　problems　because　of　the　two　sets　of　covers.　By　combining　the　NMM　with　the　AES　method,　this　paper　could　investigate　the　influence　of　soil　non-homogeneity　on　shear　band　evolution　naturally　and　mesh　independently.　The　Gauss　Seidel　projection　contraction　(GSPC)　method　is　invoked　because　it　is　simple　and　efficient　in　tackling　the　integration　difficulties　induced　by　cutting　the　tensile　part　of　the　Drucker　Prager　(DP)　model.　Finally,　this　paper　studies　the　influences　of　erratic　boulders　and　weak　soil　layer　on　the　evolution　of　shear　band　in　the　slope.