This study presents a hierarchical multiscale approach that combines the finite-element method (FEM) and the discrete-element method (DEM) to investigate tunneling-induced ground responses in coarse-grained soils. The approach considers both particle-scale physical characteristics and engineering-scale boundary value problems (BVPs) simultaneously, accurately reproducing typical tunneling-induced mechanical responses in coarsegrained soils, including soil arching and ground movement characteristics observed in laboratory tests and engineering practice. The study also unveils particle-scale mechanisms responsible for the evolution of soil arching through the underlying DEM-based RVEs. The results show that the rearrangement of microstructures and the deflection of strong contact force chains drive the rotation of macroscopic principal stress and the formation of soil arch. The microscopic fabric anisotropy direction can serve as a quantitative indicator for characterizing soil arching zones. Moreover, the effects of particle size distributions (PSD) and soil densities on ground deformation patterns are interpreted based on the stress-strain responses and contact network characteristics of DEM RVEs. These multiscale insights enrich the knowledge of tunneling-induced ground responses and the same approach can be applied to other geotechnical engineering analyses in coarse-grained soils.

Knowledge Gap: The aggregation of clay minerals-layered silicate nanoparticles-strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 mu m. Simulations investigated the impact of simulation time, platelet diameters (6 to 25 nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

The cumulative plastic deformation and damage evolution of frozen soil-rock mixtures under cyclic loading was studied by a dynamic triaxial instrument with real-time resistivity measurement function. A series of low- temperature cyclic triaxial tests were conducted under varying confining pressures (200 kPa, 500 kPa, 800 kPa), block proportions (0, 30 %, 40 %, 50 %), and dynamic stress ratios (0.4, 0.6, 0.8). The results reveal that the cumulative plastic deformation process can be divided into three stages, such as microcrack closure as the initial stage, crack steady growth as the middle stage, and rapid crack propagation until it fails as the final stage. Under the same number of cycles, the greater the dynamic stress is, the greater the cumulative plastic deformation is. Furthermore, a strong correlation is identified between the resistivity and the cumulative plastic deformation. With the increase of the number of cycles, the cumulative plastic deformation leads to the accumulation of internal damage, and the resistivity gradually increases. Thus, a damage evolution model based on resistivity damage variables is proposed. The model demonstrated an average fitting accuracy of 97.36 % with the experimental data.

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

In cold-region high-speed railway (HSR) subgrade engineering, coarse-grained soils are commonly used as frost heave prevention fillers. However, coupled water-heat migration during freeze-thaw cycles still induces frost heave. This study innovatively employs a nuclear magnetic resonance (NMR) system to elucidate the hydro-thermal transport mechanisms in coarse-grained soils during freezing. The results reveal that under identical temperature and freezing duration, high-water-content soils release substantial latent heat from pore water freezing, resulting in higher freezing zone temperatures than low-water-content soils. During freezing, unfrozen water content decreases as a power function with freezing time at different depths of soil samples, with the frozen zone experiencing the fastest water reduction, followed by the freezing front and then the unfrozen zone. Both free and bound water progressively decrease in frozen and unfrozen zones. After freeze-thaw, the change in soil pore structure leads to a decrease in bound water and an increase in free water in frozen zones, while both decrease in unfrozen zones. Furthermore, higher initial water content results in more pronounced reductions of bound water and increases of free water in frozen zones. These findings advance the understanding of hydro-thermal coupling mechanisms and provide theoretical foundations for frost damage mitigation in high-speed railway subgrades.

The moisture accumulation and freezing damage of coarse-grained fill (CGF) in high-speed railway (HSR) subgrades have been widely concerned. Based on the newly developed water-vapor-heat-mechanical coupling test apparatus, a series of soil column tests were carried out to investigate the frost heave mechanism of CGF. The results indicate that the liquid water in CGF is discontinuous and difficult to migrate to the freezing front. The primary mechanism of moisture accumulation and frost heave in CGF is vapor migration and phase transition. With increasing freeze-thaw cycles, both vapor migration and frost heave reduce. The thaw settlement of the CGF is less than the frost heave, so there is a net upward deformation in each cycle. Furthermore, the fine particle content has a prominent effect on the heat transfer and frost heave of the CGF compared to the fine particle type. Even under the condition of vapor replenishment, controlling the content of fine particles is still an important way to inhibit frost heave. Moreover, after reducing the maximum particle size of CGF, the frost heave of the sample increases. Nuclear magnetic resonance (NMR) test results show that CGF is dominated by large pores, and the freeze-thaw cycle further promote the development of large pores, providing a good channel for the migration of vapor. In conclusion, the frost heave development caused by vapor migration is slow and continuous, posing a non-negligible risk to HSR subgrades during long-term service.

Coarse-grained soil is generally used in cold-regions infrastructure to mitigate the frost damage to engineering because of its non-frost heave susceptibility; however, in certain cases, coarse-grained fill has been observed to experience frost heave under hydraulic pressure. To reveal the mechanism of hydraulic pressure on coarsegrained soil frost heave, a model was developed to describe the frost heave in coarse-grained soil, incorporating the migration of external water to ice lenses through an unfrozen water film under hydraulic pressure, then the model was validated using published results. Subsequently, based on the validated model, the influence mechanism of hydraulic pressure and fine content on coarse-grained soil frost heave were analyzed. The calculation results demonstrate that the hydraulic pressure aggravates frost heave by increasing the pore water pressure gradient in the unfrozen water film. Additionally, frost heave rate increases with fine content because of the thickening of the film, which facilitates water flow and ice segregation. Furthermore, gray correlation analysis demonstrated that the impact of hydraulic pressure on frost heave in coarse-grained soil is more significant than that of fine content. Finally, the study discusses frost damage that occurred in high-speed railway subgrade and proposes the preventive measures.

This study investigated the hydraulic and mechanical behaviors of unsaturated coarse-grained railway embankment fill materials (CREFMs) using a novel unsaturated large-scale triaxial apparatus equipped with the axis translation technique (ATT). Comprehensive soil-water retention and constant-suction triaxial compression tests were conducted to evaluate the effects of initial void ratio, matric suction, and confining pressure on the properties of CREFMs. Key findings reveal a primary suction range of 0-100 kPa characterized by hysteresis, which intensifies with decreasing density. Notably, the air entry value and residual suction are influenced by void ratio, with higher void ratios leading to decreased air entry values and residual suctions, underscoring the critical role of void ratio in hydraulic behavior. Additionally, the critical state line (CSL) in the bi-logarithmic space of void ratio and mean effective stress shifts towards higher void ratios with increasing matric suction, significantly affecting dilatancy and critical states. Furthermore, the study demonstrated that the mobilized friction angle and modulus properties depend on confining pressure and matric suction. A novel modified dilatancy equation was proposed, which enhances the predictability of CREFMs' responses under variable loading, particularly at high stress ratios defined by the deviatoric stress over the mean effective stress. This research advances the understanding of CREFMs' performance, especially under fluctuating environmental conditions that alter suction levels. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

To investigate the effects of freeze-thaw (FT) cycles on the mechanical properties of coarse-grained soil in southeastern Xizang under different moisture contents, this study focuses on coarse-grained soil from a large landslide deposit in Linzhi City, Xizang. FT cycle tests, triaxial shear tests, and numerical simulations were employed to systematically examine the comprehensive impact of varying FT cycles, moisture content, and confining pressure on the soil's mechanical characteristics. The results show that FT cycles significantly affect the stress-strain behavior of coarse-grained soil in southeastern Xizang. The degree of strain softening increased from approximately 11.6% initially to 31.2% after 15 FT cycles, with shear strength decreasing by an average of 31.8%. Specifically, cohesion decreased by 38% to 55% after 0 to 15 FT cycles, and the internal friction angle decreased by approximately 29% to 32%. Additionally, higher moisture content led to more pronounced strain softening and strength degradation, while increased confining pressure effectively mitigated these deteriorative effects. Numerical simulation results indicated that as moisture content increased from 7.6% to 11.6%, the number of FT cycles required to reach the critical instability state decreased from approximately 150 to 106, and finally to only 15, with the maximum equivalent plastic strain increasing from 0.20 to 2.47. The findings of this study provide key mechanical parameters for understanding the formation and evolution of FT landslide disasters in southeastern Xizang and lay a scientific foundation for the assessment and long-term prevention of cold-region geological hazards.

Take the reservoir landslide as an example, in addition to hydrological conditions, creep properties of soils play an important role in explaining the mechanisms behind landslide movement. Although the change of this deformation over time is small, the long-term accumulation will also bring new hidden danger to the safety control of the slope. This paper takes the shallow coarse-grained soils of Qiaotoubei landslide as the research interest, improves the test method for the deficiency of not allowing the lateral deformation of the specimen in the traditional one-dimensional compression creep test, and conducts the compression creep test of coarse-grained soils by using the modified high-pressure consolidation instrument. Based on this test data, the creep property of coarse-grained soils is analyzed and a suitable creep constitutive model is selected, that is generalized Kelvin model. Then, relevant parameters are determined and FLAC3D software is used to simulate the creep deformation of the slope deposits and the stress and deformation of the lattice beams. Finally, the coupling mechanism between coarse-grained soils creep and lattice structure was analyzed based on the comparison of the calculated results with the deformation or damage in the field. Through this study, some targeted suggestions and directions for future research are proposed for the management of reservoir deposit landslides, hoping to contribute to the operational safety of the reservoir.