Discrete element modeling (DEM) is a useful tool for linking global responses of granular materials to underlying particle-level interactions. A DEM model capable of capturing realistic soil behavior must be calibrated to a reference dataset, typically consisting of laboratory experiments. Calibration of a DEM model often requires numerous simulations as contact parameters need to be iterated upon until the simulation results satisfactorily replicate the experimentally observed behaviors. This paper presents a sensitivity investigation that examines the effects of the contact parameters on the drained triaxial compression response of a poorly-graded sand. It then introduces a calibration procedure capable of providing contact parameters that satisfactorily reproduce the results of laboratory triaxial results in a few simulations. Results show that friction and rolling resistance coefficients jointly influence the mobilized peak and critical state friction angles, secant shear modulus, maximum dilation rate, total volumetric strain, and strain softening magnitude. These parameters also influence the mode of failure at contacts and the evolution of fabric anisotropy. The influence of mu r or mu on the triaxial response and particle-level interactions is coupled, becoming more profound as the other parameter is increased. Contact stiffness is shown to influence the shear modulus and volumetric change behavior independently of mu and mu r. An algorithm that estimates values for mu and mu r needed to reproduce experimental results is developed using triaxial response parameters from experimental datasets. The performance of the proposed calibration method is demonstrated for three natural sands showing that it provides appropriate calibrated parameters for poorly graded sands with different relative densities and confined with varying effective stress magnitudes.

This study introduces a novel discrete element method (DEM) model for compacted loess, incorporating a bond rate parameter within a linear contact bond model to simulate constitutive damage behavior. This enhancement significantly improves the characterization of structural damage from repeated wet-dry cycles, offering a quantitative method for predicting damage progression. Unlike existing DEM models, our model directly uses a bond rate parameter to quantitatively describe inter-particle bond deterioration, reflecting reduced bonding strength due to pore structure development and the weakening effect of water. Rigorous calibration and validation were performed using comparative experiments. A key innovation is the systematic analysis of microscopic parameters (contact stiffness, friction coefficient, contact strength, and bond rate) and their impact on macroscopic mechanical behavior. Our findings show that decreasing the bond rate significantly reduces the macroscopic mechanical properties, providing valuable insights into the micro-macro relationship. We comprehensively evaluated prediction sensitivity to these parameters. This methodology offers a new perspective on using DEM for predicting crucial civil engineering material properties, providing a valuable reference for incorporating bond rate parameters into future modeling, particularly for long-term geotechnical material behavior under environmental degradation. The model's accurate representation of wet-dry cycle effects on loess strength improves earth structure design and safety.

Understanding the response of sand to complex loading conditions is vital for practical geotechnical engineering. Circular rotational shear is a special stress path where the magnitudes of three principal stresses vary following a circular stress trajectory in the it-plane with their directions fixed. Although experimental studies under such stress paths are limited, the discrete element method appears to be an appealing approach to examine the response of granular materials to varying complex loading paths in numerical virtual tests. This study presents comprehensive numerical simulations of granular samples subjected to a circular stress path under varying conditions, including samples prepared with different bedding-plane angles and densities and subjected to different stress ratios. Both macroscopic and microscopic behaviors are presented and interpreted. A contactnormal-based fabric tensor is adopted in a detailed analysis to measure the internal structure of the granular assembly. The fabric, strain, and strain increment tensors are decomposed with respect to the stress tensor, and the evolutions of these components are presented along with the key influential factors. The results obtained in this study provide useful physical insight for the development of constitutive models for granular soils under general loading conditions.

