[an error occurred while processing this directive]
Citation Search Quick Search
About
   » About Journal
   » Editorial Board
   » SCI IF
   » Staff
   » Contact
Authors
   » Submit an Article
   » Manuscript Tracking
   » Call for Papers
   » Scope
   » Instruction for Authors
   » Copyright Agreement
   » Templates
   » Author FAQs
   » PACS
Referees
   » Review Policy
   » Referee Login
   » Referee FAQs
   » Editor in Chief Login
   » Editor Login
   » Office Login
Browse Journal
   » Accepts
   » In Press
   » Current Issue
   » Archive
   » View by Fields
   » Top Downloaded
   » Sci Top Cited
Announcement
Links
   »
 
  --2026, 45 (S1)   Published: 01 September 2026
Select | Export to EndNote
Artiles

 2026, 45 (S1): 0-0
[PDF] 358 KB (160)    
Show Abstract

Comparative study on dynamic response and energy evolution of tuff and limestone

YE Haiwang1, 2, 3, LI Rui1, YU Yan4, SHANG Fengchuan5, DONG Yefeng6, SHI Binhong4, LEI Tao1, 2, 3*, WANG Qizhou1, 2, 3, LI Ning1, 2, 3, YU Menghao1, LI Xinping7
 2026, 45 (S1): 1-15 doi: 10.3724/1000-6915.jrme.2024.0992
[PDF] 2301 KB (171)     [HTML]
Show Abstract
In order to compare the dynamic mechanical response and energy evolution laws of tuff and limestone, the uniaxial compressive loading tests with strain rate varying from 40 s-1 to 200 s-1 on the two rocks were carried out using the 50 mm diameter split Hopkinson pressure bar (SHPB) system, and the test phenomena were reasonably interpreted from the microscopic and mesoscopic perspective by combining with X-ray diffraction and nuclear magnetic resonance (NMR) tests. Based on the statistical damage theory of Weibull distribution, the damage differences between tuff and limestone under different strain rates are analyzed. The results show that with the increase of strain rate, the dynamic compressive strength, energy absorption rate and energy consumption density of the two rocks increase, existing a significant strain rate effect. Tuff contains high-hardness quartz and feldspar minerals, and its dynamic compressive strength is higher than that of limestone. In addition, tuff has a complex mineral composition and poor internal pore connectivity, resulting in higher energy absorption rates and energy consumption density. In contrast, limestone is primarily composed of softer calcite minerals and has better pore connectivity, exhibiting greater strain rate sensitivity but lower damage levels.

Shaking table experiment and numerical simulation of acceleration amplification effect in tableland-valley complex terrain across different strata

ZHAO Shixing1, LUO Qirui1*, WANG Wei2, XIONG Feng3, XU Yao2, YANG Jiahao2, ZHANG Min4, LI Yingmin5
 2026, 45 (S1): 16-30 doi: 10.3724/1000-6915.jrme.2025.0349
[PDF] 8465 KB (140)     [HTML]
Show Abstract
 To study the complex seismic response in mountainous areas, a shaking table experimental was designed and completed with a tableland-valley combination terrain as the research object. The dynamic damage characteristics, acceleration and amplification effect distribution of the model specimen was obtained. Further, a finite element model was established to explore the differences in seismic responses between isolated and combined terrains. The main conclusions are as follows: The tableland has less movement constraints than the valley and is more prone to the “whip-lash effect” during the propagation of seismic waves, with its peak horizontal displacement being approximately 20%–50% higher than that of the valley terrain. The distribution of the horizontal ground motion amplification effect on the surface maintains a high spatial similarity with the terrain undulation. The maximum amplification coefficient reaches 3.03 at the top. Different geological structures have a significant impact on the surface seismic response, increasing it by 36.3% (amplification coefficient of 4.13). The valley terrain can reflect and scatter seismic waves, focusing the seismic waves on the surface of the tableland and intensifying the vibration (The horizontal peak acceleration at the top increased by 25.3 %–28.4% compared to that of an isolated tableland). This phenomenon is closely related to the site frequency and the spectral characteristics of the incident wave. The research results reveal the influence of terrain features and geological structures on the seismic wave propagation mechanism, providing an important basis for seismic design of mountainous area engineering.

Mechanical properties of glacier ice and ice-rock composite with different interface dip angles

QIN Changbing1, 2*, WANG Qiangdong1, LI Siyu3
 2026, 45 (S1): 31-46 doi: 10.3724/1000-6915.jrme.2025.0459
[PDF] 5783 KB (103)     [HTML]
Show Abstract
Global warming is driving frequent glacial and ice-rock avalanches on the Tibetan Plateau, posing severe threats to lives and property. To understand the underlying mechanics, this study investigated the mechanical properties of glacier ice and the effect of interface dip angle on ice-rock composites through laboratory tests. The main findings are as follows. (1) The strength of glacier ice demonstrates a significant coupling effect between temperature and confining pressure, decreasing linearly with increasing temperature and increasing linearly with rising confining pressure. (2) The mechanical behavior of ice-rock composites is influenced by the combined effects of temperature, confining pressure, and interface dip angle. As confining pressure increases from 0 to 0.5 MPa, the critical dip angle shifts from 30° to 45°. Moreover, when the temperature reaches T≥-10 ℃, the peak strength of the composites decreases significantly. (3) Under zero confining pressure, the failure mode transitions from ice bulging to composite ice-rock failure, and ultimately to interfacial slip as the dip angle increases. At confining pressures above 0.5 MPa, the failure mode gradually reverts to ice bulging. However, this effect of confining pressure diminishes when temperatures rise above -10 ℃, where interfacial slip becomes the predominant failure mode. This study presents the first systematic triaxial test investigation revealing the significant impact of interface dip angle on the overall mechanical performance of ice-rock composites, providing crucial experimental evidence for understanding the failure mechanics of ice-rock complexes in avalanche disasters.

Orthogonal test and new microscopic simulation method of similar materials of soft surrounding rock in model experiments

CUI Lan1, 2*, ZHUO Hongming1, 2, SHENG Qian1, 2, ZHENG Junjie3, DONG Youkou4, XU Chen5, FAN Leyang1, 2
 2026, 45 (S1): 47-63 doi: 10.3724/1000-6915.jrme.2024.0969
[PDF] 3903 KB (124)     [HTML]
Show Abstract
The results of the current study reveal significant differences in the effects of the proportioning scheme of similar materials used for surrounding rock on their mechanical properties. There is limited discussion on the microscopic simulations of similar materials for surrounding rock. Existing simulations often focus on single-component materials, neglecting the fact that similar materials consist of multiple components. To address this, this study utilizes a test section of a high-speed railway tunnel to prepare similar material for Grade IV surrounding rock in tunnel model tests. The experimental method is designed using an orthogonal test scheme, incorporating barite, fly ash, river sand, rosin alcohol, and machine oil as raw materials. A total of 16 sets of proportioning schemes for the similar materials are developed to conduct uniaxial compression tests and direct shear tests. Based on the results of the orthogonal tests, the gravity, uniaxial compressive strength, elastic modulus, internal friction angle, and cohesion of the 16 sets of similar materials are determined. Among these, the set of samples closest to the target parameters is selected as the similar materials for the surrounding rock in the model test. Subsequently, the influence of control factors on mechanical parameters—namely gravity, uniaxial compressive strength, elastic modulus, internal friction angle, and cohesion—is analyzed. The results from the orthogonal test serve as a crucial reference for the proportioning schemes of similar materials for surrounding rock in model tests. Furthermore, a novel microscopic numerical method is proposed, which considers particle density and radius for similar materials. This method can represent different components of the similar materials, such as barite, fly ash, river sand, and rosin. The relationship between the mechanical parameters of the model (uniaxial compressive strength and elastic modulus) and the microscopic parameters of the model (particle radius, bonding strength, and contact modulus) is established through the numerical method of uniaxial compression. This approach provides an effective numerical framework for analyzing the failure process and mechanical mechanisms of similar materials. By comparing the stress-strain curves, crack evolution patterns, and displacement distributions between the experimental and numerical methods, the deformation behavior and failure mechanisms of the similar materials for surrounding rock are discussed. As a result, it is found that increasing the proportion of barite enhances the gravity and uniaxial compressive strength of the similar materials. Additionally, the cohesion and internal friction angle of the similar materials increase with a higher proportion of rosin alcohol. The proposed new microscopic numerical method effectively simulates the entire loading and failure process of similar materials, aligning well with the experimental results.

A fracturing height propagationevaluation method for multi-lithological laminated shale reservoirs based on the characterization of equivalent macroscopic mechanical properties

LIANG Shuang1, 2, YANG Gongchen1, SUN Shuo3*, WANG Di4, WEI Jianguang2, 3, SHEN Anqi1, 2
 2026, 45 (S1): 64-80 doi: 10.3724/1000-6915.jrme.2025.0218
