In geotechnical engineering, the development of efficient and accurate constitutive models for granular soils is crucial. The micromechanical models have gained much attention for their capacity to account for particle-scale interactions and fabric anisotropy, while requiring far less computational resources compared to discrete element method. Various micromechanical models have been proposed in the literature, but none of them have been conclusively shown to agree with the critical state theory given theoretical proof, despite the authors described that their models approximately reach the critical state. This paper modifies the previous CHY micromechanical model that is compatible with the critical state theory based on the assumption that the microscopic force-dilatancy relationship should align with the macroscopic stress-dilatancy relationship. Moreover, under the framework of the CHY model, the fabric anisotropy can be easily considered and the anisotropic critical state can be achieved with the introduction of the fabric evolution law. The model is calibrated using drained and undrained triaxial experiments and the results show that the model reliably replicates the mechanical behaviors of granular materials under both drained and undrained conditions. The compatibility of the model with the critical state theory is verified at both macroscopic and microscopic scales.

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.

A two-lift gradient design for airport pavements has been proposed to mitigate the functional degradation, especially the salt-frost (S-F) damage induced by deicing slat fluids. Herein, this study focuses on elucidating the mechanism and improvement of incorporating mineral admixtures in the development of a novel S-F resistant surface concrete material, which is of great significance for delaying the functional deterioration of pavement surface in northern China. The results indicated that the filling effect and secondary hydration reaction between the fly ash (FA) and silica fume (SF) and cement hydration products results in a dense spatial network structure, effectively reducing porosity and optimizing pore structure. It was found that SF can effectively improve the frost resistance and salt corrosion resistance of cement mortar, while the influence of FA depends on its content and environmental conditions. The incorporation of FA and SF significantly enhanced the structural density of cement concrete and reduced chloride ion permeability. The improvement in impermeability is most pronounced when both FA and SF are used in combination. In addition, a fitting equation between the admixture content and chloride ion permeability has been established, demonstrating good fitting results. In non-frozen saline soil areas, a large amount of FA or SF could be incorporated; in seasonally frozen areas, the priority should be given to SF to ensure salt corrosion resistance and frost resistance. The findings of this study provide a scientific basis for sustainable airport pavement construction in northern China.

The large amount of slag generated during the construction of earth pressure balance shield (EPBS) not only incurs significant disposal costs, but also exacerbates environmental pollution. To improve the utilization of the shield slag, silty clay with additive is proposed as a slag conditioner instead of bentonite. Firstly, various macroscopic properties of the bentonite and silty clay slurries are tested. Subsequently, the relationships between the macroscopic properties of the silty clay slurries containing additives and the modification mechanism are evaluated at microscopic, mesoscopic, and macroscopic scales by using infrared spectroscopy (IR), scanning electron microscope (SEM), and Zeta potential tests, respectively. Based on these tests, reasons for variations in modification effects of different slurries are identified. The results show that addition of 3 % sodium carbonate to the silty clay can effectively improve the rheological properties of the slurry. The modification mechanism of sodium carbonate involves the formation of hydrogen bonds between water molecules and inner surface hydroxyl groups within the lattice layer of kaolinite. This process significantly enhances the rheological properties of the silty clay slurry. Furthermore, sodium carbonate alters the contact relationships between the silty clay particles, which increases viscosity and reduces permeability of the slurry. Finally, sodium carbonate increases thickness of the electrical double layer of the silty clay particles. This allows the particles to bind more water molecules, therefore improving slurry-making capacity of the silty clay. This paper presents an innovative multiscale analysis of the modification process of silty clay. The substitution of recycled silty clay for bentonite as a slag conditioner not only substantially reduces the cost of purchasing materials, but also considerably decreases the expenses associated with transportation and disposal of the soil discharged by EPBS.

Due to its particulate nature, the mechanical properties of bulk clay are determined by interparticle forces and fabrics of particle assemblies. A thorough study of the connection between properties across length scales is crucial to a fundamental understanding of the mechanisms behind the complex mechanical behavior of clays and clayey soils. This paper demonstrates the development of a multiscale constitutive model for describing the small-strain elastic properties of illite, based on the results of coarse-grained mesoscale molecular dynamic simulations for monodisperse assemblies of illite primary particles. The formulation consists of a homogenization scheme linking the potential energy of the system with an optimal parameter set describing the mesoscale fabric of the particles, and a perturbation scheme describing the change of the parameters in response to infinitesimal strains applied to the systems. The small strain elastic stiffness tensors are calculated as the second-order derivative of the potential energy with respect to the infinitesimal strain. The results from model prediction are validated against the stiffness properties interpreted from numerical simulations as well as experimental findings from prior research studies. The multiscale constitutive model is able to effectively capture the elastic properties of illite in terms of magnitude and material symmetry purely based on the information of interparticle forces and fabrics.

