Shear strength of hydrate-bearing sediment is an essential parameter for assessing landslide potential of hydrate reservoirs under exploration conditions. However, the characteristics and simulation of this shear strength under varying dissociation conditions have not been thoroughly investigated. To this end, a series of triaxial compression tests were first carried out on sediments with varying initial hydrate saturations along dissociation pathways. Combining measured data with microscale analysis, the underlying mechanism for the evolution of shear strength in hydrate-bearing sediment was studied under varying partial dissociation pathways. Moreover, a shear strength model for hydrate-bearing sediment was proposed, taking into account the hydrate saturation and the unhydrated water content. Apart from the parameters derived from the hydrate characteristic curve, only one additional model parameter is required. The proposed model was validated using measured data on hydrate sediments. The results indicate that the proposed model can effectively capture the shear strength behavior of hydrate-bearing sediment under varying dissociation paths. Finally, a sensitivity analysis of the model parameters was conducted to characterize the proposed model. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

To overcome the limitations of microscale experimental techniques and molecular dynamics (MD) simulations, a coarse-grained molecular dynamics (CGMD) method was used to simulate the wetting processes of clay aggregates. Based on the evolution of swelling stress, final dry density, water distribution, and clay arrangements under different target water contents and dry densities, a relationship between the swelling behaviors and microstructures was established. The simulated results showed that when the clay-water well depth was 300 kcal/mol, the basal spacing from CGMD was consistent with the X-ray diffraction (XRD) data. The effect of initial dry density on swelling stress was more pronounced than that of water content. The anisotropic swelling characteristics of the aggregates are related to the proportion of horizontally oriented clay mineral layers. The swelling stress was found to depend on the distribution of tactoids at the microscopic level. At lower initial dry density, the distribution of tactoids was mainly controlled by water distribution. With increase in the bound water content, the basal spacing expanded, and the swelling stresses increased. Free water dominated at higher water contents, and the particles were easily rotated, leading to a decrease in the number of large tactoids. At higher dry densities, the distances between the clay mineral layers decreased, and the movement was limited. When bound water enters the interlayers, there is a significant increase in interparticle repulsive forces, resulting in a greater number of small-sized tactoids. Eventually, a well-defined logarithmic relationship was observed between the swelling stress and the total number of tactoids. These findings contribute to a better understanding of coupled macro-micro swelling behaviors of montmorillonite-based materials, filling a study gap in clay-water interactions on a micro scale. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Structures constructed on collapsible soil are prone to failure under flooding. Agro-waste like rice husk ash (RHA) and its geopolymer (LGR), consisting of lime (L), RHA, water glass (Na2SiO3), and caustic soda (NaOH), present a potential solution to address this issue. RHA and LGR were mixed up to 16% to improve the collapsible soil. Samples were remolded at optimal water content and maximum dry density for strength and collapsible potential tests. Unconfined compressive strength, deformation modulus, and soaked California bearing ratio exhibit exponential improvement with the inclusion of LGR. Additionally, for comparison of microstructural characteristics, analyses involving energy-dispersive X-ray spectroscopy (EDAX) and scanning electron microscope (SEM) were conducted on both virgin and treated specimens. LGR resulted in the emergence of new peaks of sodium silicates and calcium silicates, as indicated by EDAX. The formation of H-C-A-S gel and H-N-A-S gel observed in SEM suggests the development of bonds among soil particles attributed to geopolymerization. SEM reveals the transformation of the inherent collapsible soil from a dispersed and silt-dominated structure to a reticulated structure devoid of micro-pores following the incorporation of LGR. A numerical model was constructed to forecast the performance of both virgin and stabilized collapsible soils under pre- and post-flooding conditions. The outcomes indicate an enhancement in the soil's bearing capacity upon stabilization with 12% LGR. The implementation of 12% LGR significantly resulted in a lower embodied energy-tostrength ratio, emissions-to-strength ratio, and relatively lower cost-to-strength ratio compared to the soil treated with 16% cement kiln dust (CKD). (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/