[PDF] 3443 KB (99)     [HTML]
Show Abstract
To accurately evaluate the fracture height propagation effect in multi-lithology laminated shale reservoirs during hydraulic fracturing, this study proposes an innovative mineral-composition-driven methodology for predicting equivalent macroscopic mechanical properties and enabling multi-parameter synergistic evaluation. First, a representative volume element (RVE) was constructed based on mineral composition. Through calibrated finite element simulations of uniaxial compression experiments, a meso-to-macro scale transition was achieved using periodic boundary conditions, allowing for the prediction and establishment of a quantitative relationship model between equivalent macroscopic mechanical properties and mineral composition. Subsequently, focusing on the shale-limestone interbedded reservoir within the Da?anzhai Formation, an integrated approach incorporating mineral-derived mechanical parameters, in-situ stress parameters, reservoir/interlayer parameters, and construction parameters was developed. The analytic hierarchy process (AHP) was employed to establish a fracture height propagation evaluation index. Parameter sensitivity charts were generated, identifying a critical index value of 0.512 as the penetration threshold for limestone interlayers. Validation was achieved by comparing calculated evaluation indices for Well A?s fracturing stages with field microseismic monitoring results, confirming the method?s accuracy. This research pioneers the integration of mineral mechanics and in-situ stress multi-source parameters for synergistic penetration evaluation, providing an effective tool for assessing and optimizing fracture height propagation in multi-lithology layered shale reservoirs.

Analytical solution of load transfer of inclined partial-cohesive anchor cable during the whole pull-out process

HAN Jianyong1, 2, ZHAO Wen3, CHEN Yang4*, BAI Xiaoyu5, JIA Pengjiao6, LI Tianliang7
 2026, 45 (S1): 81-93 doi: 10.3724/1000-6915.jrme.2025.0219
[PDF] 1401 KB (79)     [HTML]
Show Abstract
Due to the nonlinear characteristics of the pullout behavior of partially cohesive anchor cables, developing a theoretical analytical method for their load transfer mechanism poses certain challenges. This study focuses on the mechanical behavior of the anchorage interface and proposes a shear stress calculation method for the grout-soil interface. Based on the typical pullout failure characteristics of anchor cables, analytical solutions for the load transfer of inclined partially cohesive anchor cables are derived for both the elastic stage and the slippage development stage. On this basis, the distribution characteristics of axial stress and interfacial shear stress of the partial cohesive anchor cable in soil under pullout load are systematically analyzed, revealing the failure evolution throughout the pullout process. The results indicate that the grout in the free section of the partially cohesive anchor cables provides significant resistance to overall pullout performance. The axial stress in the steel strand of the anchorage section follows a hyperbolic sine distribution, while the maximum axial compressive stress of the grout occurs at the junction between the free and anchorage sections, with a tensile stress zone emerging at the bottom of the anchorage section. During the slippage development stage, the length of the interfacial plastic zone at both the free and anchorage sections generally demonstrates a linear relationship. The proposed analytical method effectively captures the shear slip behavior at the anchorage interface of inclined partially cohesive anchor cables and elucidates the corresponding load transfer mechanism, thereby providing theoretical support for understanding their pullout failure evolution and for the design of temporary reinforcement systems in slopes or excavations.

Spatiotemporal evolution of rock burst in surrounding rock with structural planes based on stress relaxation

ZHANG Xiaojun1, 2*, YANG Xiaoyao1
 2026, 45 (S1): 94-105 doi: 10.3724/1000-6915.jrme.2025.0350
[PDF] 1609 KB (85)     [HTML]
Show Abstract

Rockbursts in surrounding rocks containing structural planes present significant challenges to engineering construction. Therefore, the spatiotemporal evolution law is crucial for the prevention and control of deep underground rockbursts. To address the deficiencies in the study of the spatiotemporal evolution of structural plane-induced rockbursts, granite specimens with non-closed structural planes at various angles were prepared. Utilizing the concept of stress relaxation, uniaxial compression stress relaxation tests were conducted to analyze and discuss the evolution of rockburst fractures and the differences between immediate and time-delayed occurrences influenced by structural planes of differing angles. The results indicate that: (1) when the stress relaxation or attenuation process exhibits three stages—rapid attenuation, gradual attenuation, and stable attenuation—it is classified as a stress relaxation process, characterized by the absence of rockbursts or the occurrence of local weak rockbursts with a long time delay. Conversely, if the process shows two stages—gradual attenuation and rapid attenuation—it is classified as a stress attenuation process, characterized by immediate or short time-delayed rockbursts. (2) For a given structural plane, as the load increases, the state of the rock mass containing structural planes transitions from a stress relaxation process (without rockbursts) to a process where stress relaxation transforms into stress attenuation (with rockbursts, characterized by long time-delayed rockbursts), and eventually to a stress attenuation process (with rockbursts, characterized by short time-delayed or immediate rockbursts). (3) As the angle of the structural plane increases, the failure mechanism of rockbursts shifts from being controlled by the structural plane to being governed by strain, with the failure of the structural plane progressing from the tip to the entire structural plane. (4) Under the same peak load, rockbursts are more likely to occur in structural planes with smaller angles, and the time required for these occurrences is shorter. For structural planes with the same angle, higher peak loads result in shorter durations before rockbursts occur. The findings of this research provide a foundational understanding and guidance regarding the mechanisms, prevention, control, monitoring, and prediction of structural plane-induced rockbursts in deep underground engineering.

Capacity expansion characteristics of granite under triaxial compression

ZHU Zeqi1, 2, 3*, TIAN Kiawei1, 2, SHENG Qian1, 2
 2026, 45 (S1): 106-112 doi: 10.3724/1000-6915.jrme.2025.0895
[PDF] 1599 KB (99)     [HTML]
Show Abstract
To study the capacity expansion characteristics of rock, this research uses Dagangshan granite as a case study. Various tests, including porosity tests, triaxial compression tests, acoustic emission tests, and CT scanning, were conducted to investigate the crack propagation and capacity expansion characteristics of granite during triaxial compression. A calculation method for rock pore volume deformation was developed, along with an evolution equation for porosity during triaxial loading. Based on the principle of strain equivalence, the relationship between the total volumetric strain of the rock sample, the volumetric strain of the rock skeleton, and the volumetric strain of the pore was analyzed to derive a relationship expression between rock porosity and bulk modulus. By integrating the porosity evolution equation, a constitutive model that accurately reflects the expansion characteristics of granite was established, which was preliminarily validated through comparison with experimental results.

Wave-tracking-free calculation method for time-domain seismic response of overburden sites

LI Yanpeng1*, HU Zhiqiang2, LI Zhiyuan3, LIN Gao2
 2026, 45 (S1): 113-119 doi: 10.3724/1000-6915.jrme.2025.0382
[PDF] 526 KB (68)     [HTML]
Show Abstract
The time-domain seismic response analysis of sites is fundamental for investigating the true nonlinear seismic behavior of overburden sites. Typically, the time-domain seismic response is derived from the Fourier inverse transformation of frequency-domain solutions. This process involves wave response superposition, which inherently limits frequency-domain methods from accurately capturing the true nonlinear seismic response of the soil. In this study, we propose a recursive method for analyzing the time-domain seismic response of overburden sites without the need for wave tracing, based on the theoretical solution of the wave equation and wave propagation characteristics. This method employs the reflected and transmitted waves incident at each soil interface to construct both up-going and down-going waves at these interfaces, treating each soil interface as a wave source. The up-going and down-going waves at each soil interface for the current moment are generated by the reflected and transmitted waves from incident waves that originated from neighboring soil interfaces at the previous moment. Ultimately, these up-going and down-going waves at the soil interfaces are utilized to determine the free-field time-domain seismic response at any point within the field. The proposed method avoids tracking the reflected and transmitted waves at each soil interface, significantly reducing computational effort and simplifying the wave propagation process. Numerical examples validate the applicability and accuracy of this method, while also exploring the impact of overburden thickness and soft interlayers on seismic wave propagation in overburden sites. The results indicate that the overburden layer significantly amplifies acceleration within a specific depth range from the surface, while exhibiting minimal amplification below this depth. Furthermore, the soft soil interlayer within the overburden demonstrates a notable energy absorption effect on seismic wave energy. This proposed method establishes a foundation for the true nonlinear seismic response analysis of overburden sites.

An empirical model for unsteady creep of rocks

ZHANG Liangliang1, 2*, WU Jianghua1, CHENG Hua1, 3, WANG Xiaojian1
 2026, 45 (S1): 120-131 doi: 10.3724/1000-6915.jrme.2025.0440
[PDF] 563 KB (74)     [HTML]
Show Abstract
In response to the limited types of existing rock empirical creep models and their inability to accurately describe the unsteady creep behaviors of rocks, an improved Knothe time model for predicting dynamic surface subsidence due to coal mining is proposed. By applying a mathematical function rotation transformation, the improved Knothe time model function is rotated counterclockwise by 90°to establish a new empirical creep model for rocks. The accuracy and applicability of this model are validated using conventional triaxial compression creep test data from four types of rocks: marble, mudstone, sandy shale, and sandy mudstone. The results indicate that: (1) the characteristics of the new rock unsteady empirical creep model curve align with the typical characteristics of unsteady creep curves for rocks; the model is simple in form and consists of only three parameters; (2) the curve of the new empirical creep model closely overlaps with the creep test curves of marble, mudstone, and sandy shale, effectively describing both rock attenuation creep and constant velocity creep, as well as the nonlinear variations observed during the acceleration stage; (3) the accuracy of the new rock empirical creep model in fitting the creep test curve of sandy mudstone significantly surpasses that of the power function, logarithmic function, and exponential function empirical creep models, thus overcoming the challenges associated with describing unsteady creep behaviors. This research provides valuable insights and references for the study of rock creep models.

