Earthquakes are common geological disasters, and slopes under seismic loading can trigger coseismic landslides, while also becoming unstable due to accumulated damage caused by the seismic activity. Reinforced soil slopes are widely used as seismic-resistant geotechnical systems. However, traditional geosynthetics cannot sense internal damage in reinforced soil systems, and existing in-situ distributed monitoring technologies are not suitable for seismic conditions, thus limiting accurate post-earthquake stability assessments of slopes. This study presents, for the first time, the use of a batch molding process to fabricate self-sensing piezoelectric geogrids (SPGG) for distributed monitoring of soil behavior under seismic conditions. The SPGG's reinforcement and damage sensing abilities were verified through model experiments. Results show that SPGG significantly enhances soil seismic resistance and can detect soil failure locations through voltage distortions. Additionally, the tensile deformation of the reinforcement material can be quantified with sub-centimeter precision by tracking impedance changes, enabling high-precision distributed monitoring of reinforced soil under seismic conditions. Notably, when integrated with wireless transmission technology, the SPGG-based monitoring system offers a promising solution for real-time monitoring and early warning in road infrastructure, where rapid detection and response to seismic hazards are critical for mitigating catastrophic outcomes.

A set of direct shear tests on the soil-geotextile interface (SGI) were conducted using a temperature-controlled constant normal stiffness (CNS) direct shear apparatus. This was done in order to evaluate the effects of normal stiffness, initial normal stress, soil water content, and temperature on SGI shear behavior and microdeformation patterns. The observations indicate that all shear stress-shear displacement curves demonstrate strain-hardening characteristics, with SGI cohesion and friction angle increasing at higher normal stiffness and lower temperatures. At freezing conditions, water content significantly affects the interface friction angle, while this effect is minimal at positive temperatures. Normal stress increases with higher water content, lower temperatures, and higher normal stiffness. Shear stress initially rises with normal stress before decreases, with a more pronounced rise under sub-zero conditions. Normal stress shrinkage shows a positive correlation with normal stiffness. Micro-deformation analysis of soil particles at the interface indicates significant strain localization within the shear band, which is less pronounced under sub-zero temperatures compared to positive temperatures. These patterns of normal displacement vary across analysis points within the shear band, with the macroscopic normal displacement reflecting a cumulative effect of these microscopic variations.

The majority of existing effective stress-based constitutive models approach thermal effects through the temperature dependency of surface tension and its effects on the soil-water retention curve (SWRC) and effective stress. Experimental tests and theoretical studies, however, suggest that the temperature effect on surface tension alone is not sufficient to properly explain thermal-induced changes in the effective stress and SWRC. This study focuses on the temperature-dependent elastoplastic behavior of low plasticity unsaturated soils by developing a set of constitutive-level relations that incorporate temperature-dependent SWRC and effective stress models. These models account for the effect of temperature on the enthalpy, contact angle, and surface tension. The application of the presented constitutive relations was demonstrated and validated for low plasticity soils, specifically incorporating temperature effects into the hardening modulus, specific volume change, yield stress of the modified Cam-Clay model, and stress-strain relationships. The proposed relationships are incorporated in any effective stress-based constitutive model for modeling temperature dependency of elastoplastic response in low plasticity unsaturated soils. Employing these relationships can enhance the numerical simulation of low plasticity unsaturated soils under thermo-mechanical or other coupled processes involving temperature-dependent conditions.

