Shear failure mode and strength of artificial rock joints are studied under different levels of joint compression. A quasi-static discrete element modeling is presented that reproduces the brittle material as a dense assemblage of deformable triangular particles. Using this, plaster-made samples with saw-tooth asperities are modeled and numerical results are compared with experimental measurement. The investigation approves that the proposed DEM is successful of predicting shear response of rock joints.