Dynamic mechanical behavior of jointed rock mass under different initial stress states: experimental and numerical analysis

Xie Xiaokun1, YUAN Wei2, Shi Shaoshuai3, WANG Haojun4, 5, Rong Xiaoli1, Wu Weitao1, Hu Jie1*
 2026, 45 (S1): 132-149 doi: 10.3724/1000-6915.jrme.2025.0399
[PDF] 7360 KB (102)     [HTML]
Show Abstract
Prior to blasting excavation, the jointed surrounding rock is subjected to a complex in-situ stress field. During blasting, the combined effect of in-situ stress and dynamic loading may induce dynamic disasters such as collapse. To investigate the dynamic response of jointed rock masses under complex stress conditions, systematic tests were performed on jointed sandstone subjected to confining pressure, biaxial, and true triaxial pre-stresses. Both a conventional confining pressure SHPB system and a true triaxial SHPB device were employed to investigate the evolution of dynamic compressive strength and failure modes. Furthermore, discrete element simulations using UDEC and 3DEC were conducted to investigate crack initiation and propagation paths at different strain rates. These simulations provided deeper insights into the micro-mechanical damage mechanisms of jointed rock mass under multi-axial stress states. Experimental results further indicate that the dynamic compressive strength of jointed sandstone is highly sensitive to these multiaxial stress conditions. Axial prestress () promotes crack initiation and reduces strength, whereas lateral prestress(,) and confining pressure () effectively suppress crack propagation and enhance load-bearing capacity. The most pronounced strength increase is observed under true triaxial confinement. Failure modes vary depending on the stress constraint conditions. Under confining pressure, specimens exhibit conical compression-shear failure. Under biaxial constraint, failure is dominated by tensile fractures along multiple surfaces. Under triaxial constraint, composite failure occurs, including localized spalling at the specimen ends. Discrete element simulations reveal that stress concentration at the intersection of cross joints is the primary mechanism driving crack initiation and propagation. Additionally, the simulations clarify the relationship between microscopic damage evolution and macroscopic failure modes.

Analysis of the upper-bound limit of stability of the tunnel face considering the effect of pipe shed support

SHAN Yepeng1, 2, CAI Guoqing1, 2*, YANG Rui1, 2, LU Yiyun3, QI Zhi3, WANG Fulin3
 2026, 45 (S1): 150-161 doi: 10.3724/1000-6915.jrme.2025.0112
[PDF] 2515 KB (75)     [HTML]
Show Abstract
The pipe shed support demonstrates remarkable effectiveness in enhancing tunnel stability. The conventional geometric parameter recursion method for the elliptical cone failure mode limits computing efficiency.
Therefore, based on ensuring the similarity of projected shapes for all sliding blocks, this study accurately calculated the energy dissipation rate by reasonably considering the cross-sectional area, effective number, and elevation angle of pipe shed crossing the slip surface. The stability safety factor (Fs) was formulated through a combination of the upper-bound limit analysis and the strength reduction method, followed by its solution via nonlinear optimization.
Comparative analysis against model tests, existing theoretical solutions, and engineering cases confirmed the validity of the proposed method and elucidated the influence of the pipe shed elevation angle on Fs. The results reveal that the synergistic interaction between the shear resistance of pipe sheds and frictional forces along slip surfaces reduces the required shear strength of surrounding rock for limit equilibrium, thereby increasing the Fs. The Fs exhibits a non-linear decay trend with increasing elevation angle of the pipe shed, which is more sensitive to such changes. The results of this study contribute to refining the stability evaluation system for tunnel faces under pipe shed support.

Weathering of petroglyph carrier surfaces in Baiyin under coupled wet-dry and chemical effects

KONG Xiangchao1, 2, HE Faguo1, 2*, XIE Mingyu1, 2, GUO Sirui1, 2, LU Tengfei1, 2, MA Yunfeng1, 2,YANG Xingjie1, 2, LV Xuemeng1, 2, ZHANG Lizhi1, 2
 2026, 45 (S1): 162-175 doi: 10.3724/1000-6915.jrme.2025.0599
[PDF] 5721 KB (81)     [HTML]
Show Abstract
The Baiyin petroglyph carrier, influenced by thermodynamic cycles and hydrochemical processes, has developed a composite degradation mode characterized by fractures, bulging, warping, and peeling. Through laboratory-simulated dry-wet cycling experiments, aligned with the actual environmental conditions at the petroglyph site, a quantitative top-dripping method was employed to replicate the precipitation infiltration-lateral migration process that drives the water-rock interface reactions responsible for the chemical weathering of the petroglyphs. Systematic monitoring of the evolution of sample indicators across different cycle periods reveals the differential impact of dry-wet cycles on the weathering process of sandstone surfaces under varying environmental conditions. The results indicate that, under different environmental conditions, the number of dry-wet cycles is significantly positively correlated with the degree of degradation of the sample′s physical properties. The weathering of the carrier is primarily driven by the combined effects of expansion, salt crystallization, and dissolution, leading to changes in both the macroscopic and microscopic structure and mineral content of the sample surface. A weathering damage model is established based on the evolution of the sample′s pore rate, facilitating the quantitative characterization of surface damage under varying environmental stresses.

Characteristics of impulse waves generated by three-dimensional rock-mass landslides in open channel flow

TIAN Ye1, 2, WANG Pingyi2*, YU Tao2, REN Qingyang1, WANG Meili2, 3, HAN Linfeng2
 2026, 45 (S1): 176-187 doi: 10.3724/1000-6915.jrme.2025.0612
[PDF] 4819 KB (92)     [HTML]
Show Abstract
The generation and propagation of landslide-induced impulse waves are strongly influenced by ambient river flow. A three-dimensional physical model of fragmented, jointed rock debris landslides under hydrodynamic conditions was constructed using typical geometries and fracture patterns from the Three Gorges Reservoir area. Changes in waveform and wave-field characteristics were examined for current Froude number from 0 to 0.23. Water-entry behavior was controlled by slide shape and dip angle. Two modes were identified: bulk-sliding mode and two-stage-toppling-sliding mode. The two-stage-toppling-sliding mode produced a secondary major crest and a multi-peaked spectrum. Under a shear current, nonlinear deformation increased with increasing Fr. For waves propagating upstream, crests were strongly amplified, troughs became shallower and flatter, crest duration was lengthened, trough duration was shortened, and the waveform approached a solitary-wave-like form. For waves propagating downstream, crest duration was shortened, trough duration was lengthened, crest height increased slightly, troughs deepened markedly, and crest-trough asymmetry was enhanced. The shear current transformed the wave field from a symmetric near-circular pattern to a downstream shifted asymmetric ellipse. A wave-field deformation coefficient was introduced to quantify this transformation. Attenuation models for co-flowing and counter-flowing waves were developed to delineate upstream and downstream hazard extents. The findings provide guidance for predicting wave propagation and hazard zones in reservoir backwater reaches.

Modeling of permeability and fracture-scale nonlinear response of coal body under multiple stresses

GUAN Shiju1,2, LIU Huihui1, 2*, YU Bin1, 2, LI Tingfang1, 2, WANG Zhiji1, 2, XU Chao3
 2026, 45 (S1): 188-199 doi: 10.3724/1000-6915.jrme.2025.0346
[PDF] 795 KB (65)     [HTML]
Show Abstract
The evolution of coal permeability results from the interplay between external stress, pore pressure, and adsorption stress, with varying responses of the coal body′s fissure scale and permeability control under each stress condition. Therefore, establishing a coal permeability model that accounts for the individual stress factors and the nonlinear response of fissure scale is crucial for accurately and quantitatively characterizing permeability evolution. To this end, this paper first investigates the evolution characteristics of permeability and volumetric strain in the coal body under the influence of single stress factors: external stress, pore pressure, and adsorption stress. It then deduces the dynamic evolution model of the fracture compression coefficient corresponding to each stress and establishes a new model for coal permeability. The results indicate that: (1) Coal permeability exhibits a rapid initial decrease followed by a gradual decline as external stress and methane pressure increase, while under pore pressure influence, coal permeability shows a slow increasing trend. (2) The permeability models of coal mass under effect of external stress, pore pressure and adsorption stress respectively are established. Among them, fracture volume compressibility coefficient of coal mass under external stress follows an exponential decay relationship with external stress, fracture volume compressibility coefficient of coal mass under pore pressure exhibits an exponential growth relationship with pore pressure, and fracture volumetric strain of coal mass under adsorption stress conforms to a Langmuir-like equation with gas pressure. (3) A coal permeability model is developed that comprehensively considers the nonlinear response of fracture scale and individual stress factors. Compared to conventional permeability models, the newly established model demonstrates better alignment with experimental permeability data, providing a more accurate characterization of the influence of individual stress factors on coal permeability.