In earthquake-prone areas, mountain tunnels often suffer from seismic damage when traversing active fault zones. To capture the seismic behavior of mountain tunnel under the action of active faults motion, the rate and state friction (RSF) relation is introduced to define the stick-slip dynamic behavior of a fault. The RSF relation is implemented in the finite element methods (FEMs). Numerical simulations of triaxial patch tests indicate that the RSF method can effectively capture the stick-slip dynamics. To reproduce the seismic damage to Daliang tunnel caused by slip of the Lenglongling fault, a three-dimensional (3D) numerical model including tunnel structure and plates of the fault is established. Seismic waves triggered by fault slip are then reproduced using the model. The simulation results show that the waves are dissipated while travelling and that their amplitudes decrease with depth. The failure of the tunnel lining is captured, and its seismic responses, including the displacement and strain of the structure, are extracted for various fault strike angles. The simulations are consistent with the observations, and it indicates that the movement of the simulated tunnel structure adjacent to the fault surface is significantly greater than those in the foot wall and in the middle of the fault. This study has the potential to provide a more direct means of understanding the seismic action of infrastructure induced by earthquakes. Seismic waves are no longer needed as input to the numerical simulation and instead, the earthquakes are generated by directly modeling the stick-slip motion of the fault. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

The creep phenomenon of inelastic deformation of surrounding rock may occur under the action of deep geological stress for a long period of time, potentially resulting in large-scale deformations or even instability failure of the underground engineering. Accurate characterization of the creep behavior of the surrounding rock is essential for evaluating the long-term stability and safety of high-level radioactive waste (HLW) disposal repositories. Although the laboratory creep tests of brittle undamaged rocks, such as granite, have been extensively performed, the creep characteristics of fractured surrounding rock under the multi-field coupling environment still require attention. In this study, a series of creep experiments was conducted on Beishan granite, which was identified as the optimal candidate surrounding rock for the disposal repository in China. The effects of various factors, including inclination angle of fractures, stress conditions, temperatures, and water content, were investigated. The experimental results show that the axial total strain increases linearly with increasing stress level, while the lateral total strain, axial and lateral creep strain rates increase exponentially. The failure time of saturated specimens fractured at 45 degrees and 60 degrees is approximately 1.05 parts per thousand and 0.84 parts per thousand of that of dry specimens, respectively. The effect of temperature, ranging from room temperature to 120 degrees C, is minimal, compared to the substantial variations in strain and creep rates caused by stress and water content. The creep failure of specimens fractured at 30 degrees is dominated by rock material failure, whereas the creep failure of specimens fractured at 60 degrees is dominated by pre-existing fracture slip. At a 45 degrees fracture angle, a composite failure mechanism is observed that includes both rock material failure and pre-existing fracture slip. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

