Simulating the effect of a variable normal stiffness on shear behavior of discontinuities