Fault friction and slip stability experiments under different CO2 phase conditions

LIN Yafeng1, 2, LI Xiaying1, 2, HE Miao1, 2, XIA Yuanyou3, LI Bin1, 2, LI Qi1, 2*
 2026, 45 (S1): 200-212 doi: 10.3724/1000-6915.jrme.2025.0736
[PDF] 2798 KB (76)     [HTML]
Show Abstract
Carbon capture, utilization, and storage (CCUS) is recognized as a crucial technology for achieving significant reductions in CO2 emissions. However, the injection of CO2 into underground reservoirs may induce the reactivation of concealed or nearby faults, thereby posing a threat to storage security. To investigate the characteristics and mechanisms of fault reactivation under different CO2 phases, saw-cut sandstone fractures were utilized to simulate faults. This study examined the effects of various CO2 phases on fault frictional behavior and stability by conducting velocity-stepping shear tests under gaseous CO2, supercritical CO2, and dry conditions using a self-developed L-shaped shear-flow apparatus. The experimental data were analyzed and numerically fitted according to the rate-and-state friction (RSF) law. The results indicate that: (1) at a loading rate of 1 μm/s, all samples exhibited pronounced stick-slip behavior, and the evolution of the friction coefficient during velocity steps aligned with RSF predictions; (2) fracture surface roughness significantly decreased after shearing, with the gaseous CO2 group displaying the most substantial weakening accompanied by normal compression, while the supercritical CO2 and dry groups exhibited normal dilation; (3) the fault stability coefficient remained positive across all three conditions during velocity steps, indicating consistent velocity-strengthening behavior. These findings provide valuable experimental insights into the mechanical stability assessment and risk control associated with CO2 geological storage sites.

An approach for rock density field inversion and spatial reconstruction of mechanical parameters

CHEN Yanzhao1,,2, ZHOU Hao1, 2, YAO Xulong1, 2**, TAO Zhigang2, 3, LUAN Liming4, RONG Hui5, ZHAO Jizhong6
 2026, 45 (S1): 213-225 doi: 10.3724/1000-6915.jrme.2025.0769
[PDF] 2030 KB (105)     [HTML]
Show Abstract
As conventional methods for obtaining mechanical parameters neglect the evolution of internal damage within rock masses, rely heavily on fundamental mechanical experiments, and struggle to accurately reflect the overall mechanical response of rock masses in engineering contexts due to the challenges of obtaining localised mechanical parameters on-site, a density field and mechanical parameter inversion method based on Gardner′s coupledintegrated damage variable is proposed. Precise inversion of rock mass density is achieved by incorporating damage variables into Gardner′s empirical formula, combined with mass conservation constraints and a density range iterative correction mechanism, enabling the further derivation of rock mechanical parameters. Verification was achieved through combined experiments involving sandstone wave velocity imaging and cutting sample density measurements. The average relative error between the inverted density and the measured density was 1.13%, indicating that the two were essentially consistent. Three-dimensional spatial reconstruction of rock mechanical parameter fields is achieved by the inversion methods. The mechanical parameter field exhibits pronounced spatial coupling with the density field. High-density zones correspond to mechanically strengthened zones (featuring peak elastic moduli reaching 20.1 GPa, shear moduli 30.6% above the mean, and bulk moduli 1.9 GPa above the mean), while low-density zones correspond to damage-prone zones (where microcrack development results in damage levels 2.14 times the rock mass mean). The statistical reliability of the reconstructed mechanical parameter field was validated using the Weibull distribution. This approach overcomes the limitations of conventional empirical formulas by employing multi-physics collaborative inversion, thereby providing a quantitative basis for analysing rock mass fracture mechanisms and assessing engineering stability.

P-Wave first arrival picking in microseismic signals using a relaxed continuous minima strategy

CHEN Yifan1, WANG Weifeng1, 2, 3*, DONG Mengyang2, LIU Yixiang2, MA Yansong2, LIU Bin2
 2026, 45 (S1): 226-241 doi: 10.3724/1000-6915.jrme.2025.0624
[PDF] 2292 KB (89)     [HTML]
Show Abstract
The accurate identification of P-wave arrivals in microseismic signals is a critical step for source localization and disaster early warning. In the complex environment of underground coal mines, strong background noise and pulse interference significantly affect the precision of arrival time identification. To address this challenge, a P-wave picking method based on the Akaike Information Criterion (AIC) was proposed, which achieves high-precision identification in the presence of both white noise and pulse interference. First, a composite feature curve (C-AIC) was constructed by incorporating the first- and second-order differences of the signal into the variance calculation of the AIC curve, followed by a weighted fusion of these components. Second, to mitigate mispicking caused by global minimum criteria amidst pulse disturbances, a local-minimum-based detection strategy was developed: consecutive local minima with intervals smaller than a threshold β are grouped into segments, and the first point of the segment containing the largest number of minima is selected as the P-wave arrival. Tests conducted on 80 real microseismic signals demonstrate that the proposed method achieves high-precision P-wave picking under both white noise and pulse interference, with 66 signals (82.5%) exhibiting errors below 5 ms and a mean absolute error of only 2.2 ms. This significantly outperforms VAR-AIC (336.7 ms) and PAI-K (396 ms). Compared to other methods, AIC maintains precision under white noise but is prone to mispicking under pulse interference. In contrast, the proposed method demonstrates enhanced stability and robustness across various noise conditions, indicating substantial practical value for microseismic monitoring and hazard prevention.

An implicit hydro-mechanical coupling scheme using PFC2D with the finite volume method and its validation

ZHOU Jian1, 2*, JIANG Renhua1, PENG Xiao1, ZHANG Yimeng1, ZHANG Lei1, ZHANG Luqing3
 2026, 45 (S1): 242-255 doi: 10.3724/1000-6915.jrme.2025.0143
[PDF] 5252 KB (158)     [HTML]
Show Abstract
Understanding hydro-mechanical interactions in geological media is fundamental for addressing critical challenges in geohazard prediction, resource extraction, and crustal dynamics. While the discrete element method (DEM) coupled with fluid flow solvers has shown promise in these domains, conventional explicit coupling schemes suffer from overly restrictive time-step constraints and excessive computational costs when modeling fluid-fracture interactions. This study presents an innovative implicit coupling framework that synergistically combines the finite volume method (FVM) for solving Stokes flow equations with DEM-based mechanical analysis using two-dimensional particle flow code (PFC2D). The hybrid approach features: (1) An FVM implementation in Python for efficient implicit solution of fluid pressure fields; (2) DEM-based mechanical response calculation; and (3) Tight coupling through iterative information exchange between fluid and solid domains. Three benchmark simulations validate the proposed methodology and its high efficiency: (1) Steady-state seepage analysis demonstrates linear pore pressure distribution along flow paths and Darcy-law compliant flow rates proportional to hydraulic gradients; (2) Transient consolidation modeling of soil pillars reveals excellent agreement with Terzaghis analytical solution across temporal evolution of pore pressure dissipation; (3) Hydraulic fracture propagation simulations capture both viscosity-dominated and toughness-dominated regimes consistent with KGD model predictions, while uniquely resolving fracture-induced pressure oscillations and rock vibration dynamics; (4) The calculation time step size in the three cases is 3–4 orders of magnitude larger than that of the original method, leading to a significant improvement in computational efficiency. The demonstrated capabilities in modeling coupled flow-deformation-fracture processes position this framework as a powerful computational tool for geological engineering studies.

Discrete element simulation of the dynamic compaction for rock-filled embankment considering block stone gradation and shape

XU Ping1, 2, LUO Yanghao1, QIAO Shifan3*, DONG Hui1, 2, CHEN Bin1, 2, HE Jiaxiong1
 2026, 45 (S1): 256-270 doi: 10.3724/1000-6915.jrme.2025.0551
[PDF] 10332 KB (80)     [HTML]
Show Abstract
To analyze the differences in the dynamic compaction process and effects of rock-filled embankments under the conditions of heavy hammer with low drop height and light hammer with high drop height, a discrete element model was developed that accounts for the shape of block stones and particle breakage. By comparing and analyzing the displacement, stress response, and particle breakage characteristics of the rock-filled embankment under the two conditions, the construction parameters and effective reinforcement ranges for dynamic compaction were determined. The results indicate that: (1) The analysis of tamping times, crater depth, and heave deformation characteristics under both conditions shows high consistency with field test measurements, demonstrating that the established discrete element model for dynamic compaction can effectively analyze the reinforcement effect of rock-filled embankments; (2) The critical vertical displacement area and horizontal displacement area of the block stone fillers, caused by the impact energy of dynamic compaction, exhibit pear-shaped and apple-shaped profiles, respectively. The effective vertical reinforcement depths for dynamic compaction using heavy hammer with low drop height and light hammer with high drop height are approximately 4.9 m and 4.5 m, respectively, while the effective lateral reinforcement widths are approximately 3.9 m and 3.7 m, respectively; (3) The impact energy from dynamic compaction leads to significant particle fragmentation within the block stone fillers, with the fragmentation development process occurring as follows: initial state→local fragmentation→through fragmentation→overall fragmentation; (4) The reinforcement effect of dynamic compaction using heavy hammer with low drop height is superior to that of light hammer with high drop height. The heavy hammer with low drop height not only effectively reinforces the block stone fillers in the shallow areas of the embankment but also significantly improves the compaction effect in the deeper areas. These research findings provide a basis for the design of construction parameters and the evaluation of reinforcement effects in dynamic compaction projects involving rock-filled embankments in mountainous areas.

