Simulating a Weak Rock Mass by a Constitutive Model

This study presents a constitutive model for simulating the behavior of weak rock masses under various stress conditions, including the effects of pore pressure and temperature. Addressing the limitations of existing models in accurately representing the complex anisotropic response of these materials, the model utilizes Monte Carlo simulations to integrate stress anisotropy, pore pressure effects, and deviatoric stress states. This approach aims to capture the impact of geological factors such as foliation and jointing on the mechanical behavior of weak rock masses, which are often characterized by low strength and high deformability. Five rock types (claystone, mudstone, sandstone, shale, and siltstone) were simulated, generating 1000 cases per type with variability modeled using Weibull distributions. Statistical validation, employing the Kolmogorov–Smirnov test and Q–Q plots, demonstrated a strong agreement between simulated and experimental data. The results suggest that the proposed model can effectively predict deformation patterns in weak rock masses, offering potential applications in mining, geothermal energy extraction, and other engineering projects involving these challenging geological formations.