This paper introduces a thermo-hydro-mechanical (THM) framework to model thaw consolidation in permafrost regions. By integrating internal energy degradation functions and a modified Cam-Clay model within a phase-field damage framework, the model focuses on simulating the simultaneous effects of phase change and particle rearrangement. The model integrates two distinct phase-field variables with the modified Cam-Clay plasticity framework. One phase-field variable monitors pore phase composition, while the other captures particle rearrangement. These variables are directly coupled to the constitutive model, providing critical data for updating the stress-strain relationship by accounting for particle rearrangement-induced softening and hardening effects due to volumetric deformation. The model converges to the modified Cam-Clay model when there is no phase change. This approach addresses a significant gap in existing models by capturing the associated microstructural evolution and plastic softening in thaw-sensitive soils. Validation efforts focus on experimental scenarios assessing both the mechanical impacts of thaw consolidation and the dynamics of phase transitions, particularly emphasizing latent heat effects. The results demonstrate the proposing model's capability of handling complex behaviors of permafrost under thaw conditions, confirming its potential for enhancing infrastructure resilience in cold regions.

Frozen mixed soils are widely distributed in the strata and slopes of permafrost regions. This paper aims to study the strength criterion and elastoplastic constitutive model for frozen mixed soils from micro to macro scales. Based on the knowledge of mathematical set theory and limit analysis theory, the support function of frozen soils matrix is derived. The concept of local equivalent strain is proposed to solve the problem of nonuniform deformation caused by rigid inclusions in frozen mixed soils. According to the nonlinear homogenization theory and the Mori-Tanaka method in micromechanics, the strength criterion of frozen mixed soils is established, which can consider coarse particle contents. By introducing the concepts of equivalent yield stress and equivalent plastic deformation, the elastoplastic constitutive model is proposed by the associated flow rule, which can also consider the influence of coarse particle contents. Finally, using the data in the literature, the proposed strength criterion and elastoplastic constitutive model for frozen mixed soil are verified, respectively. The effects of coarse particle contents on the mechanical properties of frozen mixed soils are discussed.

Solid particles may experience different kinds of cohesive forces, which cause them to form agglomerates and affect their flow in multiphase systems. When such systems are simulated through computational fluid dynamics (CFD) programs, appropriate modelling tools must be included to reproduce this feature. In this review, these strategies are addressed for various systems and scales. After an introduction of the different forces (van der Waals, electrostatic, liquid bridge forces, etc.), the modelling approaches are categorized under three methodologies. For diluted slurries of very fine particles, many researchers succeeded with pseudo-single phase approaches, employing a model for the non-Newtonian rheology. This was especially popular for sludges in anaerobic digestions or certain types of soils. In other cases, continuum-based approaches seem to be more adequate, including cohesiveness in the kinetic theory of granular flows or the restitution coefficient. Geldart-A particles experiencing van der Waals forces are the primary focus of such studies. Finally, when each particle is modelled as a discrete element, the cohesive force can be directly specified; this is especially widespread for the wet fluidization case. For each of these approaches, a general overview of the main strategies, achievements, and limits is provided.

Sandy cobble soil is a composite made of soil matrix and cobbles, and the estimation of its shear strength always requires expensive large-scale experiments. The strength of the sandy cobble soil exhibits macroscopic anisotropy with respect to the direction of the major stress due to the observed dominant distribution of the cobble dip angle. In the present paper, a numerical homogenization procedure for anisotropic strength identification of the sandy cobble soils is established, which can take into account the influencing factors of the size, shape, and inclination of the cobbles and the mesoscopic strength of the soil-rock interface. To consider the condition of plain strain, the particle size distribution of the cross of the stratum is derived based on the fractal theory and the transformation method of Walraven. The mesostructure of the sandy cobble soils is randomly produced using ellipses to model the cross of the cobbles. An iterative procedure is utilized to represent the major stress orientation-dependent macroscopic strengths. The results are validated against the data from indoor experiments and global mesoscopic computations. It is shown that the macroscopic strength of the sandy cobble mixtures can be accurately determined and the iterative multiscale limit analysis method is reliable and efficient. Parameter analysis is finally conducted to discuss the effect of the mesoscopic properties on the macroscopic strength.

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.