Acoustic emission (AE) offers the potential to monitor and interpret soil-pipe interaction behavior by sensing particle-scale interactions. However, application of AE is limited by gaps in understanding related to how particle-scale interactions influence AE activity. Discrete element method (DEM) simulations of buried pipe uplift with energy tracking were performed and compared with experimental mechanical, displacement, and AE measurements, to ensure realistic behavior was captured by the modeling approach. A parametric investigation was then performed to evaluate the influence of pipe displacement direction and pipe diameter on plastic energy dissipation, and hence AE. Trends of dissipated plastic energy and measured AE with stress level (via burial depth) and pipe velocity were analogous. Relationships were quantified (R2 ranging from 0.74 to 0.98) between AE, dissipated plastic energy, and pipe velocity. Measured AE and dissipated plastic energy were linked with a general expression, comprising increments of friction (sliding and rolling), damping, and damage energies. Sliding friction energy accounted for >80% of the total dissipated energy on average during buried pipe deformation. Exemplar relationships were established between dissipated energy, pipe movement direction, embedment ratio, and mobilized soil volume (R2 values ranging from 0.92 to 0.97). A conceptual framework for interpreting buried pipe behavior using AE monitoring was presented.

This study investigates the influence of the soil-structure interaction (SSI) on the seismic performance of structures, focusing on the effects of foundation size, soil type, and superstructure height. While the importance of SSI is well recognized, its impact on structural behavior under seismic loads remains uncertain, particularly in terms of whether it reduces or amplifies structural demands. A simplified dynamic model, incorporating both the mechanical behavior of the soil and structural responses, is developed and validated to analyze these effects. Using a discrete element approach and the 1940 El Centro earthquake for validation, the study quantitatively compares the response of soil-interacting structures to those with fixed bases. The numerical results show that larger foundation blocks (20 m x 20 m and 30 m x 30 m) increase the seismic response values across all soil types, causing the structure to behave more like a fixed-base system. In contrast, reducing the foundation size to 10 m x 10 m increases the flexibility of structures, particularly buildings built on soft soils, which affects the displacement and acceleration response spectra. Softer soils also increase natural vibration periods and extend the plateau region in regard to spectral acceleration. This study further finds that foundation thickness has a minimal impact on spectral displacement, but structures on soft soils show more than a 15% reduction in spectral displacement (SD) compared to those on hard soils, indicating a dampening effect. Additionally, increasing the building height from 7 to 21 m results in a more than 20% decrease in SD for superstructures with natural vibration periods exceeding 2.4 s, while taller buildings with longer natural vibration periods exhibit opposite trends. Structures built on soft soils experience larger foundation-level displacements, absorbing more seismic energy and reducing earthquake accelerations, which mitigates structural damage. These results highlight the importance of considering SSI effects in seismic design scenarios to achieve more accurate performance predictions.

To achieve the loading of the stress path of hard rock, the spherical discrete element model (DEM) and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations. Furthermore, based on the deep tunnel of China Jinping Underground Laboratory II (CJPL-II), the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed, and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective. The results indicated that the stress-strain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress, accompanied by sudden changes in strain rate. Stress rotation induces spatially directional deformation, resulting in fractures of different degrees and orientations, and increasing the degree of deformation anisotropy. The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress, as well as its initial level is significant and positive. The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity. After transient unloading, both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize. This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on the epsilon 1-epsilon 3 1-epsilon 3 plane. Transient unloading will induce a tensile stress wave. The ability to induce fractures due to changes in principal stress magnitude, orientation and rotation paths gradually increases. The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude, which is determined by the magnitude and rotation of principal stress. A high tensile strain rate is more likely to induce fractures under low minimum principal stress. (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/).

The evaluation of triggering potential and induced deformations is crucial for liquefaction assessment of geosystems constructed on, or comprised of, coarse-grained soils. A series of cyclic constant volume direct simple shear (DSS) tests were simulated using the discrete element method (DEM) to investigate the influence of gradation on the pre- and post-liquefaction behavior of coarse-grained soils, including the resistance to initial liquefaction, accumulation of shear strains, and evolution of fabric. The DEM simulations were conducted on isotropically consolidated cubical specimens composed of non-spherical particles with a coefficient of uniformity (CU) C U ) range between 1.9 and 6.4. At similar relative densities, specimens with broader gradations yield lower liquefaction triggering resistance than poorly-graded specimens. However, at the same initial state parameter, the more broadly-graded specimens have higher liquefaction triggering resistance. Regardless of the parameter of state used for comparison, the specimens with broader gradation accumulate post-triggering shear strains at smaller rates than the poorly-graded specimens. The presence of coarser particles in the broadly-graded specimens promotes stable packing and enhanced interlocking by forming a substantially large number of contacts, which limits the particle movement and produces lower post-liquefaction strain accumulation. In contrast, the finer particles contribute minimally to the overall specimen stability. Additionally, wider gradations increase the anisotropy in the contact forces, with the strong forces aligning in the major principal stress direction and the weak forces providing resistance against buckling. These insights contribute to a better understanding of the pre- and post-liquefaction triggering behavior in granular soils and the implication for the design of resilient infrastructure built on, or comprised of, coarse-grained broadly-graded soils.

