The service performance of frozen soil is one of the important factors that needs to be considered in designing and assessing the safety of artificial ground freezing projects. We conducted shear tests on ice-containing frozen soil and assessed soil performance and damage characteristics of the ice-frozen soil interface. On the basis of experimental results, we further investigated the damage of ice-containing frozen soil numerically using the finite-discrete element method. Experimental and numerical results show that temperature, the normal load, and moisture content are the primary factors influencing the mechanical properties of the ice-frozen soil interface. The effects of these parameters on shear strength, shear modulus, cohesion, and angle of internal friction were analyzed and discussed. There was a transition from ductile to brittle behavior at the ice-frozen soil interface with decreasing temperature. Transition occurred at higher temperatures in soils with higher moisture content. Because ice and sand differ in terms of stiffness, fractures appeared first at the ice-frozen sand interface. Under continued loading, the specific form of damage and maximum load-bearing capacity varied as a function of the location of the maximum shear stress zone and the ice in the soil. Our research findings provide valuable theoretical insights for the design and evaluation of the safety of artificial ground freezing engineering projects.

The mechanical properties and failure characteristics of soil-rock mixtures (SRMs) directly affect the stability of tunnels constructed in SRMs. A new SRM modelling method based on the combined finite-discrete element method (FDEM) was proposed. Using this new SRM modelling method based on the FDEM, the mechanical characteristics and failure behaviour of SRM samples under uniaxial compression, as well as the failure mechanism of SRMs around a tunnel, were further investigated. The study results support the following findings: (1) The modelling of SRM samples can be achieved using a heterogeneous rock modelling method based on the Weibull distribution. By adjusting the relevant parameters, such as the soil-rock boundaries, element sizes and modelling control points, SRM models with different rock contents and morphologies can be obtained. (2) The simulation results of uniaxial compression tests of SRM samples with different element sizes and morphologies validate the reliability and robustness of the new modelling method. In addition, with increasing rock content, VBP (volumetric block proportion), the uniaxial compressive strength and Young's modulus increase exponentially, but the samples all undergo single shear failure within the soil or along the soil-rock interfaces, and the shear failure angles are all close to the theoretical values. (3) Tunnels in SRMs with different rock contents all exhibit X-shaped conjugate shear failure, but the fracture network propagation depth, the maximum displacement around the tunnel, and the failure degree of the tunnel in the SRM roughly decreases via a power function as the rock content increases. In addition, as the rock content increases, such as when VBP = 40%, large rocks have a significant blocking effect on fracture propagation, resulting in an asymmetric fracture network around the tunnel. (4) The comparisons of uniaxial compression and tunnel excavation simulation results with previous theoretical results, laboratory test results, and numerical simulation results verify the correctness of the new modelling method proposed in this paper.

Purpose - The purpose of this paper is to propose a new combined finite-discrete element method (FDEM) to analyze the mechanical properties, failure behavior and slope stability of soil rock mixtures (SRM), in which the rocks within the SRM model have shape randomness, size randomness and spatial distribution randomness. Design/methodology/approach - Based on the modeling method of heterogeneous rocks, the SRM numerical model can be built and by adjusting the boundary between soil and rock, an SRM numerical model with any rock content can be obtained. The reliability and robustness of the new modeling method can be verified by uniaxial compression simulation. In addition, this paper investigates the effects of rock topology, rock content, slope height and slope inclination on the stability of SRM slopes. Findings - Investigations of the influences of rock content, slope height and slope inclination of SRM slopes showed that the slope height had little effect on the failure mode. The influences of rock content and slope inclination on the slope failure mode were significant. With increasing rock content and slope dip angle, SRM slopes gradually transitioned from a single shear failure mode to a multi-shear fracture failure mode, and shear fractures showed irregular and bifurcated characteristics in which the cut-off values of rock content and slope inclination were 20% and 80 degrees, respectively. Originality/value - This paper proposed a new modeling method for SRMs based on FDEM, with rocks having random shapes, sizes and spatial distributions.

