The discrete element method (DEM) has demonstrated significant advantages in simulating soil-tool interaction and an appropriate contact model notable affected the simulation accuracy. The accuracy of numerical simulation is compromised due to the variations in soil properties when tillage implements are employed in clay-moist soil conditions. This study aims to establish a discrete element model of clay-moist soil based on the Edinburgh Elasto-Plastic Adhesion (EEPA) contact model. Calibration tests using a combination of direct shear tests and cone penetration tests were conducted to identify sensitive parameters that need to be calibrated in the model and analyze the effects of each parameter. The results indicated that contact plasticity ratio and surface energy had significant influence on representing the mechanical properties of clay-moist soil. Then, by utilizing scanning technology to acquire furrow shape data, soil bin test was conducted to validate the reliability of the calibration parameters. Using sensitive parameters as variables, the actual value of clay-moist soil with a moisture content of 33 % as the target value obtained from experimental tests. The optimal combination was: the coefficient of static friction of 0.45, the coefficient of rolling friction of 0.18, and the surface energy of 27.95 J.m-2, the contact plasticity ratio of 0.59. The relative error between the simulated draft force value and the actual measured value was 7.98 %, and the relative errors in the furrow type parameters did not exceed 5 %. The accuracy of the calibration results was verified through comparative analysis of simulation and empirical results. This study provides a scientific approach for employing DEM in modeling clay-moist soil-tool interaction.

Microbial Induced Calcium Carbonate Precipitation (MICP), recognized as a low-carbon and environmentally sustainable consolidation technique, faces challenges related to inhomogeneous consolidation. To mitigate this issue, this study introduces activated carbon into uranium tailings. The porous structure and adsorption capacity of activated carbon enhance bacterial retention time, increase the solidification rate, and promote the growth and distribution of calcium carbonate, resulting in more uniform consolidation and improved mechanical properties of the tailings. Additionally, a novel independently developed grouting method significantly enhances the mechanical strength of the tailing sand samples. To perform a micro-scale analysis of the samples, distinct activated carbon-tailings DEM models are constructed based on varying activated carbon dosages. Physical experiments and parameter calibration are employed to investigate the micro-mechanical properties, such as velocity field and force chain distribution. Experimental and simulation results demonstrate that incorporating activated carbon increases the calcium carbonate production during the MICP process. As the activated carbon content increases, the peak stress of the tailings initially rises and then declines, reaching its maximum at 1.5 % activated carbon content. At 100 kPa confining pressure, the peak stress is 2976.91 kPa, 1.23-1.59 times that of samples without activated carbon and 6.08-7.86 times that of unconsolidated samples. Micro-scale motion analysis reveals that particle movement is predominantly axial at the ends and radial near the central axis. The initial direction of the primary force chains aligns with the loading direction. Following failure, some primary force chains dissipate, while new chains form, predominantly along the axial direction and secondarily in the horizontal direction. Compared with samples without activated carbon, those containing activated carbon exhibit more uniform force chain distribution, higher stress levels, and greater peak stress. This study offers a novel approach to enhance the stabilization and solidification efficiency of MICP and establishes a DEM model that provides valuable insights into the structural deformation and micro-mechanical characteristics of MICPcemented materials.

In cold regions, the strength and deformation characteristics of frozen soil change over time, displaying different mechanical properties than those of conventional soils. This often results in issues such as ground settlement and deformation. To analyze the rheological characteristics of frozen soil in cold regions, this study conducted triaxial creep tests under various creep deviatoric stresses and established a corresponding Discrete Element Method (DEM) model to examine the micromechanical properties during the creep process of frozen clay. Additionally, the Burgers creep constitutive model was used to theoretically validate the creep deformation test curves. The research findings indicated that frozen clay primarily exhibited attenuated creep behavior. Under low confining pressure and relatively high creep deviatoric stress, non-attenuated creep was more likely to occur. The theoretical model demonstrated good fitting performance, indicating that the Burgers model could effectively describe and predict the creep deformation characteristics of frozen clay. Through discrete element numerical simulations, it was observed that with the increase in axial displacement, particle displacement mainly occurs at both ends of the specimen. Additionally, with the increase in creep deviatoric stress, the specimen exhibits different deformation characteristics, transitioning from volumetric contraction to expansion. At the same time, the vertical contact force chains gradually increase, the trend of particle sliding becomes more pronounced, and internal damage in the specimen progresses from the ends toward the middle.

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.

The delayed breakage of particles significantly affects the long-term mechanical properties of rockfill materials. This study examines the effects of particle strength dispersion on the distribution of time-dependent strength using fracture mechanics and probabilistic methods. Subsequently, the distribution of normalized maximum contact force (NMCF), defined as the ratio of the maximum contact force to instantaneous strength, for specimens with uniform particle size is derived using extreme value theory and Discrete Element Method (DEM). Based on this analysis, the probabilities of delayed breakage in rockfill specimens over various time intervals are calculated using a joint probability delayed breakage criterion. The feasibility of the proposed method is validated by comparing theoretical calculation with DEM triaxial creep simulation results that accounted for particle breakage. The findings offer innovative tools and theoretical insights for understanding and predicting the particle delayed breakage behavior of rockfill materials and for developing macro-micro creep crushing constitutive models.