This study presents an investigation into the mechanical behavior of geogrid-reinforced ballast subjected to cyclic loading focusing on the macro- and micromechanical features of the geogrid-ballast interaction mechanism. Key areas of interest include the effects of geogrid placement depth, aperture size, and stiffness on the motion of ballast particles, formation of contact force chains, and energy dissipation. A three-dimensional discrete element model, calibrated with experimental data, simulates ballast box tests performed on 300-mm-thick ballast layers reinforced by geogrids placed at depths ranging from 50 to 250 mm below the tie. The findings reveal that geogrids located within the upper 150 mm of the ballast layer significantly reduce tie settlement by minimizing particle movement, creating well-connected soil structures, and decreasing energy dissipation. Upon identifying 150 mm as the optimal geogrid placement depth, a parametric study evaluates the impact of the geogrid aperture size (A) and stiffness on the behavior of geogrid-reinforced ballast. The geogrid aperture size (A) is varied to give aperture size to ballast diameter (D) ratios ranging from 1.09 to 2.91, while the geogrid's stiffness ranges from 9.54 to 18.00 kN/m. Results indicate that A/D ratios greater than or equal to 1.45 are required for geogrids to perform satisfactorily, while stiffness appears to wield a negligible influence on the response of geogrid-reinforced ballast.

The current study adopts a micromechanical approach to explore the nature of stress transmission in wet granular materials. First, we derive the discrete form of the capillary stress tensor obtained from homogenization to show the virial nature of capillarity through the application of point -wise capillary forces in Discrete Element Modeling (DEM). Furthermore, the non -spherical character of the capillary stress tensor is highlighted through a series of DEM triaxial simulations. Contrary to common thinking, the capillary stress tensor has indeed both mean and deviatoric components due to the underlying micromechanical aspects. Relevant key dimensionless parameters are identified to evaluate the relative magnitude of the capillary stress to the externally applied and contact (intergranular) stresses, thus determining the specific conditions under which the contribution of the deviatoric part becomes considerable. In addition, a DEM simulation of a simple shear test is performed to confirm the anisotropy (non -sphericity) of capillary stress tensor. Finally, the effective nature of the contact stress in the sense of Terzaghi for the constitutive behavior of wet granular materials is investigated via a DEM stress probing analysis. Results suggest that a single contact stress variable - germane to an effective stress - cannot relate to strain for the constitutive law in triphasic condition.

Reinforcing calcareous sands with geogrids is a potentially effective method for large-scale geotechnical constructions in coastal lands. The breakable nature of polygonal calcareous sands determines the complex particlegeogrid interactions. A three-dimensional numerical model of geogrid reinforced calcareous sand (GRCS) was established to investigate the potential mechanical laws based on the discrete element method (DEM), and the reasonableness of the numerical model was verified by comparing with the indoor triaxial test. It follows that the macro-microscopic mechanical behavior of GRCS under the influence of aperture size and tensile resistance of geogrids was further investigated via effective DEM simulations. The presented results show that the decreased aperture size and increased tensile resistance are beneficial to enhance the macro-mechanical properties of GRCS, including strength, internal friction angle and pseudo cohesion. Particle crushing is mainly affected by shear strain and confining pressure. The bulging deformation of GRCS is partially suppressed due to the confining effect of geogrids. Besides, the source of strength enhancement of GRCS is revealed based on the microscopic particlegeogrid interactions, and the calculation method of horizontal and vertical additional stresses in the reinforced soil element considering the effects of tensile resistance and aperture size is further established.