Distinct element simulation in jointed rock masses is largely dependent upon the joint constitutive equation used. This paper describes the differences between the Barton-Bandits (BB) and the Mohr-Coulomb (MC) joint constitutive models for the stability analysis of the jointed rock slopes. The BB model, which allows the modelling of the dilation accompanying shear, predicts results very similar to the present condition of slopes. Consequently the 10 cm thick shotcrete was proposed for the reinforcement of those slopes. The MC model, however, in which the dilation angle is constant, is relatively insensitive to the behaviors of joints.