Localized rock failures, like cracks or shear bands, demand specific attention in modeling for solids and structures. This is due to the uncertainty of conventional continuum-based mechanical models when localized inelastic deformation has emerged. In such scenarios, as macroscopic inelastic reactions are primarily influenced by deformation and microstructural alterations within the localized area, internal variables that signify these microstructural changes should be established within this zone. Thus, localized deformation characteristics of rocks are studied here by the preset angle shear experiment. A method based on shear displacement and shear stress differences is proposed to identify the compaction, yielding, and residual points for enhancing the model's effectiveness and minimizing subjective influences. Next, a mechanical model for the localized shear band is depicted as an elasto-plastic model outlining the stress-displacement relation across both sides of the shear band. Incorporating damage theory and an elasto-plastic model, a proposed damage model is introduced to replicate shear stressdisplacement responses and localized damage evolution in intact rocks experiencing shear failure. Subsequently, a novel nonlinear mathematical model based on modified logistic growth theory is proposed for depicting the shear band's damage evolution pattern. Thereafter, an innovative damage model is proposed to effectively encompass diverse rock material behaviors, including elasticity, plasticity, and softening behaviors. Ultimately, the effects of the preset angles, temperature, normal stresses and the residual shear strength are carefully discussed. This discovery enhances rock research in the proposed damage model, particularly regarding shear failure mode. (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/

This paper investigates the mechanical response of coral sand under particle breakage using a hierarchical multiscale model combining the discrete element method (DEM) and the finite element method (FEM). This DEM-FEM model links the microscopic interaction mechanisms to macroscopic phenomena such as strain localization and failure. A cohesive contact model was first utilized to simulate compaction bands in the DEM and construct a cohesive assembly with smaller particles distributed around a larger particle to better simulate the grinding and angular breakage of coral sand. A representative volume element (RVE) that includes particle breakage was then constructed and analyzed under periodic boundary conditions. DEM analysis was performed, and the results were compared with triaxial compression test data obtained from the literature, demonstrating that the constructed RVE effectively represents the mechanical properties of coral sand. The constructed RVE was used for hierarchical multiscale simulations, which showed good agreement with existing triaxial testing of coral sand. Finally, by setting a larger cohesive force, the constructed coral sand particles were prevented from breakage, and comparative analysis revealed that particle breakage weakens the mechanical properties of coral sand. Furthermore, different shapes of coral sand particles were constructed, and RVE and hierarchical multiscale simulations of triaxial tests were performed. The results indicated that the triaxial tests of long strip-shaped coral sand particles exhibit higher peak values compared to spherical coral sand particles. Additionally, a double porosity model of coral sand was constructed to analyze the impact of internal porosity on soil mechanical properties. The results showed that the presence of internal porosity significantly weakened the mechanical properties of coral sand. These findings highlight the significant impact of particle breakage and shape on the mechanical behavior of coral sand, offering important insights for engineering applications.

This paper deeply couples the exponential-type nonlinear strain softening with the anisotropic method of microstructure tensor combined stress invariants, proposing an effective strength formula that reflects the anisotropy evolution of soil. Furthermore, an expression for the anisotropy ratio k of strength as an equivalent plastic strain-related variable is derived. For natural clay, this evolution of strength anisotropy is incorporated into the Mohr-Coulomb-matched Drucker-Prager (MC-matched DP) yield criterion within the Cosserat continuum framework, resulting in a more refined soil constitutive model. The main strength parameters required for this model can be conveniently obtained based on conventional soil tests, and the model functionality can be degraded through parameter adjustments. The detailed procedure of stress updating algorithm and the elastoplastic tangent modulus matrix are provided for the constitutive integration. Through the finite element implementation, the superiority of the model is demonstrated compared with existing literature. Also, a biaxial compression example is systematically analyzed to prove that the model can effectively reflect the sensitivity of soil to loading direction. Moreover, the evolution of the shear band morphology, particle rotation in the shear band, and the anisotropy degree presented by the model are consistent with previous experimental studies and discrete element method (DEM)-related literature results. Furthermore, the proposed model effectively addresses numerical convergence issues and mesh size dependence usually encountered in classical models during the simulation of strain localization occurred in the soil.

Accelerating creep before catastrophic failure commonly follows a power-law velocity-acceleration relationship, with the exponent typically near 2 but often evolving from 1 to 2 at a certain point, indicating a dynamic transition. The underlying mechanisms, however, remain unclear. Here we investigate this transition by monitoring the slip displacement of clayey soil during fluid-injection creep experiments. This transition is discontinuous in the first run but becomes continuous in the initially pre-sheared sample. Using a regularized rate-and-state friction model, we explicitly examine the relationship between the exponent and the frictional properties of the soil. This model describes the dynamic transition, with the exponent evolving from 1 to 2 across a broad range of frictional parameters. Furthermore, by incorporating idealized shear localization processes, the model qualitatively reproduces the shear-history-dependent transition. Our study demonstrates that a combination of structural evolutions and frictional properties may explain slow and fast slips observed in various shear systems.

The current understanding of the geotechnical behavior of lunar in-situ resources, particularly lunar regolith (LR), is significantly limited due to its scarcity. To address this gap, this research utilized the morphological characteristics of LR particles obtained from the Chang'E-5 (CE-5) mission to construct numerical simulants using the discrete element method (DEM). This approach was then employed to investigate the mechanical properties of LR. Firstly, high-definition lunar particle images from the CE-5 mission were selected to capture the morphological characteristics and grain size distribution. These morphological characteristics were linked with the rolling resistance parameter and incorporated into the three-dimensional (3D) micromechanical contact model. Additionally, a flexible boundary condition was employed in the triaxial simulation to ensure the evolution of strain localization. The relative particle translation gradient (RPTG) concept was utilized to capture the onset and development of strain localization during the shear process. The results indicated that the numerical lunar simulants can effectively reproduce the mechanical response of LR. Furthermore, at the particle scale, particle shape characteristics play a crucial role in particle rotation and translation during the shear process. This study may establish a foundation for lunar resource exploration and utilization techniques.

The stability performance of the frozen curtain formed under standpipe freezing is closely associated with the weak zone penetrated by thermal gradient-related fracture (TGF). The TGF-rich zone further affects the liquid phase flow when the frozen curtain is thawed. However, there is a lack of studies on the TGF-rich zone within the frozen curtain. To address this gap, a simplified and practical 2D bonded particle model-based numerical simulation strategy was developed to identify the possibility of acquiring field characteristics of the TGF-rich zone by conducting numerical tests on samples considering size effects. The results, validated by the experiment, indicated that the influence of size on crack localization zone was comparable to that of the parameter gradient but had a weaker characteristic on crack orientation, which represents the orientation of TGF. In particular, the characterization result of the TGF-rich zone using crack localization zone in the simulation closely matched that using lateral strain localization zone both in simulation and experiment. Regarding the size effects of the TGF-rich zone revealed in the simulation, the estimated field length of the TGF-rich zone accounted for approximately 30% of the zone width characterized by a horizontal thermal gradient, with maximum orthotropic deformation occurring at about 10% of the zone width. These observations validate the existence of TGF within the frozen curtain and contribute to the development of a precise grouting technique to mitigate subsidence within soil deposits subjected to freeze-thaw.

Shear strain localization refers to the phenomenon of accumulation of material deformation in narrow slip zones. Many materials exhibit strain localization under different spatial and temporal scales, particularly rocks, metals, soils, and concrete. In the Earth's crust, irreversible deformation can occur in brittle as well as in ductile regimes. Modeling of shear zones is essential in the geodynamic framework. Numerical modeling of strain localization remains challenging due to the non-linearity and multi-scale nature of the problem. We develop a numerical approach based on graphical processing units (GPU) to resolve the strain localization in two and three dimensions of a (visco)-hypoelastic-perfectly plastic medium. Our approach allows modeling both the compressible and incompressible visco-elasto-plastic flows. In contrast to symmetric shear bands frequently observed in the literature, we demonstrate that using sufficiently small strain or strain rate increments, a non-symmetric strain localization pattern is resolved in two- and three-dimensions, highlighting the importance of high spatial and temporal resolution. We show that elasto-plastic and visco-plastic models yield similar strain localization patterns for material properties relevant to applications in geodynamics. We achieve fast computations using three-dimensional high-resolution models involving more than 1.3 billion degrees of freedom. We propose a new physics-based approach explaining spontaneous stress drops in a deforming medium. Strain localization is the accumulation of strain in narrow regions of rocks and other materials like metals, soils, and concrete, occurring at different scales. The strength of most geomaterials, particularly rocks, is strongly pressure-dependent, with strength increasing with increasing pressure. We developed efficient numerical algorithms using High-Performance Computing (HPC) and graphical processing units (GPUs) to model strain localization in 2D and 3D for applications in geodynamics and earthquake physics. Unlike previous models, our method reveals non-symmetrical patterns by using very small strain increments, highlighting the need for high-detail modeling. We found that elasto-plastic and visco-plastic models show similar strain patterns for relevant materials. Our method also achieves fast, detailed computations with over 1.3 billion variables and offers a new explanation for sudden stress drops in deforming materials. We resolve material instability during deformation resulting in a non-symmetric pattern of strain localization We demonstrate the similarity in patterns of strain localization between frictional and time-dependent plasticity models We achieve fast numerical simulations in high-resolution model setups in three dimensions involving more than 500 million degrees of freedom