Prediction of elastic modulus in heterogeneous rocks through fusion of multi-source images and artificial intelligence

PENG Lei1, 2,LI Mingyao1, 3*,LIU Dejun1, 3,ZUO Jianping1, 3
 2026, 45 (S1): 271-281 doi: 10.3724/1000-6915.jrme.2025.0627
[PDF] 3194 KB (93)     [HTML]
Show Abstract

Accurate measurement of the elastic modulus of rock is a key component in geotechnical engineering construction design and stability evaluation. However, in deep high-stress environments and brittle formations with well-developed weak structures, obtaining intact and standardized rock samples for direct measurement of elastic modulus is often challenging. This study proposes a method for predicting elastic modulus by leveraging artificial intelligence to build a bridge between microstructural characteristics and elastic modulus, requiring only microstructure images as input to efficiently infer the modulus. The implementation process consists of two stages: first, an improved DeepLabv3+ model is employed to perform semantic segmentation experiments on a multi-source image dataset, generating microstructural characterization images; second, a RockNet-E model is constructed to establish a mapping relationship between the characterization images and the elastic modulus, enabling accurate prediction. Research results show that the improved DeepLabv3+ model achieves mean intersection over union (IoU) of 87.7%, 97.5% and 87.9% for quartz, feldspar, and mica, respectively, outperforming mainstream models such as DeepLabv3+, PSPNet, and UNet. On the validation set, the RockNet-E model achieves mean absolute error (MAE), root mean square error (RMSE), mean absolute percentage error (MAPE) and R-squared (R2) values of 1.20 GPa, 1.51 GPa, 1.83% and 0.61, respectively. The elastic modulus prediction errors are all below 5%, and the predicted value distribution matches the experimental results. This study provides a low-cost and highly efficient approach for intelligent, image-driven prediction of rock mechanical parameters.

Rock mass structural characterization based on machine vision and grout take prediction

OUYANG Shaoming1, DING Changdong1, 2*, ZHOU Liming1, DING Xiang3, LUO Rong1, HU Dawei2, 4
 2026, 45 (S1): 282-296 doi: 10.3724/1000-6915.jrme.2025.0531
[PDF] 6167 KB (89)     [HTML]
Show Abstract
Aiming at “black-box” technical bottleneck characterized by the inherent concealment and evaluation challenges in underground cavern grouting, this study proposes an integrated framework combining deep learning and machine learning for refined characterization of rock mass structures and quantitative assessment of grouting performance. The results indicated: (1) By employing high-resolution image acquisition and 3D point-cloud reconstruction technologies to construct a dataset of the surrounding rock, an efficient fracture identification and segmentation method based on the CBAM-UKAN model was developed. The proposed method achieved IoU and Dice coefficients of 0.76 and 0.86, respectively, on the validation set. (2) Through the integration of morphological closing, skeleton extraction, and a multi-scale geometric quantification algorithm chain, key fracture parameters such as length, width, and area were extracted, providing essential inputs for the rock mass structural characterization. (3) Comparative analysis of four ensemble learning models identified the simulated annealing-optimized Gradient Boosting Regression Tree (SA-GBRT) as the best-performing predictive model, achieving an R² value of 0.79 on the test set. The proposed approach enables reliable prediction of unit ash consumption based on rock mass fracture characteristics and grouting construction parameters, thereby providing a quantitative basis for pre-evaluation of seepage-control performance. These findings offer significant practical value for the intelligent optimization and scientific decision-making of grouting operations in underground engineering.

Influence of tight sandstone sliding friction characteristics on hydraulic fracture propagation mechanism:A case study from the Chang 7 Member, Ordos Basin

WU Tao1, YUAN Xu2, 3, 4, ZHANG Feifei2, 3, 4*
 2026, 45 (S1): 297-308 doi: 10.3724/1000-6915.jrme.2025.0249
[PDF] 3244 KB (68)     [HTML]
Show Abstract
The frictional properties of rocks exert a critical control on the propagation pathways and morphological evolution of hydraulic fractures during fracturing stimulation. In response to the lack of quantitative characterization of frictional properties in the fracture propagation mechanisms of tight sandstone, this study takes the Chang-7 tight sandstone in the Ordos Basin as the research subject. By integrating laboratory friction experiments with numerical simulation of fracture propagation, it investigates the influence mechanisms of rock surface roughness, mineral composition, and fluid wetting conditions on frictional properties. Furthermore, it quantitatively analyzes the governing effects of the friction coefficient and injection rate on fracture propagation paths and complexity. The results indicate that: In fluid-saturated environments, clay minerals undergo hydration and expansion, which enhances interfacial lubrication and significantly reduces shear slip resistance. Hard mineral particles form localized strong contacts at the interface, prone to inducing stress concentrations and leading to pronounced fluctuations in the friction coefficient. The content of clay minerals is a key factor determining the frictional regime of tight sandstone. By controlling the content of clay inhibitors in slick water, interfacial frictional stability can be effectively maintained, facilitating a transition in frictional behavior from stick-slip to stable sliding. As the friction coefficient increases, the ability of hydraulic fractures to intersect natural fractures first strengthens and then weakens. Maintaining the friction coefficient around 0.4 yields the greatest fracture extension length and complexity. This research contributes to revealing the formation mechanisms of fracture networks under the combined effects of fracturing fluid, hydraulic fractures, and natural fractures in tight sandstone reservoirs.

Prediction of rainfall induced landslides and debris flows based on hydrological models

LUO Jie, WANG Hairong, KOU Boxiao, TANG Zhihui*
 2026, 45 (S1): 309-329 doi: 10.3724/1000-6915.jrme.2025.0296
[PDF] 5255 KB (125)     [HTML]
Show Abstract
Rainfall plays a crucial role in triggering these disasters. Traditional studies utilizing statistical analysis and machine learning to forecast these disasters primarily rely on historical data, which results in models with limited generalizability. To enhance the simulation of rainfall-induced disaster processes, this study focuses on Muping Town and develops an integrated hydrological model based on the Macro Dynamic Propagation Model (MDPM) to simulate rainfall-runoff evolution in the region. By incorporating raster data, the study area is transformed into a network structure, abstracting the rainfall propagation process as a network fluid flow problem. Additionally, based on the mechanisms underlying rainfall-induced disasters, the model is divided into five computational modules: surface runoff propagation, groundwater propagation, vegetation interception, infiltration, and runoff erosion. The rainfall simulation results are represented using water argument and groundwater argument data. The study designed four rainfall scenarios for simulation: single-peak rainfall (peak intensity of 100 mm/h), single-peak rainfall (peak intensity of 200 mm/h), constant rainfall (50 mm/h), and constant rainfall (100 mm/h). The model's simulation results are visualized and compared with historical disaster data from the study area. The findings indicate that the model effectively captures the formation processes of natural phenomena, such as groundwater accumulation and surface runoff generation. It vividly illustrates the temporal and spatial evolution of groundwater and surface runoff during rainfall, successfully identifying high-risk zones susceptible to landslides and debris flow disasters within the study area. This has practical significance for the prevention and control of geological disasters.

Excess interstitial water pressure effect in bedding-jointed rock mass slopes during strong earthquakes

LIN Qinghua, JU Nengpan*, XIE Mingli, HE Yuhao, WANG Zhaoyu, FENG Changyun, ZHOU Shuyu, SHU Zihao
 2026, 45 (S1): 330-342 doi: 边坡工程;强震区;岩质滑坡;间隙水压力;滑坡启动机制
[PDF] 3023 KB (69)     [HTML]
Show Abstract
Wenchuan Earthquake: Water-saturated, interlayer-fractured jointed rocks triggered high-speed landslides during intense seismic activity, driven by excess interstitial water pressure. To elucidate the initiation mechanism of water-saturated, steeply dipping bedding slopes under the coupling of vibration and interstitial water pressure, this study employed shaking table and PFC-FLAC tests. Results: (1) Under 0.8 g seismic waves, the peak ground acceleration (PGA) amplification of the model’s weak layer reached 1.2–1.27, with interstitial water pressure surging by 100 Pa (excess triggered). The model failed under 1.8 g, 20 Hz sine vibration (failure threshold). (2) Vibration caused concentrated excess interstitial water pressure within the weak layer, leading to incompatible deformation due to differences in rock properties. Interstitial water pressure was transiently generated and gradually accumulated; the phenomenon of “vibration compression-tension alternation” (induced by strong quakes) resulted in transient excess pressure. (3) Numerical tests summarized the laws of fracture propagation and deformation, confirming the interaction between water and rock (expansion and excess effects) during strong seismic events. (4) The mechanism of landslides in layered jointed rocks is described as follows: “strong quake → incompatible deformation-induced damage and dilatancy in the weak layer → alternating expansion and excess effects → triggering and accumulation of excess pressure → high-speed landslide.”