During the excavation of large-scale rock slopes and deep hard rock engineering, the induced rapid unloading serves as the primary cause of rock mass deformation and failure. The essence of this phenomenon lies in the opening-shear failure process triggered by the normal stress unloading of fractured rock mass. In this study, we focus on local-scale rock fracture and conduct direct shear tests under different normal stress unloading rates on five types of non-persistent fractured hard rocks. The aim is to analyze the influence of normal stress unloading rates on the failure modes and shear mechanical characteristics of non-persistent fractured rocks. The results indicate that the normal unloading displacement decreases gradually with increasing normal stress unloading rate, while the influence of normal stress unloading rate on shear displacement is not significant. As the normal stress unloading rate increases, the rocks brittle failure process accelerates, and the degree of rocks damage decreases. Analysis of the stress state on rock fracture surfaces reveals that increasing the normal stress unloading rate enhances the compressive stress on rocks, leading to a transition in the failure mode from shear failure to tensile failure. A negative exponential strength formula was proposed, which effectively fits the relationship between failure normal stress and normal stress unloading rate. The findings enrich the theoretical foundation of unloading rock mechanics and provide theoretical support for disasters prevention and control in rock engineering excavations. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Thermal damage mechanisms are crucial in reservoir stimulation for enhanced geothermal system (EGS). This study investigates the thermal damage mechanisms in granite samples from the Gonghe Basin, Qinghai, China. The granite samples were heated to 400 degrees C and then cooled in air, water, or liquid nitrogen. The physical and mechanical properties of the thermally treated granite were evaluated, and microstructural changes were analyzed using a scanning electron microscope (SEM) and computed tomography (CT). The results indicate that cooling with water and liquid nitrogen significantly enhances permeability and brittleness while reducing P-wave velocity, strength, and Young's modulus. Specifically, liquid nitrogen cooling increased granite permeability by a factor of 5.24 compared to the untreated samples, while reducing compressive strength by 13.6%. After thermal treatment, the failure mode of the granite shifted from axial splitting to a combination of shear and tension. Microstructural analysis revealed that liquid nitrogen-cooled samples exhibited greater fracture complexity than those cooled with water or air. Additionally, acoustic emission (AE) monitoring during damage evolution showed that liquid nitrogen cooling led to higher cumulative AE energy and a lower maximum AE energy rate, with numerous AE signals detected during both stable and unstable crack growth. The results suggest that liquid nitrogen induces a stronger thermal shock, leading to more significant thermal damage and promoting the development of a complex fracture network during EGS reservoir stimulation. This enhances both the heat exchange area and the permeability of the deep hot dry rock (HDR) in EGS reservoirs. The insights from this study contribute to a deeper understanding of thermal damage characteristics induced by different cooling media and provide valuable guidance for optimizing deep geothermal energy extraction. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Aiming at mitigating the high risks associated with conventional explosive blasting, this study developed a safe directional fracturing technique, i.e. instantaneous expansion with a single fracture (IESF), using a coal-based solid waste expanding agent. First, the mechanism of directional fracturing blasting by the IESF was analyzed, and the criterion of directional crack initiation was established. On this basis, laboratory experiments and numerical simulations were conducted to systematically evaluate the directional fracturing blasting performance of the IESF. The results indicate that the IESF presents an excellent directional fracturing effect, with average surface undulation differences ranging from 8.1 mm to 22.7 mm on the fracture surfaces. Moreover, during concrete fracturing tests, the stresses and strains in the fracturing direction are measured to be 2.16-3.71 times and 8 times larger than those in the non-fracturing direction, respectively. Finally, the IESF technique was implemented for no-pillar mining with gob-side entry retaining through roof cutting and pressure relief in an underground coal mine. The IESF technique effectively created directional cracks in the roof without causing severe roadway deformation, achieving an average cutting rate and maximum roadway deformation of 94% and 197 mm, respectively. These on-site test results verified its excellent directional rock fracturing performance. The IESF technique, which is safe, efficient, and green, has considerable application prospects in the field of rock mechanics and engineering. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Rock masses are often exposed to dynamic loads such as earthquakes and mechanical disturbances in practical engineering scenarios. The existence of underground caverns and weak geological structures like columnar jointed rock masses (CJRMs) and interlayer shear weakness zones (ISWZs) with inferior mechanical properties, significantly undermines the overall structural stability. To tackle the dynamic loading issues in the process of constructing subterranean caverns, a programmable modeling approach was utilized to reconstruct a large-scale underground cavern model incorporating ISWZs and columnar joints (CJs). By conducting dynamic simulations with varying load orientations, the analyses focused on the failure patterns, deformation characteristics, and acoustic emission activity within the caverns. Results revealed that the failure modes of the underground caverns under dynamic loading were predominantly tensile failures. Under X-direction loading, the failed elements were mainly distributed parallel to the CJs, while under Y-direction loading, they were distributed parallel to the transverse weak structural planes. Furthermore, the dynamic stability of the overall structure varied with the number of caverns. The dual-cavern model demonstrated the highest stability under X-direction loading, while the single-cavern model was the least stable. Under Y-direction loading, the cavern stability increased with the number of caverns. Importantly, different weak structures affected the dynamic response of caverns in different ways; the CJRMs were the primary contributors to structural failure, while ISWZs could mitigate the rock mass failure induced by CJs. The findings could offer valuable insights for the dynamic stability analysis of caverns containing CJRMs and ISWZs. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

The deterioration of soft rocks caused by freeze-thaw (F-T) climatic cycles results in huge structural and financial loss for foundation systems placed on soft rocks prone to F-T actions. In this study, cementtreated sand (CTS) and natural soft shale were subjected to unconfined compression and splitting tensile strength tests for evaluation of unconfined compressive strength (UCS, qu), initial small-strain Young's modulus (Eo) using linear displacement transducers (LDT) up to a small strain of 0.001%, and secant elastic modulus (E50) using linear variable differential transducers (LVDTs) up to a large strain of 6% before and after reproduced laboratory weathering (RLW) cycles (-20 degrees C-110 degrees C). The results showed that eight F-T cycles caused a reduction in qu, E50 and Eo, which was 8.6, 15.1, and 14.5 times for the CTS, and 2.2, 3.5, and 5.3 times for the natural shale, respectively. The tensile strength of the CTS and natural rock samples exhibited a degradation of 5.4 times (after the 8th RLW cycle) and 2.7 times (after the 15th RLW cycle), respectively. Novel correlations have been developed to predict Eo (response) from the parameters quand E50 (predictors) using MATLAB software's curve fitter. The findings of this study will assist in the design of foundations in soft rocks subjected to freezing and thawing. The analysis of variance (ANOVA) indicated 95% confidence in data health for the design of retaining walls, building foundations, excavation in soft rock, large-diameter borehole stability, and transportation tunnels in rocks for an operational strain range of 0.1%-0.01% (using LVDT) and a reference strain of less than 0.001% (using LDT). (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).