The macroscopic mechanical properties of granular systems largely depend on the complex mechanical responses of force chains at the mesoscopic level. This study offers an alternative to rapidly identify and predict force chain distributions under different stress states. 100 sets of gradation curves that effectively represent four typical continuous gradation distributions are constructed. Numerical specimens corresponding to these gradation curves are generated using the discrete element method (DEM), and a dataset for deep neural network training is established via biaxial compression numerical simulations. The relationship between particle distribution characteristics and force chain structure is captured by the Pix2Pix conditional generative adversarial network (cGAN). The effectiveness of the generated force chain images in reproducing both particle gradation and spatial distribution characteristics is verified through the extraction and analysis of pixel probability distributions across different color channels, along with the computation of texture feature metrics. In addition, a GoogLeNet-based prediction model is constructed to demonstrate the accuracy with which the generated force chain images characterize the macroscopic mechanical properties of granular assemblies. The results indicate that the Pix2Pix network effectively predicts and identifies force chain distributions at peak stress for different gradation

The mechanical behaviour of soil subject to shear loading or deformation is typically considered either completely drained or undrained. Under certain conditions, these drained and undrained scenarios can represent boundaries on the allowed volumetric strain. There is growing interest in exploring the response under intermediate conditions where partial drainage is allowed, particularly in the development of new approaches to mitigate the risk of liquefaction induced failure and the design of off-shore structures. This study uses the discrete element method (DEM) to investigate the effect of partial drainage conditions on the mechanical behaviour of spherical assemblies. Samples with different interparticle friction values are isotropically compressed and then subjected to undrained, drained, and partially drained triaxial shearing. The partially drained conditions are simulated in the DEM samples by applying a controlled volumetric strain that is a fraction of the drained volumetric strain. Results on loose samples indicate that allowing drainage enhances peak shear resistance and can also prevent liquefaction. Moreover, dense samples show a substantial increase in shear resistance when small changes in drainage and volumetric strain take place. The peak stress ratio and the stress ratio at the phase transformation point are insensitive to the drainage level. There is a linear correlation between the state parameter and the drainage level at the peak stress ratio and the phase transformation point. This observation could be used to trace partially drained stress-paths and could also aid the development of uncoupled constitutive models that account for drainage effects.

The physical and mechanical characteristics of saline soil are significantly influenced by salt content, with macro- and mesoscopic mechanical properties closely correlated. This study investigates the strength and deformation behaviors of saturated saline sand through indoor triaxial shear testing under varying confining pressures and salt contents. The key innovation lies in developing a coupled finite element and discrete element analysis model to simulate the mesoscopic behavior of saline sand under triaxial shear stress state. Flexible boundary conditions were applied, and appropriate contact models for salt-sand interactions were selected. By adjusting mesoscopic parameters, stress-strain curves and variations in porosity, coordination number, particle displacement, and contact force chains were analyzed. The study further explores shear band development and shear failure mechanisms by examining relative particle displacement and the breaking of contact force chains. Additionally, the influence of salt particle size on the overall strength of the DEM model was assessed. The findings provide valuable insights into the internal structural changes of saline sand during shear deformation, contributing to a better understanding of its mechanical behavior in engineering applications.

This paper analyzed the influence of the inherent anisotropy of sand on active and passive arching by simulating the trapdoor emplying the discrete element method (DEM). The inherent anisotropy is reflected by the bedding plane angle alpha of particles. The granular material constitutive responses are captured on representative volume elements (RVEs). A new modeling method is employed to prepare particle specimens, aiming to obtain a more uniform soil model. The results indicate that the discrete element method can simulate the influence of the inherent anisotropy of granular material on the evolution of soil arching. An asymmetric arching evolution phenomena is observed in the alpha other than 0 degrees or 90 degrees cases, which leads to obvious asymmetric deformation and stress distribution in the soil. As the filling height increases, this phenomenon becomes more and more obvious. From a microscopic perspective, the reorientation of the contact normal fabric caused by particle rotation is the main reason for the differences in soil arching evolution with different alpha. This study provides a theoretical basis for predicting ground deformation failure caused by underground engineering activities and changes in surrounding environmental conditions.

Construction and Demolition Wastes (CDW) serves as an effective filler for highway subgrades, demonstrating commendable performance characteristics. The efficient utilization of CDW not only contributes to environmental sustainability but also yields significant economic benefits. This study employs discrete element simulation to develop a triaxial sample model comprising particles with four distinct levels of sphericity. By varying the combinations of sphericity, brickconcrete ratio, and void ratio, triaxial simulation tests are conducted, and the critical state soil mechanics framework is applied to fit the critical state line (CSL) of the samples. The results indicate that sphericity, brick-concrete ratio, and void ratio substantially influence the macroscopic mechanical properties of CDW. Notably, as sphericity increases, the peak deviatoric stress of the samples decreases, and significant volume deformation occurs. The slope of the CSL in the q-p ' plane diminishes, while the slopes of both forms of the CSL in the e-log p ' plane increase. Furthermore, a decrease in the brick-concrete ratio enhances the anti-deformation and compressive capacities of the samples. As the brick-concrete ratio decreases, both the slopes and intercepts of the CSL in the e-log p ' plane exhibit an upward trend. Conversely, an increase in the void ratio leads to a reduction in the overall strength and anti-deformation capacity of the specimens, an increase in the compressibility of the specimen volume, an elevation of the CSL slope on the q-p ' plane, and a gradual increase in both the slope and intercept of the semilogarithmic form of the CSL on the e-log p ' plane, as well as a gradual increase in the slope of the power-law form of the CSL.