Spatial response of blast vibration in deep dual-track tunnel

WANG Changjian1, 2,YU Chong1, 2*,MA Yongan1, 2,LI Haibo1, 2,LIU Yaqun1, 2
 2026, 45 (S1): 343-353 doi: 10.3724/1000-6915.jrme.2025.0570
[PDF] 6020 KB (74)     [HTML]
Show Abstract
During the blasting excavation of deep-buried tunnels, the dynamic disturbance caused by the blasting load and the transient unloading is the key factor for the instability of the surrounding rock of the tunnel. This study takes a double-track tunnel in the western region as the background and uses the finite difference method to explore the spatial distribution characteristics of blasting vibration velocity on the tunnel contour surface, the superposition effect of ground stress and dynamic disturbance of blasting load, and the influence of stress adjustment under plastic deformation on blasting vibration velocity. The research results show that: (1) During blasting excavation, the blasting vibration velocity on the contour surface of the adjacent tunnel shows significant spatial differences. The peak particle vibration velocity (PPV) on the blast facing side is significantly greater than that on the back side. The PPV is the largest at the middle side wall on the blast facing side and the smallest at the arch top. (2) The initial in-situ stress level is positively correlated with PPV, and the vibration velocity induced by transient unloading under high in-situ stress increases significantly. (3) While using the elastoplastic constitutive model, the surrounding rock undergoes plastic deformation. After stress adjustment, the stress at the excavation face and the pilot tunnel is less than the initial ground stress. The spatial distribution characteristics of the blasting vibration velocity at key positions are consistent with the elastic constitutive model. The PPV changes with the ground stress in a different way. The increase in vibration velocity caused by transient unloading is lower than that of the elastic constitutive model.

Influence of roughness on shear behavior of red sandstone joint infilling with red clay

HU Yunjin1, 2, PAN Jianbo1, 2, PENG Jun1, 2 *, WANG Hanglong3, YUE Enxi1, 2, ZHU Jiahui1, 2
 2026, 45 (S1): 354-364 doi: 10.3724/1000-6915.jrme.2025.0588
[PDF] 3005 KB (92)     [HTML]
Show Abstract
The shear behavior of naturally-infilled rock joints is critical to the stability evaluation of engineering rock masses. This study investigates the effect of roughness on the shear behavior of infilled red-sandstone joints. Three Barton standard profiles are fabricated on natural red sandstone by a CNC engraving technique, and a series of direct shear tests are performed under normal stresses ranging from 0.50 to 1.50 MPa to examine the strength and failure mechanism. The results show that: (1) the shear stress–displacement curves can be divided into an elastic stage, a plastic (hardening) stage, and a slip stage. Increase of roughness transforms the stress-displacement response from a non-peak type to a clear peak-type with a more pronounced stress drop, indicating a brittle behavior at failure. (2) Increasing normal stress leads to a marked increase of the shear strength of infilled joints. When the normal stress rises from 0.50 to 1.50 MPa, the strength increases by 271%, 298%, and 191% for low, medium, and high roughness joints, respectively. The shear strength of infilled joints does not invariably increase with increasing roughness and depends on both joint morphology and asperity characteristic. (3) The presence of infilling material significantly modifies the shear mechanism, and the failure mode gradually shifts from sliding to asperity shearing as roughness and normal stress increase. (4) Binary image–based damage analysis of the joint surface indicates that the damage coefficient of the high-roughness joint is 3.459% under normal stress of 1.50 MPa, which is markedly higher than those of the low- and medium-roughness joints (0.507% and 0.765%). Damage and abrasion zones concentrate near asperity summits. Increasing normal stress induces stress concentration, enlarges the abrasion area, and significantly intensifies surface damage.

Permeability reduction and microscopic mechanism of sandy gravel soil induced by wollastonite carbon sequestration

LI Junhong, CHEN Zhimin*, ZHANG Runlong, ZHANG Kaikai, LU Shuai, SHI Haobo,TAN Chenglong
 2026, 45 (S1): 365-376 doi: 10.3724/1000-6915.jrme.2025.0634
[PDF] 3387 KB (73)     [HTML]
Show Abstract
Aiming at the problem of tunnel water leakage and in response to the national "dual-carbon" initiative, this study proposes to reinforce sandy gravel soil using CO2-induced silicate minerals to reduce its permeability, thereby achieving the goal of preventing water leakage. Through investigations, mineral particle size, dosage, reaction pH value, and carbonation time were selected as influencing factors. Orthogonal tests were conducted to analyze the law governing how each factor affects the permeability of sandy gravel soil, and further exploration was made into the significance and variation law of the impact of interactions between factors on permeability. In addition, the influence of different catalysts on the reaction conversion rate was studied, and microscopic tests such as X-ray diffraction (XRD) and scanning electron microscopy (SEM) were integrated to reveal the evolution law of the internal microstructure of the samples and the mechanism of the carbonation reaction.The test results show that carbon sequestration by wollastonite can significantly reduce the permeability of sandy gravel soil, with the permeability coefficient being reduced by up to three orders of magnitude. the permeability coefficient tends to increase gradually as the mineral particle size decreases; it decreases with the increase of mineral dosage and carbonation time, and increases as the reaction pH value rises. The optimal carbonation mode is determined as follows: mineral particle size of 75 μm, 30% dosage, pH = 4, and carbonation duration of 24 hours. Among all factors, mineral dosage has the most significant impact on the permeability coefficient, and the interaction between mineral particle size and mineral dosage exerts a significant influence on the permeability coefficient. Acetic acid demonstrates the best catalytic effect, with a conversion rate ranging from 25.08% to 29.04%, followed by sodium bicarbonate catalyst. Carbonation promotes the formation of carbonate crystals (CaCO3), which can effectively fill the internal pores of the samples and bond soil particles. However, most of these crystals adhere to the surface of the soil, and the number of crystals that play a bonding role on soil particles is relatively small—this prevents the formation of a cross-network skeleton. Therefore, the reduction range of the permeability coefficient of sandy gravel soil remains limited.

Rebound and recompression deformation characteristics of soft soil foundation in levee back-cutting

LIU Jie1, LIU Yong2, LU Zheng1*, LI Zhiyong3, YAO Hailin1
 2026, 45 (S1): 377-388 doi: 10.3724/1000-6915.jrme.2025.0148
[PDF] 1991 KB (46)     [HTML]
Show Abstract
The results indicate that the rebound deformation rate of the soil and the deformation during reloading are primarily associated with the unloading ratio and the reloading ratio, respectively. The rebound deformation of soft soil surrounding the pipeline, induced by the reverse trenching construction of existing dike conveyance pipelines, exhibits a “reverse pot-bottom shape” characterized by larger deformation in the middle and smaller deformation toward the sides. A higher foundation compression modulus results in reduced pipeline deformation and vertical displacement of the surrounding soil. As pipeline stiffness increases, the ability of the pipeline and surrounding soil to deform in coordination deteriorates, leading to a corresponding increase in the maximum bending moment experienced by the pipeline. For pipelines with joints, the vertical displacement demonstrates a pronounced slope mutation at the joint locations, and the pattern of vertical displacement transitions from the “reverse bell shape” of jointless pipelines to a “funnel shape”. The maximum displacement of pipelines with joints surpasses that of jointless pipelines, and the vertical displacement of pipelines increases linearly with the number of joints, exhibiting a slope discontinuity at the maximum displacement. Increased saturation of the foundation soft soil causes the maximum rebound and recompression deformation of the foundation soft soil to exhibit a linear or exponential increasing trend. The proposed empirical relationship between the rebound modulus and the unloading ratio can effectively and accurately predict the rebound deformation of soft soil foundations in reverse trenching, providing theoretical guidance for the design, construction, and maintenance of buried pipelines in dike engineering.

Analysis of shear strength indices of rooted soil based on Bayesian theory

FU Jiangtao1, 2, LI Huatan3, 4*, GUO Hong1, 2, LI Xiaokang5, WU Jie6, LI Shuang6, GAO Yuan6,ZHANG Yunru1, 2, ZHENG Haiyun1, 2
 2026, 45 (S1): 389-401 doi: 10.3724/1000-6915.jrme.2025.0616
[PDF] 1666 KB (45)     [HTML]
Show Abstract
To investigate the variability and correlation of shear strength indices (cohesion (c) and internal friction angle (?)) of rooted soil in alpine grasslands, this study selected alpine grasslands with varying degrees of degradation in a county in Qinghai as the testing site. The shear strength indices of the rooted soil samples collected from this site were determined through indoor direct shear testing. A joint probability distribution model, based on linear regression and Bayesian analysis, was constructed to describe the variation of cohesion with the internal friction angle. The model parameters were estimated using the Markov Chain Monte Carlo (MCMC) method. The results indicated a significant negative linear correlation between cohesion and the internal friction angle of the rooted soil. Additionally, the samples generated by the joint distribution model grounded in Bayesian theory demonstrated a high consistency with the measured data, with the relative deviation of the mean internal friction angle being only 0.34%. It was found that the transformation uncertainty error  significantly affected parameter correlation, with the correlation coefficient decreasing from -0.95 to -0.43 as   increased from 0.5 kPa to 3.0 kPa. This study confirmed that the Bayesian method effectively quantifies the uncertainty of shear strength indices of rooted soil, providing a scientific basis for evaluating soil anti-erosion capacity in grassland areas and for formulating ecological restoration measures.

