Heterogeneity is an inherent component of rock and may be present in different forms including mineral heterogeneity,geometrical heterogeneity,weak grain boundaries and micro-defects.Microcracks are usually observed in crystalline rocks in two forms:natural and stress-induced;the amount of stressinduced microcracking increases with depth and in-situ stress.Laboratory results indicate that the physical properties of rocks such as strength,deformability,P-wave velocity and permeability are influenced by increase in miciocrack intensity.In this study,the finite-discrete element method(FDEM)is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks using the proposed micro discrete fracture network(μDFN) approach.The characteristics of the microcracks required to create μDFN models are obtained through image analyses of thin sections of Lac du Bonnet granite adopted from published literature.A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-μDFN models indicate that micro-heterogeneity has a profound influence on both the mechanical behavior and resultant fracture pattern.An increase in the microcrack intensity leads to a reduction in the strength of the sample and changes the character of the rock strength envelope.Spalling and axial splitting dominate the failure mode at low confinement while shear failure is the dominant failure mode at high confinement.Numerical results from simulated compression tests show that microcracking reduces the cohesive component of strength alone,and the frictional strength component remains unaffected.Results from simulated Brazilian tests show that the tensile strength is influenced by the presence of microcracks,with a reduction in tensile strength as microcrack intensity increases.The importance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and its effects on the mechanisms leading to spalling damage near an underground opening are also discussed. Heterogeneity is an inherent component of rock and may be present in different forms including mineral heterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks are usually observed in crystalline rocks in two forms: natural and stress-induced; the amount of stressinduced microcracking increases with depth and in-situ stress. Laboratory results indicate that the physical properties of rocks such as strength, deformability, P-wave velocity and permeability are influenced by increase in miciocrack intensity. In this study, the finite element method ( FDEM) is used to model microcracks heterogeneity by introducing into a model sample sets of microcracks using the proposed micro discrete fracture network (μDFN) approach. These characteristics of the microcracks required to create μDFN models are obtained through image analyzes of thin sections of Lac du Bonnet granite adopted from published literature. A suite of two-dimensional laboratory tests including uniax ial, triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data. The FDEM-μDFN models indicate that micro-heterogeneity has a profound influence on both the mechanical behavior and resultant fracture pattern. Ann increase in the microcrack intensity leads to a reduction in the strength of the sample and changes the character of the rock strength envelope. Parting and axial splitting dominate the failure mode at low confinement while shear failure is the dominant failure mode at high confinement. Numerical results from simulated compression tests show that microcracking reduces the cohesive component of strength alone, and the frictional strength component remains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced by the presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. The importance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure en velope and its effects on the mechanisms leading to spalling damage near an underground opening are also discussed.