Home Latest Articles
Latest Articles
  • A.S. Devapriya, T. Thyagaraj
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 4108-4124.

    Landfill facilities around the world are designed to protect the environment and public health by using impermeable liner systems that isolate the waste and leachate produced from the waste. However, the functionality of liners has been reported to be significantly compromised by environmental loading due to the seasonal climatic and physico-chemical changes that alter their volume deformation and hydraulic characteristics. Bentonite admixed natural soils are employed as liner materials if they meet the hydraulic conductivity requirement in their as-compacted state. However, limited studies addressed the effects of wet-dry cycles combined with chemical contamination on the volumetric and hydraulic behaviour of bentonite admixed natural soils. In this study, Indian red soil was ameliorated with 10%, 20%, and 30% bentonite by weight, and the mixtures were subjected to alternate wetting and drying cycles using distilled water, 0.4 M NaCl, and 0.4 M CaCl2 solutions. All red soil-bentonite specimens met the hydraulic conductivity design criterion of 1 × 10-7 cm/s in their as-compacted states. However, significant variation in hydraulic behaviour was observed at the end of the wet-dry cycles, particularly with chemical contamination. The microstructural examination through scanning electron microscopy (SEM) and mercury intrusion porosimetry (MIP) revealed an increase in macropores volume with wet-dry cycles and increase in the induced osmotic suction, which was found to be a key factor influencing the hydraulic conductivity.

  • Zhuoyu Duan, Dongxing Wang, Yong Zou, Hao Luo
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 4081-4093.

    To realize the soil reinforced through the carbonation of ternary binder under ambient pressure and mild conditions, the present study introduces triethanolamine (TEA), which serves as an effective carbonation accelerator. Through the unconfined compressive strength (UCS) test, the soft soil solidified with ternary eco-binder consisting of ground granulated blast-furnace slag (GGBS), metakaolin (MK), and calcium carbide residue (CCR), subjected to carbonation, is investigated. The effect of TEA on the carbonation of soil is evaluated by the UCS and the CO2 mineralization. This study clarifies the influence factors, including the initial water content, TEA dosage, binder constituent ratio, and content. The optimal binder constituent ratio for the strength growth and carbonation efficiency of carbonated soil is approximately 4:4:2 for GGBS, CCR, and MK, respectively. The incorporation of TEA at a low dosage (<0.15 %) enhances the strength of carbonated soil, whereas the high dosages impair the strength. The synergistic effect of TEA and carbonation further improves the strength and compressibility of soil. The soil with 1.5 % TEA carbonated for 7d exhibits a 44.8 % increase in strength compared to that without TEA, which is attributed to a 2.2-fold increase in carbonation efficiency. The addition of TEA accelerates the ion dissolution and CO2 dispersion, promoting the carbonation reaction in soft soil. Calcite and aragonite precipitate during carbonation, contributing to the strength development of soil. The carbonates phase difference and the pore structure density with different TEA dosages are also demonstrated to be the strength influence factors.

  • Yang Xiao, Chang Zhao, Qingyun Fang, Xiang He, Jian Chu, Hanlong Liu
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 4023-4033.

    Effective sealing of geological fractures is essential for subsurface stability and mitigating environmental risks such as groundwater contamination and inefficient CO2 sequestration. Enzymatically Induced Carbonate Precipitation (EICP) offers a promising bio-mediated approach due to its ability to fill and seal fractures. However, real-time precipitation patterns and clogging behavior of EICP under varying fracture and flow conditions remain poorly understood. This study employs a transparent fracture model with visualization to systematically investigate the effects of fracture aperture, flow conditions, and surface roughness on EICP-mediated sealing. Results indicate that fractures with narrower apertures promote tortuous finger-like flow paths, while wider-aperture fractures show more uniform deposition, with fewer but wider preferential flow paths. An appropriate injection rate around 1 mL/min ensures uniform precipitation and effective clogging, avoiding inlet clogging at lower rates (0.1 mL/min) and flushing effect reducing deposition at higher rates (10 mL/min). Additionally, rough fractures exhibit higher precipitation efficiency and greater permeability reduction, driven by their irregular surface geometry, which creates more deposition sites and complex flow compared to smooth fractures. Image processing reveals that precipitation patterns in rough fractures match closely with aperture distribution, compared to more concentrated deposition in smooth fractures. These findings provide insights for optimizing EICP-mediated fracture sealing, with implications for groundwater protection and geotechnical practices.

  • Sheng'ao Jia, De'an Sun, Pan Chen
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3934-3948.

    The soil-water retention and soil shrinkage characteristics are both crucial constitutive relations for unsaturated soils. Although existing research has explored the correlation between these two characteristics to some extent, the underlying mechanisms remain inadequately investigated. To investigate the correlation between the soil-water retention and soil shrinkage behavior, a series of soil-water retention and soil shrinkage tests is performed on compacted clays over a wide suction range (0-367 MPa). The test results show that the pore water in compacted clays is first expelled from large pores in low suction range. The drainage of pore water at low suctions is predominantly responsible for the phase of structural shrinkage in the soil shrinkage curve. The consistency between the characteristic transitional water contents in the soil shrinkage curve (SSC) and the inflection points in the soil-water retention curve (SWRC) is identified for all the compacted clays. The bimodal pore-size distributions (PSDs) of different clayey soils are obtained using the mercury intrusion porosimetry. The bimodal pore-size distribution characterization is the intrinsic factor in shaping the bimodal morphology in the SWRC over a wide suction range. The low proportion of micropores in clays is responsible to the indistinct zero-shrinkage stage of the SSC. The microstructure measured by the scanning electron microscope indicates the manifestation of aggregation effects during desaturation process. The results demonstrate that soil shrinkage is primarily caused by the contraction of inter-aggregate pores, rather than the evolution of intra-aggregate pores. The findings can greatly enhance the understanding of the soil-water retention and mechanical behavior of compacted clays in varying water content conditions.

  • Mingming Ren, Manchao He, Jie Hu, Hongru Li, Yuxiang Ding, Xinhao Miao, Hongyi Zhang
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3722-3741.

    Accurate extraction of rock mass discontinuity parameters is crucial for stability assessment and engineering safety. High-resolution remote sensing facilitates automated extraction, but its effectiveness relies heavily on precise normal estimation to ensure geometric reliability. Conventional methods struggle to preserve sharp features such as edges and corners, thereby reducing accuracy. To address this, we propose a normal estimation method based on local geometric adjustment that enhances feature extraction while maintaining sharp geometries. The approach consists of four steps: (1) classifying points, (2) applying normal and axial projections, (3) fitting segmentation lines via least squares, and (4) refining normals by optimizing local neighborhoods. The proposed method was evaluated on computer-aided design (CAD) models, real objects, and rock mass point clouds, and benchmarked against eight representative algorithms, including principal component analysis (PCA), 2-Jet PCA, Voronoi-based PCA, PCPNet, neural gradient function (NeuralGF), low rank representation (LRR), normal estimation via shifted neighborhood (NSN) and pair consistency voting (PCV). Experimental results demonstrate that our method achieves superior accuracy and efficiency, significantly improving structural plane extraction and ensuring better preservation of sharp geometric features.

  • Ningjie Li, Xinli Hu, Jian Wang, Junxiang Huang, Hongchao Zheng, Wei Li
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3675-3687.

    Traditional deterministic numerical simulation often has a poor prediction performance for landslide-induced wave run-up (LIWR) hazards, as it neglects the effects of uncertainty. The limitation for efficiently quantifying the uncertainties in primary parameters remains largely unsolved. In this study, we propose a probabilistic evaluation method, integrating the adaptive Kriging (AK) metamodel method and probability density evolution method (PDEM) based on generalized F-discrepancy. A Taylor expansion-based adaptive design strategy is applied to construct the global AK model over representative points generated by generalized F-discrepancy, thereby approximating the numerical physical response (i.e., maximum LIWR). Using these approximate responses, the PDEM is used to compute the exceedance probabilities that LIWR heights exceed elements at risk based on a construction of virtual time, and then a probabilistic criterion is introduced to classify hazard zones. The proposed method is demonstrated via two examples: Example Ⅰ, which possesses risk element (building), and Example Ⅱwith water-level variations. The results indicate that the proposed method has an acceptable performance (showing a 1.7 % difference in exceedance probability compared to Monte Carlo simulation with 50,000 samples) with low computation cost (requiring 284 deterministic analyses). For two specific scenarios in this study, the wave induced by the landslide exhibits a solitary-like leading wave. The proposed probabilistic method provides promising prospects for quantifying LIWR uncertainties, and is helpful for direct, efficient, and low-cost quantification assessment of cascading hazards.

  • Zhibo Zhang, Zhanbo Cheng, Xunchang Fei, Xueyu Geng, Kangda Wang, Ziwen Yuan, Zihou Liu, Wei Wu
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 4064-4080.

    Incineration bottom ash (IBA) holds attractive potential as a construction material, yet its shear behavior under cyclic loading remains insufficiently understood. This study comprehensively characterizes the monotonic and cyclic simple shear behavior of Singapore-derived IBA under constant volume conditions, with particular emphasis on its reuse potential in dynamic load-bearing applications. Key findings reveal that: (1) The material exhibits marked strain-hardening characteristics, demonstrating a density-dependent friction angle increment from 38.3° (loose state) to 42.5° (dense state). (2) Mechanical performance shows strong dependence on Si-Ca-Fe/Al ternary chemical composition and particle gradation characteristics. (3) Distinct failure modes emerge under different loading conditions - liquefaction dominates under unidirectional cyclic simple shear (UDCSS) conditions at low cyclic stress ratios (CSRs) and confining pressures, while bidirectional cyclic simple shear (BDCSS) loading induces cyclic mobility failure at elevated CSR levels, with corresponding cyclic resistance ratios (CRRs) showing a 30 % reduction in BDCSS compared to UDCSS configurations. (4) Pore pressure ratio (Ru) evolution follows a triphasic pattern: liquefaction failures exhibit rapid Ru acceleration in initial and tertiary phases (terminal Ru > 0.9), contrasting with cyclic mobility failures characterized by decaying Ru growth rates and lower terminal Ru values. (5) Notably, the established correlation between CRR and normalized shear wave velocity (Vs1) aligns closely with that of sand-gravel mixture with 5 % fines, which demonstrates the comparable cyclic load-bearing capacity of IBA to that of conventional construction materials. The study highlights the effect of load direction, particle size, and mineralogy in design applications and supports IBA's suitability for reuse in infrastructure subjected to dynamic loads.

  • Mingxi Chen, Zihan Fu, Feng Xiong, Jie Jiang, Qinghui Jiang
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3779-3803.

    Timely identification of accelerating precursors and performing reliable time-to-failure analysis are the key components in the management of slope failure risks. This study focuses on rock slope failures and proposes a framework for online identification of accelerating precursors and dynamic probabilistic prediction of failure time grounded in Bayesian inference. By integrating the Bayesian online change-point detection (BOCD) method with a typical dimensionless trend (TDT) model, the BOCD-TDT algorithm is first developed for online identification of acceleration events and their corresponding onset of acceleration (OA). Subsequently, a Bayesian approach is employed to estimate the parameters of the inverse velocity (INV) method, enabling the dynamic probabilistic prediction of slope failure time while quantifying observational and model uncertainties across different accelerating deformation stages. Building on this, the influence of starting point (SP) selection, trend update (TU), and multi-data fusion on prediction reliability is evaluated, and a novel decision criterion for impending slope failure is proposed. The feasibility of the proposed methods is then validated using 73 rock slope failure cases. Results show that using INV data, the BOCD-TDT algorithm can reliably identify acceleration events and the corresponding OA. In time-to-failure analysis, the reliability of dynamic failure predictions can be enhanced by incorporating both observational and model uncertainties corresponding to the deformation stages into the Bayesian prediction model, along with TU detection and multi-data fusion. The proposed failure probability criterion provides valuable guidance for the identification of impending failure and the establishment of ultimate alert thresholds.

  • Jiasen Liang, Xueming Du, Lei Wang, Hongyuan Fang, Xiaohua Zhao, Bin Li, Kejie Zhai, Shanyong Wang
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3703-3721.

    Water-rich sand layers are frequently encountered as adverse geological conditions during underground construction. Polymer slurry grouting has been widely recognized as an effective technique for reducing permeability and enhancing the stability of such strata. In this study, a mathematical model is established to describe the diffusion behavior of polymer slurry in porous media under dynamic water conditions and is further validated through laboratory experiments. The theoretical formulation of the slurry permeation process is developed based on Darcy's law, the Hagen-Poiseuille flow principle, and the physicochemical characteristics of the slurry. The derivation primarily focuses on analyzing the dynamic response of the slurry under the influence of water flow, considering the effects of flow velocity, grouting pressure, and sand-layer porosity on diffusion behavior. To verify the proposed model, a visualized grouting simulation system was designed to observe the diffusion process of polymer slurry in water-rich sand layers. The results demonstrate that slurry diffusion is significantly affected by grouting pressure, porosity, and water flow velocity. The observed staged diffusion characteristics, dynamic evolution patterns, and directional effects are in good agreement with theoretical predictions. Furthermore, the average relative deviations between the theoretical and experimental results for diffusion pressure and diffusion distance are both less than 25 %, confirming the reliability of the proposed model. Additionally, this study identifies distinct differences in slurry diffusion between porous and void media. In porous media, slurry propagation encounters greater hydraulic resistance, leading to rapid pressure attenuation and a limited diffusion range. Conversely, diffusion in void media occurs more smoothly due to the continuous cavity structure, resulting in slower pressure decay and a substantially larger diffusion radius. These findings elucidate the mechanisms governing slurry diffusion under dynamic water conditions and provide a theoretical basis for optimizing grouting parameters and improving construction efficiency in water-bearing strata.

  • Hua Xu, Runfang Sun, Tianbin Li, Qixiang Yan, Yaozhong Cui
    Journal of Rock Mechanics and Geotechnical Engineering. 2026, 18(5): 3832-3851.

    Tunnel portal sections have historically been more susceptible to earthquakes than other components, exhibiting significant seismic damage. However, critical seismic behaviors of portal sections remain unrevealed owing to insufficient consideration of actual topography. Moreover, the extent of asymmetric seismic responses induced by topography remains unclear, which is essential for seismic design. To overcome these limitations, this study replicated the actual geological conditions of a tunnel portal section, including the portal slope, topography, slope and tunnel supports, and the often-overlooked portal wall using large-scale shaking table tests. The asymmetric seismic responses and their impact ranges identified in the experiments were validated through numerical simulations. The results revealed that the seismic damage to the slope is attributable to the presence of the tunnel, with slope acceleration near the tunnel portal increasing by 20 %-40 % compared to slopes without a tunnel. Additionally, the tunnel facilitates seismic wave propagation in specific directions, leading to further seismic damage across the portal section. Portal walls, being exposed structures, are susceptible to higher seismic strain and acceleration than tunnel linings and thus warrant increased attention. Importantly, the asymmetric seismic response was found to vary based on different sides and influence ranges. Within 15 m of the portal, the tunnel was dominated by the open-side asymmetric response of acceleration, strain, and displacement. From 15 to 35 m range, the seismic response of the mountainside was more pronounced, exhibiting increased seismic earth pressure and stress on the right sidewall. Seismic earth pressure and stress diminished within the 35-75 m range and steadily decreased beyond 75 m. This enhanced understanding of seismic behaviors facilitates the targeted establishment of future seismic fortifications based on these classified ranges.