Field test for behavior of pre-bored PHC pile in deep soft soil layers

ZHOU Jiajin1, 2, ZHANG Rihong1, 2*
 2026, 45 (S1): 402-410 doi: 10.3724/1000-6915.jrme.2025.0613
[PDF] 916 KB (44)     [HTML]
Show Abstract
In this paper, the bearing performance of pre-bored PHC piles embedded in deep soft soils was analyzed through two groups of field tests. The diameters of the six test piles were 0.75 m and 0.9 m, respectively, and distributed fiber optic sensors were installed along the pile shafts to measure the variation in axial force during the loading process. The test results indicated that when the 0.75 m pre-bored PHC piles were loaded to the designed ultimate capacity of 3 520 kN, the pile head displacements ranged from 2.48% to 2.56% of the pile diameter (D). Similarly, the pile head displacement of the 0.9 m pre-bored PHC piles, when loaded to the ultimate capacity of 5280 kN, also ranged from 2.48% to 2.56% of D. The applied load at the pile head was predominantly supported by shaft resistance when the piles were loaded to their designed ultimate capacities in this study. The grouting stage during the installation of pre-bored PHC piles significantly enhanced the frictional capacity of the pile-soil interface, with the measured ultimate skin friction of the pre-bored PHC piles ranging from 1.04 to 1.18 times that of bored piles across different soil layers. Furthermore, the CPT method was utilized to calculate the ultimate skin friction of pre-bored PHC piles in various soil layers, and the relationship between the ultimate skin friction and cone tip resistance (qc) was established based on the field test results.

Intelligent subsidence prediction of super large anchor caisson foundations during construction based on multi-source data fusion

XU An1, WEI Xiangping1, JIANG Fan1, YAN Hezhong1, DONG Xuechao2, 3*, GUO Mingwei2, 3
 2026, 45 (S1): 411-422 doi: 10.3724/1000-6915.jrme.2025.0253
[PDF] 1485 KB (52)     [HTML]
Show Abstract
Abnormal subsidence of open caisson foundations may induce risks such as delayed sinking and sudden sinking during construction, and thus real-time and accurate prediction of subsidence is essential to ensure construction safety. An intelligent subsidence prediction model is proposed for open caisson sinking based on multi-source data fusion and the Extra Trees Regressor (ETR) algorithm. By integrating multi-source monitoring data, including soil pressure and structural stress, the model extracts nonlinear relationships between data features and subsidence to achieve dynamic prediction. Based on the construction of the super-large open caisson foundation for the north anchorage of the Zhangjingao Yangtze River Bridge’s North Channel Bridge, two single-source prediction models (based on soil pressure on the cutting edges or structural stress) and two multi-source fusion models (feature fusion and decision fusion). The results demonstrate that the multi-source decision fusion model achieves optimal performance, with a root mean square error of 0.003 85 and a coefficient of determination of 0.88. The model requires an average computation time of 1.98 s per prediction cycle, and the average relative error of cumulative subsidence prediction is only 0.19%, which meets both accuracy and timeliness requirements for engineering applications. Sensitivity analysis reveals that the model exhibits high stability and robustness, as its performance remains insensitive to variations in the number of base decision trees. Furthermore, a graded early-warning mechanism is proposed, integrating prediction results with engineering thresholds to establish a closed-loop “monitoring-warning-feedback” management system, providing reliable technical support for risk control during open caisson construction.

The shape effects of prefabricated vertical drains and smear zones on consolidation

HUANG Chaoxuan1*, LEI Guohui2, WANG Jun1, GU Tingwei1, DENG Yuebao3
 2026, 45 (S1): 423-436 doi: 10.3724/1000-6915.jrme.2025.0329
[PDF] 1999 KB (50)     [HTML]
Show Abstract
The smear zone created by the installation of a prefabricated vertical drain (PVD) is typically considered elliptical in practice, whereas current theoretical methods assume it to be circular. This assumption impacts the accuracy of consolidation calculations. Utilizing the consolidation theory of equal strain for a vertical drain foundation and the assumption of an elliptical cylinder for both the PVD and smear zones, analytical solutions for excess pore water pressure and the average degree of consolidation under complex loads are derived. An example analysis is conducted using finite element numerical simulations and laboratory test data to validate this method. The results indicate that the ratio of the equivalent diameter of the elliptical smear zone to the equivalent diameter of the area equivalence decreases as the oblateness of the elliptical smear zone increases. When the elliptical oblateness is substantial, the ratio gradually diminishes. Conversely, the equivalent diameter ratio increases with an increase in the width of the PVD. Under conditions of equal smear zone areas, a flatter smear zone results in a smaller adverse effect and a faster consolidation rate of the foundation soil. A fitting formula for the equivalent diameter of the oblate elliptical smear zone is provided for commonly used prefabricated vertical drains in engineering applications, offering a technical reference for related engineering design and research.

Thermal conductivity of lunar regolith simulants with varying water-ice content

GUO Liang1, PAN Pengzhi1, 2, DONG Yi1, 2*
 2026, 45 (S1): 437-448 doi: 10.3724/1000-6915.jrme.2025.0956
[PDF] 2392 KB (51)     [HTML]
Show Abstract
To meet the demands for water ice detection and resource utilization in the lunar polar regions, this paper systematically investigates the thermal conductivity of simulated lunar regolith (IRSM-L1) under vacuum conditions. Utilizing the transient plane source method, the effects of dry density, vacuum pressure, temperature, and water ice saturation on thermal conductivity were analyzed. The results indicate that at vacuum pressures below 10?² Pa, dry density does not significantly influence thermal conductivity. Conversely, within the pressure range of 1 to 10? Pa, thermal conductivity markedly increases with rising pressure, primarily due to enhanced heat transfer through molecular collisions of gas. The temperatures impact on thermal conductivity is modulated by gas pressure: under extremely high vacuum conditions, thermal conductivity increases with temperature. However, at higher vacuum pressures, gas conduction becomes the dominant mechanism, with thermal conductivity primarily governed by pressure, increasing as vacuum pressure rises. By comparing samples prepared using the wet-mixing method and the ice-powder method, the critical role of water ice morphology on thermal conductivity was revealed. In the ice-powder method, ice particles are distributed in point contacts, resulting in a gradual linear increase in thermal conductivity with saturation, dominated by a series conduction mechanism. In contrast, the wet-mixing method exhibits a three-stage nonlinear growth pattern, characterized by a sharp rise in thermal conductivity within the 20% to 80% saturation range due to the formation of continuous ice-bridge networks, where the conduction mechanism aligns more closely with the geometric mean model. These findings provide important experimental evidence for analyzing the thermal properties of water ice-bearing lunar regolith in polar regions and for identifying the occurrence state of water ice.

Particle size effect investigation on soil arching in sands

HAO Shuai1, LUO Weiping2, ZHENG Zaidi3, XUE Dongjie4, LIU Xueyan1*
 2026, 45 (S1): 449-464 doi: 10.3724/1000-6915.jrme.2025.0445
[PDF] 6253 KB (62)     [HTML]
Show Abstract
The soil arching effect is a crucial factor to consider in tunnel excavation, bridge pile foundations, slope cutting, and other construction projects. It is influenced by factors such as trapdoor width, soil cover depth, and soil strength. However, particle size also plays a significant role in the soil arching effect and requires further investigation. This study employs cubic and hexagonal particle packing models and utilizes the discrete element method to explore the evolution of the soil arching effect under plane strain conditions for single-sized sand in both loose and dense states. The study analyses the regulatory mechanisms of particle size on the ground reaction curve, shear band shape, loosening zone range and shape, and the contact force network.
Results suggest that: (1) under the same conditions, the maximum normalized arch stress, ultimate normalized stress, and maximum normalized arch displacement decrease with increasing particle size; (2) the shear band shape varies with particle size: large and medium-sized particles (d = 3–6 mm) induce tower-like shear failure above the trapdoor, while smaller particles (d = 1 mm) exhibit more complex failure due to the interaction of horizontal and triangular shear bands; (3) regarding the loosening zone, dense sand with large particles (d = 6 mm) displays only a triangular loosening zone, whereas medium and small particles (d = 1–3 mm) show both triangular loosening and inward-contracting funnel-shaped settlement. In contrast, loose sand tends to fail over a larger area, exhibiting a contracting funnel-shaped loosening pattern; (4) as particle size decreases, the average contact force between particles diminishes, while contact density increases. The reduction in contact force weakens the efficiency of load transmission between particles, which is not conducive to the formation of the soil arch, whereas the increase in contact density enhances the cohesiveness of the soil, allowing more particles to contribute to the evolution of the soil arch.

Calculation and application of three-dimensional foundation bearing capacity for large-scale open caissons