During the operational phase of pumped storage hydropower stations, rockfill materials within the dam experience cyclic loading and unloading due to water level fluctuations. This cyclic behavior can result in the accumulation of irreversible deformation, posing a substantial threat to dam safety. However, there is an absence of a constitutive model capable of accurately capturing the low-frequency multi-cycle hysteresis behavior of rockfill materials due to the constraints of conventional laboratory test methods. In this study, we employed the combined finite and discrete element method (FDEM) to investigate the mechanical characteristics of rockfill materials and develop an improved constitutive model capable of effectively capturing their hysteretic behavior. The results demonstrate the FDEM accurately reproduces the mechanical behavior of rockfill material under shear and cyclic loading and unloading conditions. The hysteresis loop exhibits a discernible densification trend with increasing cycles. And the variations in elastic modulus and strain primarily occur within the initial five cycles. The plastic strain increment exhibits a strong positive correlation with stress level, while its relationship with confining pressure is comparatively less pronounced. The proposed constitutive model successfully captures the complex low-frequency multi-cycle hysteresis characteristics of rockfill material with few parameters, showing substantial potential for practical applications.

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering. The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles. Direct shear tests were conducted on cubic slate samples with foliation angles of 0 degrees, 30 degrees, 45 degrees, 60 degrees, and 90 degrees. The effect of foliation angles on failure patterns, acoustic emission (AE) characteristics, and shear strength parameters was analyzed. Based on AE characteristics, the slate failure process could be divided into four stages: quiet period, step-like increasing period, dramatic increasing period, and remission period. A new empirical expression of cohesion for layered rock was proposed, which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments. The comparative analysis demonstrated that the new expression has better prediction ability than other expressions. The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method (FDEM), and it was found to align well with the experimental results. Considering both computational efficiency and accuracy, it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations. To balance the relationship between the number of elements and the simulation results in the direct shear simulations, the recommended element size is 1 mm. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

In this study, the axial swelling strain of red-bed mudstone under different vertical stresses are measured by swell-under-load method, and the microstructure of mudstone after hygroscopic swelling is studied by mercury intrusion porosimetry (MIP). The weakening coefficient and Weibull distribution function are introduced into the coupling model of mudstone moisture diffusion-swelling deformation-fracture based on finite-discrete element method (FDEM). The weakening effect of moisture on mudstone's mechanical parameters, as well as the heterogeneity of swelling deformation and stress distribution, is considered. The microcrack behavior and energy evolution of mudstone during hygroscopic swelling deformation under different vertical stresses are studied. The results show that the axial swelling strain of mudstone decreases with increase of the vertical stress. At low vertical stresses, moisture absorption in mudstone leads to formation of cracks caused by hydration-induced expansion. Under high vertical stresses, a muddy sealing zone forms on the mudstone surface, preventing further water infiltration. The simulation results of mudstone swelling deformation also demonstrate that it involves both swelling of the mudstone matrix and swelling caused by crack expansion. Notably, crack expansion plays a dominant role in mudstone swelling. With increasing vertical stress, the cracks in mudstone change from tensile cracks to shear cracks, resulting in a significant reduction in the total number of cracks. While the evolution of mudstone kinetic energy shows similarities under different vertical stresses, the evolution of strain energy varies significantly due to the presence of different types of cracks in the mudstone. The findings provide a theoretical basis for understanding the hygroscopic swelling deformation mechanism of red-bed mudstone at various depths. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

Anthropogenic climate change threatens water storage and supply in the periglacial critical zone. Rock glaciers are widely distributed alpine aquifers with slower response to temperature increases, that provide the summer water flow of many alpine streams. Knowing the extent and makeup of rock glaciers is necessary to evaluate their potential for water supply. We used non-invasive methods, integrating geological, geomorphological, meteoro-logical, and geophysical information to characterize the internal structure and hydrology of the Upper Camp Bird rock glacier (UCBRG) located on level 3 of Camp Bird Mine in Ouray, Colorado, and assessed the applicability of two electromagnetic induction systems in this highly heterogeneous landform with a history of anthropogenic activity. The time-domain (G-TEMTM) system achieved deep subsurface penetration (similar to 100 m) and realistic modeling of the internal structure of the UCBRG: a shell of volcanic rock fragments (< 3 m thick; 1-100 Ohm-m), a meltwater component (10(2)-10(3) Ohm-m), located between 50 and 100 m near the toe (subpermafrost flow), and 1-30 m in the soundings farthest from the toe (suprapermafrost flow within the active layer), and a frozen component (permafrost 50-80 m thick; 10(3)-10(6) Ohm-m). The frequency-domain system, however, was highly susceptible to local environmental conditions, including anthropogenic objects (i.e., mine carts, lamp posts, tunnel tracks, etc.) and was unable to resolve UCBRG's internal makeup. The non-invasive methodology and general conceptual framework presented here can be used to characterize other alpine aquifers, contributing to the quantification of global water resources, and highlighting the importance of preserving rock glaciers as storage for critical water supply in the future.

在气候变暖的背景下，冰川运动逐渐活跃，发生冰崩灾害的可能性和危险性随之增加，研究冰川失稳过程与机理具有重要的现实意义。基于有关文献中的人工制备多晶冰室内静力三轴试验结果，采用FDEM方法建立多晶冰试样的数值模型标定参数、分析破坏机理，进行了不同条件下西藏阿里地区阿汝冰川失稳过程的二维数值模拟。通过FDEM方法计算冰体从连续到离散的变形特征，描述出阿汝冰川的从稳态开始、裂缝的产生及发展、局部滑移、崩解和碰撞反弹等破坏过程。同时，发现冰川和基岩的连接界面摩擦减小比冰川断裂增加对冰川失稳的影响更大。FDEM方法可以应用于冰川失稳破坏研究，在研究非连续问题领域也有相当大的潜力。