The　three-dimensional　numerical　manifold　method　(3DNMM)　method　is　further　enriched　to　simulate　wave　propagation　across　homogeneous/jointed　rock　masses.　For　the　purpose　of　minimizing　negative　effects　from　artificial　boundaries,　a　viscous　non-reflecting　boundary,　which　can　effectively　absorb　the　energy　of　a　wave,　is　firstly　adopted　to　enrich　3DNMM.　Then,　to　simulate　the　elastic　recovery　property　of　an　infinite　problem　domain,　a　viscoelastic　boundary,　which　is　developed　from　the　viscous　nonreflecting　boundary,　is　further　adopted　to　enrich　3DNMM.　Finally,　to　eliminate　the　noise　caused　by　scattered　waves,　a　force　input　method　which　can　input　the　incident　wave　correctly　is　incorporated　into　3DNMM.　Five　typical　numerical　tests　on　P/S-wave　propagation　across　jointed/homogeneous　rock　masses　are　conducted　to　validate　the　enriched　3DNMM.　Numerical　results　indicate　that　wave　propagation　problems　within　homogeneous　and　jointed　rock　masses　can　be　correctly　and　reliably　modeled　with　the　enriched　3DNMM.