LI Jiahang1, 2, GUO Mingwei1, 2*, LU Zheng1, 2, DONG Xuechao1, 2, ZHENG Qinggang3, JIANG Fan4,PENG Wei3, YIN Dongya5
 2026, 45 (S1): 465-475 doi: 10.3724/1000-6915.jrme.2025.0676
[PDF] 1740 KB (54)     [HTML]
Show Abstract
The calculation formulas for foundation bearing capacity recommended by the current standards are typically derived based on the plane strain conditions of strip footings or obtained from small-scale load tests. However, there are still certain limitations in the calculation of foundation bearing capacity for super-large open caissons. Based on the spatial elastic stress state of a semi-infinite soil mass under a vertical circular load, a simplified three-dimensional bearing capacity formula is derived. For non-circular foundations, conversion is performed based on the equivalent area of the foundation base. For rectangular foundation types with different aspect ratios (), the finite element analysis method is used to reveal the variation law of the three-dimensional bearing capacity of the foundation. The calculation errors of the bearing capacity of rectangular foundations with different aspect ratios were clarified. The results show that, for a constant base area, the calculation error of the bearing capacity of the simplified formula increases with the increase of the aspect ratio of the rectangular foundation. When , the calculation error is less than 5%. When , the error is less than 16%. When , the error approaches or even exceeds 20%. The calculation results of the large caisson foundation project of Changtai Yangtze River Bridge show that the allowable bearing capacity of the silty clay layer at the bottom of the open caisson obtained by the simplified formula is 2.93 MPa, which is about 23.1% higher than the value calculated by the railway code. Moreover, field settlement monitoring of the open caisson indicates that the foundation is generally in the elastic deformation stage, which to some extent verifies the rationality of the results. This study clarifies the calculation errors of the three-dimensional bearing capacity formula for rectangular foundations with various aspect ratios. The findings enrich the theoretical calculation of bearing capacity for large-sized deep foundations, and provide an important reference for the design of similar large rectangular, rounded rectangular and rounded-end open-caisson engineering.

Analytical solution for seismic response of large diameter, double-layer concrete-lined shafts in unsaturated soil under vertically-incident P-waves

YANG Jiale1,DI Honggui1*,YAO Yanming2,ZHOU Shunhua1,ZHANG Xiaohui1,HE Chao1
 2026, 45 (S1): 476-489 doi: 10.3724/1000-6915.jrme.2025.0124
[PDF] 1307 KB (41)     [HTML]
Show Abstract
Based on the theory of three-dimensional continuum mechanics, an analytical model was developed to investigate the seismic response of large-diameter, double-layer concrete-lined shafts in unsaturated soil subjected to vertically incident P-waves. The three-dimensional dynamic governing equations for unsaturated soils were established, taking into account capillary pressure and a saturation-dependent dynamic shear modulus. The double-layer shaft lining was idealized as two parallel Rayleigh-Love rods to account for their transverse inertial effects.
Initially, the free-field response of the unsaturated soil was determined using one-dimensional wave propagation theory. Subsequently, the wave equations of the soil were solved through the potential function method, operator decomposition, and the separation of variables. A closed-form solution for the seismic response of the shaft was derived by applying the boundary conditions of the coupled soil-shaft system. The model proposed in this study was validated through comparison with existing results and followed by a parametric sensitivity analysis to evaluate the influence of soil saturation, shaft geometry, and material properties on the seismic response. The results indicated that the seismic response of large-diameter, double-layered shafts is significantly influenced by the excitation frequency: the low-frequency behavior is governed by the static characteristics of the shaft, while the high-frequency response is dominated by the soil-structure dynamic interaction. As soil saturation increases, both the natural frequency of the system and the kinematic amplification factor at the shaft head decrease. Additionally, the transverse inertial effect of the double-layer lining has a negligible impact within the primary frequency range of vertical seismic loading. The proposed analytical framework provides a theoretical basis for the seismic design and safety assessment of shafts in unsaturated ground.

Estimating method of retaining pile and surface deformation of loess pile-supported foundation pit#br#

WANG Xiaoyong1, 2*, ZHENG Jianguo1, 3, YU Yongtang2, MAO Jidong4, ZHANG Zhiyi2, YANG Lei2, LIANG Yi2
 2026, 45 (S1): 490-503 doi: 10.3724/1000-6915.jrme.2025.0333
[PDF] 5114 KB (59)     [HTML]
Show Abstract
This study addresses the critical challenge of accurately predicting retaining pile deformation and ground surface settlement of pile-braced foundation pits located in loess areas. Based on the data collected from 19 projects in Xi?an, the horizontal deformation of retaining piles and the ground surface settlement were systematically analyzed. By comparing the deformation characteristics of foundation pits in loess and soft clay, a modified Gaussian curve and a skewed distribution curve were proposed to better characterize the settlement profiles. By systematically integrating the effects of pile insertion ratio, the piles? flexural rigidity, and the bracing stiffness, a comprehensive stiffness coefficient  for the entire support structure was established. The findings reveal that the pile deformation exhibits a concave profile, with the maximum displacement located at the depth of 0.55–0.75 times the depth of the foundation pit (H). Horizontal deformation diminishes and approaches zero beyond a depth of 1.5H. The maximum ground surface settlement is observed at a distance of 0.5H from the pit?s edge, with the zone of influence extending to approximately 2H. Compared with traditional methods, the comprehensive stiffness coefficient proposed in this paper significantly enhances the accuracy of pile displacement prediction, which improves the fitting determination coefficient by 1.5 times. These findings offer valuable insights for optimizing the design and safety assessment of pile-braced foundation pits in loess areas.

Numerical investigations on cone penetration test based on MPDEM

WANG Junhao1, XU Wenjie1*, FU Jianbao2, LI Bin3
 2026, 45 (S1): 504-514 doi: 10.3724/1000-6915.jrme.2025.0334
[PDF] 5114 KB (54)     [HTML]
Show Abstract
The cone penetration test (CPT) is widely used in in-situ soil testing due to its operational simplicity and high reliability. To investigate the deformation and failure process of the soil and the interaction mechanism between the soil and probe during penetration, the advantages of the material point method (MPM) in simulating large deformation problems of geomaterials and the discrete element method (DEM) in contact calculations are fully combined in this study, which overcomes the limitations of traditional single numerical methods. Using the self-developed GPU-parallelized MPM-DEM coupling software CoSim-MPDEM, the entire process of CPT is analyzed. Based on Jardine?s laboratory tests, the contact parameters between the probe and soil are inverted, and their influence on cone tip resistance is explored. By comparing the variation of tip resistance with depth obtained from numerical and physical tests, the feasibility and advantages of MPDEM in numerical analysis of CPT are validated. Through the numerical tests, the deformation and failure mechanism of the surrounding soil during the penetration process is systematically investigated. The evolution of the displacement and stress fields inside the soil, as well as the distribution of contact forces and plastic zones are explored, which reveals the penetration mechanism of CPT. The results show that the contact forces between the probe and soil exhibit a non-uniform distribution, with significant stress concentration and plastic failure zones around the probe, where the plastic zones extend to 3 times the probe diameter. It is demonstrated that the proposed MPDEM algorithm achieves accurate contact calculations between large-deformation soil and structures, which provides valuable insights for soil-structure interaction analysis.

Wetting deformation and scale effect of rockfill materials based on the breakage packing concept

王 行,沈超敏*,马双星,刘斯宏,毛航宇
 2026, 45 (S1): 515-528 doi: 10.3724/1000-6915.jrme.2025.0407
[PDF] 4188 KB (59)     [HTML]
Show Abstract
The coupling mechanism between wetting-induced breakage and scale effects in rockfill materials remains insufficiently understood, and predictive methods for wetting deformation considering scale effects are still lacking. In this study, based on the breakage-packing concept, the volumetric deformation induced by wetting-related particle breakage is decomposed into two independent subproblems: the particle breakage evolution and the relationship between breakage and void ratio, providing a framework for unified consideration of wetting breakage and scale effects. Uniaxial compression wetting tests were conducted on rockfill specimens with different representative particle sizes under both single-line and double-line wetting paths. The results demonstrate that the breakage-packing concept is applicable to coupling analysis of wetting deformation and scale effects; after parameter normalization, the breakage index-plastic work relationships under varying particle sizes, moisture states, and wetting paths converge to a single curve, while the breakage index-packing void ratio relationships show invariance with respect to these factors. Finally, incorporating the Weibull scaling law for single-particle breakage, a mathematical volume-stress relationship model was developed to account for the coupled effects of wetting and scale. This model requires only one set of physically meaningful parameters, can be extended to prototype-scale materials, and provides a theoretical basis for predicting deformation of ultra-high rockfill dams during initial reservoir impoundment.
Copyright © 2005-2015 Edited and Published by Editorial Office for Chinese Journal of Rock Mechanics and Engineering
Sponsored by : Chinese Society for Rock Mechanics and Engineering 
Published by : Science Press 
Address: Wuhan 430071, China  Tel: (027)87199250  Fax: (027)87199250  Email: rock@whrsm.ac.cn
DESIGNED BY: Beijing Magtech Science & Technolgy Development Co.,Ltd.
鄂公网安备 42010602003581号