Rock is a very complex and heterogeneous material, containing structural flaws due to geologic generation process. Because of those structural flaws, deformation and failure of rock when subjected to differential compressive stresses is non-linear. To simulate the non-linear behavior of rock, mechanical crack models, that is, sliding and shear crack models have been used in several studies. In those studies, non-linear stressstrain curves and various behaviors of rock including the changes of effective elastic moduli (E₁, E₂, ν₁, ν₂, G₂) due to crack growth were simulated (Kemeny, 1993; Jeon, 1996, 1998). Most of the studies have mainly focused on the verification of the mechanical crack model with relatively less attempt to apply it to practical purposes such as numerical analysis for underground and/or slope design. In this study, the validity of mechanical crack model was checked out by simulating the non-linear behavior of rock and consequently it was applied to a practical numerical analysis, finite element analysis commonly used.