In rocks or rock-like materials, the constituents, e.g. quartz, calcite and biotite, as well as the microdefects have considerably different mechanical properties that make such materials heterogeneous at different degrees. The failure of materials subjected to external loads is a cracking process accompanied with stress redistribution due to material heterogeneity. However, the latter cannot be observed from the experiments in laboratory directly. In this study, the cracking and stress features during uniaxial compression process are numerically studied based on a presented approach. A plastic strain dependent strength model is implemented into the continuous numerical tool—Fast Lagrangian Analysis of Continua in three Dimensions (FLAC^3D), and the Gaussian statistical function is adopted to depict the heterogeneity of mechanical parameters including elastic modulus, friction angle, cohesion and tensile strength. The mean parameter μ and the coefficient of variance (h_cv, the ratio of mean parameter to standard deviation) in the function are used to define the mean value and heterogeneity degree of the parameters, respectively. The results show that this numerical approach can perfectly capture the general features of brittle materials including fracturing process, AE events as well as stress-strain curves. Furthermore, the local stress disturbance is analyzed and the crack initiation stress threshold is identified based on the AE events process and stress-strain curves. It is shown that the stress concentration always appears in the undamaged elements near the boundary of damaged sites. The peak stress and crack initiation stress are both heterogeneity dependent, i.e., a linear relation exists between the two stress thresholds and h_cv. The range of h_cv is suggested as 0.12 to 0.21 for most rocks. The stress concentration degree is represented by a stress concentration factor and found also heterogeneity dominant. Finally, it is found that there exists a consistent tendency between the local stress difference and the AE events process.

Forty-six ancient Tibetan star-shaped towers and a village are located on a giant slope, which would be partially flooded by a nearby reservoir currently under construction. Ground survey, boreholes, and geophysical investigations have been carried out, with results indicating that the slope consists of loose deposit with a mean thickness of approximately 80 m in addition to an overlying bedrock of micaceous schist and phyllite. Ground survey and Interferometric Synthetic Aperture Radar (InSAR) indicated that the slope is experiencing some local deformations, with the appearance of cracks and occurrence of two small landslides. Through using borehole logs with the knowledge of the regional geological background, it can be inferred that the loose deposit is a result of an ancient deep-seated translational landslide. This landslide was initiated along the weak layer of the bedding plane during the last glaciation in the late Pleistocene (Q₃) period, which was due to deep incision of the Dadu River at that time. Although it has not shown a major reaction since the ancient Tibetan star-shaped towers have been built (between 200 and 1600 AD), and preliminary studies based on geological and geomorphological analyses incorporated with InSAR technology indicated that the landslide is deformable. Furthermore, these studies highlighted that the rate of deformation is gradually reducing from the head to the toe area of the landslide, with the deformation also exhibiting relationships with seasonal rainstorms. The state of the toe area is very important for stabilizing a landslide and minimizing damage. It can be expected that the coming impoundment of the reservoir will increase pore pressure of the rupture zone at the toe area, which will then reduce resistance and accelerate the deformation. Future measures for protection of the slope should be focused on toe erosion and some bank protection measures (i.e., rock armor) should be adopted in this area. Meanwhile, some long-term monitoring measures should be installed to gain a deep understanding on the stability of this important slope.